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ABSTRACT 


The  objective  of  the  STAR  project  is  to  test  and  space 
qualify  a  new  continuous  cycle  cryogenic  refrigeration 
system  for  cooling  of  sensors  and  electronics  which  is 
based  upon  the  newly  discovered  thermoacoustic  heat  pumping 
effect.  The  new  refrigerator  has  no  sliding  seals,  a  cycle 
frequency  of  about  300  Hz,  and  uses  acoustic  resonance  to 
enhance  the  overall  power  density  and  efficiency.  This 
thesis  is  concerned  specifically  with  the  design  and 
testing  of  the  thermodynamic  element  (or  stack),  which  is 
responsible  for  the  thermo-acoustic  power  conversion,  and 
the  testing  of  binary  inert  gas  mixtures  as  working  fluids. 
Using  the  refrigerator’s  coefficient  of  performance 
relative  to  the  ideal  Carnot  coefficient  of  performance  as 
a  measure  of  efficiency,  we  have  achieved  a  93%  improvement 
over  previous  designs. 
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I.  INTRODUCTION 


A.  BACKGROUND 

This  being  the  second  m  a  series  of  theses,  covering 
different  aspects  of  the  Space  Thersoacoust ic  Refrigerator 
(STAR)  project,  a  great  deal  of  the  introductory  aatenal 
has  already  bee*,  written.  In  order  for  this  thesis  to  be 

complete,  the  introductory  material  that  is  felt  to  be 
necessary  will  be  repeated  here.  With  the  exception  of 
minor  corrections.  Chapter  I  and  sections  A. 1  and  A. 2  of 
Chapter  II  have  been  reproduced  from  the  work  of 
Fitzpatrick  (1988). 
i.  History 

Thermoacoustics  can  generally  be  described  as  the 
study  of  the  interaction  between  heat  and  sound.  Scientific 
interest  m  this  interaction  is  not  new.  Lord  Rayleigh 
(1878  and  1945)  discussed  various  qualitative  examples  of 
the  production  of  sound  by  heat.  In  one  of  these  examples 
he  describes  investigations  by  Sondhauss  ( 1850)  of  an 
effect  observed  by  glassblowers ?  heating  a  bulb  of  glass  at 
the  en  •  of  a  hollow  rod  occasionally  led  to  the  generation 
of  sound.  Faraday  (1818)  demonstrated  that  this  effect 
occurred  with  gases  other  than  air.  Sondhauss  proved  that 
the  vibration  of  the  glass  itself  did  not  generate  the 
sound,  but  he  offered  no  explanation  as  to  what  did.  In  his 
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description  of  Sondhauss*  work,  Rayleigh  stressed  the 
importance  of  the  phase  difference  between  temperature  and 
particle  motion.  Rijkm  (1859)  built  an  open  tube  with  a 
wire  mesh  inside,  When  the  mesh  was  heated,  the  tube 
produced  a  sound  of  impressive  intensity.  The  functioning 
of  this  oscillator  is  slightly  different  since  it  requires 
steady  gas  flow  for  its  operation.  A  similar  sound 
production  effect  was  observed  by  Tacoms,  et  al.  (1949)  in 
hollow  tubes  immersed  in  liquid  helium  (4.2*K).  The  Tacoms 
oscillations  were  addressed  quantitatively  by  Yazaki,  et 
al.  (1989).  The  work  described  above  deals  primarily  with 
prime  movers.  These  are  devices  that  convert  a  temperature 
gradient  into  sound  energy,  be  call  this  the  classical  era 
of  thermoacoustics. 

Interest  in  thernoacoust ics  was  renewed  when  the 
idea  of  the  use  of  acoustical  devices  as  refrigerators 
appeared.  Gifford  and  Longsworth  (1966)  described  the 
pumping  of  heat  along  a  surface  caused  by  a  periodic  change 
in  the  pressure  of  the  adjoining  gas.  Such  a  change  can  be 
produced  by  the  oscillations  of  a  sound  wave.  In  their 
experiments  Herkli  and  Thonann  (1974)  explored  the  heating 
and  cooling  effects  on  the  wall  of  a  gas-filled  resonant 
tube.  They  found  that  heat  was  transported  from  a  region 
near  the  velocity  antinode  (or  maximum  value)  of  the 
sound’s  standing  wave  to  regions  near  the  adjacent  pressure 
antinodes.  Thermoacoustic  prime  movers  have  also  received 
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recent  Attention.  Kempton  (1976)  invest igated  the  excess 
noise  of  aeroengines  above  that  predicted  by  theory.  He 
deteramed  that  the  sound  uas  produced  by  unsteady  heat 
transfer.  Each  of  these  latter  three  croups  of 
experimenters  used  some  theory  for  comparison,  but  it  uas 
mostly  qualitative.  Hone  provided  the  complete  quantitative 
theory  that  would  explain  their  experimental  observations. 

The  theoretical  breakthrouch  in  the  thermoacoustics 
field  uas  made  by  Nicklaus  Rott  (1969,  1974,  1975  and 
1989).  His  theory  combined  the  basic  principles  from 
physics,  thermodynamics  and  acoustics  to  quantitatively 
describe  the  effects  found  in  both  types  of  thermoacoustic 
devices:  prime  movers  and  heat  pumps.  Prime  movers  use  a 
temperature  gradient  to  create  sound,  as  in  the  effects  dis¬ 
cussed  by  Sondhauss,  Rayleigh,  Taconis  and  Kempton.  Heat 
pumps,  or  refrigerators,  use  the  oscillating  pressure  of  a 
sound  uave  to  produce  a  temperature  difference,  as 
described  by  Gifford  and  Longsuorth  and  Herkli  and  Thomann. 
Rott \  described  the  effect  found  by  Merkli  and  Thomann  as 
thermoacoustic  streaming.  It  is  this  effect  that  makes 
thermoacoustic  refrigeration  possible. 

Inspired  by  Ceperley’s  (1979)  traveling  uave  Ster¬ 
ling  Cycle  heat  engine  and  Rott’s  quantitative  theory,  the 
team  of  Wheatley,  Hofler,  Swift,  Migliori  and  Garrett 
(1982,  1983a,  1983b,  1985  and  1986)  developed  a  series  of 
thermoacoust ic  experiments  at  Los  Alamos  National 
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Laboratory  in  Non  Haxico.  Thty  investigated  the  basic 
thermoacoustic  effect*  in  both  prime  movers  and 
refrigerators  and  compared  thair  experimental  results  to 
Kott’s  theory. 

2.  XhlliS&IUUU £& 

We'll  digress  from  history  here  to  explain  the 
thermodynamic  distinction  between  pnae  aovers  and 
refrigerators.  This  discussion  follows  Sears  and  Salinger 
(1975).  A  prime  mover  receives  heat  at  a  high  temperature, 
Tjj  (the  hot  reservoir),  does  work  on  its  surroundings,  and 
rejects  heat  at  a  lower  temperature,  Tg  (the  cold 
reservoir),  as  shown  in  Figure  I-l.a.  The  first  law  of  ther¬ 
modynamics  tells  us  that  the  heat  flow  (Q^)  from  the  hot 
reservoir  must  equal  the  work  done  on  the  surroundings  ill) 
plus  the  heat  flow  (Qg)  to  the  cold  reservoir: 

QH=W+Qg.  (1-1) 

The  second  law  of  thermodynamics  states  that  in  every 
process  the  entropy  of  the  universe  must  either  remain 
constant  or  increase.  The  entropy  change  for  the  isothermal 
heat  transport  to/from  the  reservoirs  is  defined  as  the 
heat  flow  divided  by  the  reservoir  temperature.  For  a  prime 
mover  the  second  law  can  be  written: 

iSH<iSg  or  QH/THiQC/TC*  (1-2) 

The  thermal  efficiency  of  the  prime  mover  is  defined  as  the 
ratio  of  the  work  output  to  the  heat  input: 

*[  =U/Qh=  ( QH-Qg  )  /Qh=  1-  ( Qg/QH ) .  (  I  -3 ) 
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(a) 


Figure  1-1.  The  Thermodynamics  of  (a)  Prime  Movers 
and  (b)  Refrigerators. 


Comparing  equations  1-2  and  1-3  gives  us  a  limit  on  the 
efficiency  based  on  the  temperatures  of  the  reservoirs: 

il-(Tc/TH).  (1-4) 

This  limit  is  known  as  the  Carnot  efficiency. 

In  a  refrigerator  the  process  is  basically  reversed 
Uee  figure  I-l.b. ).  Heat  (Qc)  is  removed  from  the  low 


temperature 

reservoir, 

work 

(W) 

is  received 

from  the 

surroundings. 

and  heat 

•V 

IS 

delivered  to 

the  higher 

temperature  reservoir.  The  expression  for  the  first  law  is 
the  same  as  for  the  prime  mover.  Equation  1-1.  For  the 
refrigerator  it  means  the  heat  flow  into  the  hot  reservoir 
must  equal  the  heat  flow  from  the  cold  reservoir  plus  the 
work  done  on  the  system.  The  second  law  for  this  case  is: 

ASh>ASc  or  QH/TH>QC/TC.  (1-5) 

The  efficiency  parameter  for  the  refrigerator  is  called  the 
coefficient  of  performance  (COP).  It  is  defined  as  the 
ratio  of  the  heat  flow  from  the  cold  reservoir  to  the  work 
input  to  the  refrigerator,  or: 

CQP=QC/W=QC/ ( Qh-Qc ) <TC/ ( Th-Tc ) .  ( I -6 ) 

The  limit  is  known  as  the  Carnot  coefficient  of 
performance. 

The  work  in  this  thesis  deals  solely  with  refriger¬ 
ators.  Even  though  prime  movers  are  mentioned  from  time  to 
time,  the  primary  focus  of  the  rest  of  our  discussion  will 
be  refrigerators. 
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According  to  Hofler  (1986),  the  early  thermoacoustic 
work  at  Los  Alamos  focused  on  experimental  refrigerating 
engines.  The  performance  of  these  engines  fell  short  of 
expectations,  leading  to  simple  experiments  on  basic 
thermoacoustic  effects  and  a  proof-of-pr inc iple  experiment 
on  thermoacoustic  refrigeration.  Hofler  then  applied  the 
Rott  theory  to  the  experimental  systems  and  solved  the 
resulting  equations  numerically.  For  his  doctoral 
dissertation  from  the  University  of  California,  San  Diego, 
Hofler  designed  and  constructed  a  completely  functional 
thermoacoustic  refrigerator.  He  also  made  accurate 
measurements  of  its  thermodynamic  efficiency,  and  used  this 
efficiency  to  make  comparisons  to  the  Rott  theory.  After 
receiving  his  doctorate,  Hofler  came  to  the  Naval 
Postgraduate  School  (NPS)  as  a  post-doctoral  fellow,  and 
brought  his  prototype  refrigerator  with  him. 

The  purpose  of  this  thesis,  in  conjunction  with  the 
work  of  several  other  students,  is  to  modify  Hofler’s 
refrigerator  design  in  order  to  improve  its  overall 
efficiency  and  make  it  suitable  for  space  cryocooler 
applications.  The  basic  design  of  the  thermoacoustic 
refrigerator  is  shown  schematically  in  Figure  1-2.  The 
driver  (A),  which  produces  the  sound  waves,  is  coupled  to 
the  resonator  (D)  using  bolts,  via  a  reducer  cone  (B)  and 
bellows  (C).  Inside  the  resonator  is  a  stack  of  plastic 


plates  (£)  and  their  heat  exchangers  (F).  The  heat 
exchangers  allow  heat  to  be  removed  from  the  hot  end  of  the 
stack  and  applied  to  the  cold  end.  It  is  the  interaction 
between  the  sound  waves  and  the  plastic  plates  that 
produces  a  temperature  difference  across  the  plates  and 
allows  heat  to  be  pumped  from  the  cold  heat  exchanger  to 
the  hot  heat  exchanger,  A  brief  description  of  efficiency 
would  be  useful  here.  Figure  1-3  shows  a  diagram  of  the 
chain  of  energy  transformations  which  occur  in  a 
thermoacoustic  refrigerator.  There  are  three  conjugate 
variable  pairs  (V  and  I,  P  and  U,  and  AT  and  Qq)>  and  two 
transformations  (electr ical-to-acoust ical  and  acoust ical-to- 
thermal ) . 

We  start  with  electrical  input  power  to  the  driver, 
which  can  be  calculated  (using  r.m.s.  values)  with: 

=IVcose,  (1-7) 

e  iec 

where  I  is  the  current  into  the  electroacoustic  driver,  V 
is  the.  voltage  across  the  driver  terminals,  d  is  the  phase 
angle  between  I  and  V,  and  cos$  is  known  as  the  "power 
factor.  "  An  electrodynamic  driver  converts  the  electrical 
power  to  acoustic  power  (the  first  transformation). 
Analogous  to  the  previous  definition  of  the  electrical 
input  power,  the  acoustic  power  is  given  by: 

W=PUcos  (1-8) 

where  P  is  the  acoustic  pressure,  U  is  the  volumetric 
velocity,  which  equals  the  product  of  particle  velocity  and 
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Figure  1-3.  Chain  of  Energy  Transforiat ions 
in  a  Theraoacoustic  Refrigerator 


resonator  cross-sectional  Area?  and  *  is  the  phase  angle 
between  thee.  The  efficiency  of  this  electroacoustic  trans¬ 
formation  is  given  by: 


^EA“W/Pelec* 


(1-9) 


In  the  second  transformation,  the  refrigerator  con¬ 
verts  the  acoustic  power  to  a  temperature  gradient  (  T)  and 
heat  flow  (Q).  As  discussed  previously,  the  efficiency  for 
this  transition  is  given  by  the  coefficient  of  performance: 


COP=:Qc/W. 


(1-10) 


The  . overall  thermoacoustic  efficiency  is  therefore  the 
product  of  the  efficiencies  for  the  two  transformations: 


^TA^EA'  (C0P)a^C/Pelec' 


( 1-11) 


Improvements  in  the  overall  refrigerator  efficiency 
can  thus  be  divided  into  two  distinct,  though  interrelated, 
efforts.  The  first  is  the  subject  of  this  thesis: 
modifications  to  improve  the  coefficient  of  performance. 
The  second  effort  involves  modifications  to  the  driver  to 
optimize  the  electroacoustic  efficiency  and  is  addressed  by 
LT  Michele  Fitzpatrick,  U.S.  Coast  Guard,  in  her  master’s 
thesis. 


£.  SPACE  THEFMQ ACOUSTIC  REFRIGERATOR  (STAR) 

1 .  Motivation 

In  addition  to  improving  overall  efficiency,  our 
experimental  refrigerator  will  be  designed  and  built  with 
the  intention  of  flying  it  on  the  space  shuttle  as  part  of 

II 


the  National  Aeronautics  and  Space  Administration *s  (NASA) 
Get  Away  Special  (GAS)  program.  As  ot'  3  February  1988, 
there  is  a  Memorandum  of  Agreement  between  the  Naval  Post¬ 
graduate  School  (NPS)  and  the  Air  Force  which  provides 
funding  for  the  space  flight  and  which  assigns  the  NASA 
payload  number  G-337  to  this  project.  The  reason  for 
testing  the  thermoacoustic  refrigerator  in  space  is  that 
the  immediate  applications  of  the  STAR  are  expected  to  be 
the  cooling  of  electronics,  high  T^  superconductors  and 
infrared  detectors  in  space. 

There  are  currently  two  primary  methods  for  cooling 
satellite  systems  in  space:  (1)  evaporation  of  expendable 
cryogens  (liquid  helium,  nitrogen,  ammonia-methane  and 
solid  hydrogen,  etc.),  and  (2)  closed  cycle  refrigerators 
(Stirling  cycle,  Vuilleumier  cycle,  etc.),  which  involve 
large  reciprocating  piston  masses  operating  at  low 
frequency.  CUalker  (1983)  and  Smith,  et  al.  i  1984)3  The 
disadvantages  of  these  current  cooling  methods  are  their 
short  lifetimes  (expendable  cryogens)  and  high  vibration 
levels  and  low  reliability  (closed  cycle  refrigerators). 
The  advantages  of  the  thermoacoustic  refrigerator  over 
these  and  other  systems  lie  in  its  simplicity  and 
reliability.  The  STAR  has  no  sliding  seals,  relatively 
efficient  performance,  low  vibration  levels  and  (hopefully) 
a  longer  life  span. 
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The  thermoacoust ic  refrigerator  needs  to  be 
surrounded  by  a  vacuum  and  insulating  aaterial  to  minimize 
heat  loss  due  to  thereal  conduction  and  radiation.  The 
vacuum  of  space  will  provide  the  best  insulation  to 
determine  the  refrigerator's  ultimate  efficiency.  The 
absence  of  gravity  reaoves  the  possibility  of  thermal 
convection  as  an  additional  nuisance  heat  transport  effect. 

2.  Get  Away  Special  (GAS)  Program 

NASA's  GAS  program  allows  small,  self-contained  pay- 
loads  to  fly  on  the  Space  Shuttle  in  GAS  canisters  for 
relatively  low  cost  ($10,000)  (Get  Away  Special  Team, 
1984).  The  GAS  can  is  five  cubic  feet  m  volume  and  can 
house  a  payload  of  up  to  200  pounds.  Each  payload  must 
contain  its  own  electrical  power,  control,  data  acquisition 
and  storage  facilities,  etc.  The  Shuttle's  astronauts  will 
turn  a  switch  on  or  off  at  designated  times  during  the 
flight,  but  are  otherwise  not  involved  with  the  experiment. 

3.  Shared  Subsystems  from  NASA  G-313 

Another  group  of  NPS  students  and  faculty  (Boyd,  et 
al.,  1987)  have  taken  advantage  of  the  GAS  program  to 
erasure  the  resonant  acoustic  modes  of  the  shuttle  payload 
bay  and  the  ambient  acoustic  environment  produced  as  a 
result  of  main  engine  and  booster  operation  during  launch. 
Their  experiment  is  titled  "The  Space  Shuttle  Cargo  Bay 
Vibroacoust ics  Experiment"  and  is  designated  by  NASA  as 
payload  G-313.  Several  subsystems  that  were  developed  for 
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NASA  G-313  will  be  used  by  NASA  G-337  (STAX) .  A  schematic 
of  the  STAX  in  its  GAS  can  is  shown  in  Figure  1-4.  One  of 
the  systems  borrowed  from  NASA  G-313  is  the  controller 
system  used  to  run  the  experiment  and  record  the  data.  The 
recorder  system  consists  of  the  INTEL  model  BFK  5V75 
magnetic  bubble  momory  module  and  an  NSC  800  micro¬ 
processor-based  controller.  Two  other  NPS  students,  LT 
Charles  B.  Cameron,  USN  and  CPT  Ronald  Byrnes,  USA,  will  be 
designing  the  analog  electronics  and  software  to  integrate 
these  systems  into  the  STAR  experiment  for  their  master's 
theses.  A  block  diagram  of  the  electronics  is  shown  in 
Figure  1-5. 

Another  borrowed  system  is  the  power  supply,  which 
consists  of  Gates  brand  lead  acid  battery  cells  (five 
ampere-hour,  two  volts  each).  These  gelled  electrolyte 
batteries  are  ideal  for  the  STAR  due  to  their  high  power 
density,  low  cost  and  the  low  level  of  outgassing  during 
discharge  cycles.  NASA  G-313  used  a  one  layer  battery  of 
63  cells  providing  680  watt-hours  of  energy  and  weighing 
about  80  pounds  (including  the  cells*  support  structure). 
We  will  be  using  two  battery  layers  with  as  many  cells  as 
ue  can,  up  to  136.  This  will  give  us  a  minimum  of  680  watt- 
hours  and  a  maximum  of  1,360  watt-hours  of  available 
electrical  energy,  and  a  total  battery  weight  between  80 
and  160  pounds.  The  number  of  cells  we  can  use  will  depend 
on  the  total  weight  of  the  driver-resonator  assembly  and 
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Figure  1-5.  Block  Diagraa  of  the  Electronics 
for  the  Space  Theraoacoust lc  Refrigerator 
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its  avxilliary  equipment  (vacuum  can,  gaa  reservoir,  ate.) 
and  alactronics,  as  we  have  to  stay  within  tha  299  pound 
1 lait. 

Tha  usa  of  tha  GAS  can  laposas  certain  restric¬ 
tions.  Since  wa  are  using  batteries  to  supply  the  power, 
the  refrigerator  has  to  be  energy  efficient.  Also,  the  GAS 
can  setup  requires  tha  STAS  to  ba  compact  and  lightweight. 
Figure  1-6  shows  a  photograph  of  Hofler’s  prototype 
refrigerator.  This  setup,  which  includes  parts  that  will  ba 
unnecessary  in  space,  is  approximately  six  feat  in  height. 
In  comparison,  the  maximum  payload  height  for  tha  GAS  can 
is  28.25  inches,  or  lass  than  2.5  feet.  These 
considerations  played  a  major  role  on  the  choice  of 
equipment  for  and  the  design  of  the  STAS, 

C .  SCOPE 

This  thesis  is  primarily  experimental  in  nature.  The 
purpose  being  to  improve  the  thermodynamic  efficiency  of 
Hofler’s  prototype  refrigerator.  Theory  is  presented  only 
as  necessary  to  provide  the  reader  with  an  understanding  of 
the  thermoacoust ic  refrigeration  process  and  to  explain  why 
certain  modifications  were  made.  For  this  reason,  no  direct 
comparison  of  overall  refrigerator  performance  to  theory  is 
presented.  Hofler  (1986)  presents  a  quantitative  comparison 
of  the  refrigerator’s  performance  to  the  detailed  Rott 
theory. 
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Figure  1-6.  Photograph  of  Hofler's 
Prototype  Refrigerator 
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Chapter  II.  begins  with  a  fairly  in  depth  discussion  of 
thermoacoustic  theory,  to  provide  the  above  mentioned 
understanding  and  to  illustrate  possible  areas  for  the 
improvement  of  efficiency.  This  is  followed  by  the  details 
of  how  and  why  specific  modifications,  such  as  gas  mixtures 
and  stack  designs,  were  made. 

Chapter  III.  provides  all  of  the  information  necessary 
for  the  construction  of  a  thermoacoustic  stack. 

Chapter  IV.  provides  a  discussion  of  the  overall 
experimental  apparatus  and  procedures.  We  have  basically 
used  the  prototype  system  of  Hofler  (1986)  for  our 
experiments.  The  material  in  sections  A  and  B  of  this 
chapter  are  derived  primarily  from  Hofler’s  doctoral 
dissertation  (1986),  with  emphasis  on  the  feu  areas  that  we 
have  changed. 

Chapter  V,  presents  the  methods  used  for  analysis  of  our 
data  and  the  overall  results  of  our  experiments. 

Chapter  VI.  provides  conclusions  and  recommendations  for 
further  development. 
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II.  THEORY 


A.  THERMOACOUSTIC  THEORY 

As  stated  in  Chapter  I.,  the  first  two  sections  (A.  1. 
and  A. 2)  of  this  Chapter  are  essentially  reproduced  from 
the  master’s  thesis  of  LT  Michele  Fitzpatrick  (1988). 

Thermoacoustic  theory  has  been  developed  in  detail  by 
Rott  (  1969,  1974,  1975  and  1980)  and  adapted  to  the  thermo¬ 
acoustic  refrigerator  by  Wheatley  et  al.  (1982,  1983a, 
1983b,  1985  and  1986),  Wheatley  and  Cox  (1985),  Hofler 
(1986)  and  Suift  (1989). 

The  space  thermoacoustic  refrigerator  (STAR)  basically 
consists  of  an  acoustic  driver  producing  sound  uaves  in  a 
resonant  tube  (see  Figure  1-2).  This  resonator  is  filled 
with  a  mixture  of  helium-xenon  gas  (*12.5  atomic  X  xenon) 
pressurized  to  ten  atmospheres,  and  contains  a  stack  of 
plastic  plates.  This  section  of  the  Chapter  will  present  a 
qualitative  model  for  the  thermoacoustic  heat  pumping 
process  followed  by  a  quantitative  development  for  heat  and 
work  flow  at  a  single  plate.  This  is  followed  by  a  more 
rigorous  development  for  the  heat  flow,  work  flow  and 
coefficient  of  performance  for  a  stack  of  plates  that  are 
much  shorter  than  one  quarter  wavelength.  Section  A. 
concludes  with  a  discussion  of  some  of  the  changes  that 
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could  improve  the  coefficient  of  performance  of  the 
refrigerator. 

1 .  Qualitative  Picture 

Consider  a  parcel  of  gas  that  moves  back  and  forth 
along  one  of  the  plates  at  the  acoustic  frequency  (see 
Figure  II-l).  As  it  moves,  the  jp-as  parcel  will  experience 
changes  in  temperature  and  volume.  Part  of  the  temperature 
changes  come  from  the  adiabatic  compression  and  expansion 
of  the  gas  by  the  sound  pressure,  and  part  as  a  consequence 
of  the  local  temperature  of  the  plate  itself.  A  temperature 
gradient  may  develop  along  the  plate  as  a  result  of  the 
operation  of  the  refrigerator.  The  temperature  and  volume 
changes  can  be  described  by  six  separate  steps  (the  changes 
in  volume  are  indicated  by  changes  in  the  size  of  the 
square  parcel  in  the  figure). 

Assume  the  plate  is  at  a  uniform  temperature,  T.  The 
parcel  starts  at  position  1  (X=0)  with  a  temperature  of  T. 
The  sound  wave  moves  the  parcel  to  the  right  to  position  2. 
The  parcel  has  undergone  adiabatic  compression  and  its 

+  4. 

temperature  is  now  T  .  Since  the  temperature  of  the  gas  is 
higher  than  that  of  the  plate,  heat  (Q)  will  flow  from  the 
gas  parcel  to  the  plate.  The  parcel’s  volume  is  decreased 
and  its  temperature  lowered  to  T+.  This  heat  flow  also 
causes  the  plate’s  temperature  to  increase  at  the  position 
X+.  Parcel  position  3  is  actually  in  the  same  location  as 
position  2,  but  is  displaced  vertically  in  the  figure  for 
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PLATE 


Figure  I ! - 1 


Gas  Parcel  Diagram.  (Note:  Displacement 
In  the  Vertical  Direction  Is  Shown  for 
C 1 ar 1 ty- -Parce 1  Moves  in  Horizontal 
Dl rect ion  Only .  ) 
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clarity.  The  parcel  now  moves  to  position  4  and  expands 
adiabat ical  ly  to  a  new  temperature  of  T*~.  A  repeat  of  this 
process  puts  the  parcel  in  position  5  with  a  temperature  of 
T  .  Now  the  temperature  of  the  gas  is  lower  than  that  of 
the  plate  and  heat  (QM  flows  from  the  plate  to  the  gas 
parcel,  expanding  the  parcel  and  raising  its  temperature  to 
T  .  This  heat  flow  causes  the  plate’s  temperature  to 
decrease  at  position  X  .  The  parcel  moves  to  the  right, 
again  under  adiabatic  compression,  and  we  are  back  where  we 
started,  in  position  1  with  temperature  T  (ai^pted  from 
Wheatley,  et  al.,  1985). 

This  gas  parcel  cycle  i3  repeated  at  the  resonator’s 
operating  frequency.  For  the  STAR  this  frequency  is  between 
approximately  250  and  600  cycles  per  second,  depending  on 
the  gas  types  and  mixture  (specifically  the  speed  of  sound 
in  the  gas)  and  the  length  of  the  resonator. 

Notice  that  the  temperature  of  the  gas  parcel  at 
position  X=0  is  different  depending  on  which  direction  the 
gas  parcel  is  moving.  It  is  this  phase  shift  in  temperature 
relative  to  motion  that  produces  the  thermoacoustic  effect, 
as  ue  will  show  in  the  calculations  that  follow. 

It  is  also  important  to  note  that  in  order  for  heat 
to  flow  between  the  gas  parcel  and  the  plate,  the  parcel 
must  be  vertically  located  within  about  one  thermal 
penetration  depth  of  the  plate.  The  thermal  penetration 
depth  (6k)  j.s  the  distance  that  heat  can  diffuse  through 
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the  fluid  during  a  tine  1A>  ,  where  *  is  the  acoustic 
angular  frequency. 

Wheatley,  et  al.  (1986)  describe  thermoacoust ic 
engines  as  consisting  of  long  trains  of  these  gas  parcels, 
all  about  a  thermal  penetration  depth  from  the  plate.  The 
parcels  draw  heat  from  the  plate  at  one  extreme  of  their 
oscillatory  motion  and  deposit  heat  at  the  other  extreme. 
Adjacent  heat  flows  cancel,  except  at  the  ends  of  the 
plate.  The  net  result  is  that  an  amount  of  heat,  Q,  is 
passed  from  one  end  of  the  plate  to  the  other.  It  should  be 
noted  that  the  above  discussion  assumes  zero  viscosity. 

2.  Single  Plate 

Swift  (1989)  illustrated  the  basic  principles  of  the 
thermoacoustic  engine  by  examining  a  simple  example  of  a 
single  plate  in  a  gas-filled  tube  such  as  we  described 
above.  In  this  example  the  acoustic  and  thermodynamic 
effects  are  nearly  distinct.  In  the  absence  of  the  plastic 
plate,  the  plane  standing  wave  sustained  in  the  tube  has 
well  known  acoustic  properties.  It’s  pressure  and  x- 
direction  velocity  (along  the  length  of  the  tube)  are  given 
as  (using  peak  values): 

p ^=P^s in ( x/* )  ( 1 1  —  1 ) 

and  u,  =  i ( P  A/ f_a)cos( x/a)  ( 1 1—2 ) 

1  A  '  m 

where  P^  is  the  pressure  amplitude  at  the  pressure 
antinode,  *  =  */2f=aA>  is  the  reduced  wavelength,  a  is  the 
wavelength,  a  is  the  speed  of  sound  in  the  gas,  w  is  the 
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angular  frequency  and  is  the  mean  density  of  the  gas. 

The  subscript  1  indicates  the  first  order  tine  oscillation, 
and  the  subscript  a  indicates  the  aean  value.  The  total 
pressure  can  be  given  as: 

P-Pn+p!®1*1-  ( 1 1-3) 

A  similar  expression  may  be  used  to  represent  the  total 
velocity.  Following  Swift,  we  will  assume  that  the  first 
order  in  the  acoustic  amplitude  is  sufficient  for  all 
acoustic  and  thermodynamic  variables  (pressure,  velocity, 
temperature,  density  and  entropy). 

The  sound  wave  is  adiabatic  and  has  an  oscillatory 
temperature  given  by: 

V'V'VyPr  ,n‘4) 

where  p=- ( d f /dT )  / p  is  the  isobane  thereel  expansion 

r  * 

coefficient  and  cp  is  the  isobaric  (constant  pressure)  heat 
capacity  per  unit  mass.  Note  that  T^  and  p^  are  in  phase. 
For  ideal  gases,  thermodynamics  tells  us  that: 

T  (V?  c  =(y-l )T  /yP»,  ( 1 1-5) 

a  '  o  p  n  b 

where  *  =c  Vc„  is  the  ratio  of  isobaric  to  isochoric 

p*  v 

specific  heats.  *  equals  5/3  for  monatomic  gases  and  is 
smaller  but  still  greater  than  one  for  other  gases. 
Combining  Equations  1 1 —4  and  1 1 —5  gives  us: 

T. /T  =C (J-l)/y]pi /p  .  ( I 1-6) 

in  1  B 

The  introduction  of  the  plastic  plate  into  the 
standing  wave  modifies  the  original,  unperturbed  tempera¬ 
ture  oscillations.  This  modification  is  due  to  the  heat 
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flow  between  the  gas  and  the  plate,  as  described 
previously.  The  teaperature  is  modified  in  both  magnitude 
and  phase  for  gas  about  a  thermal  penetration  depth  away 

from  the  plate.  According  to  Swift,  this  results  in  two 
important  effects:  1)  a  time  averaged  heat  flow  near  the 
surface  of  the  plate,  along  the  direction  of  acoustic 

vibration,  and  2)  the  generation  or  absorption  of  real 

acoustic  power  near  the  surface  of  the  plate.  In 

refrigerators  the  effect  is  an  absorption  of  acoustic 
power.  In  prime  movers  the  acoustic  power  is  generated  by 
the  temperature  gradient  present  in  the  plate. 

There  are  several  assumptions  we  will  make  to 

simplify  calculations,  again  following  Swift.  We  assume 
that  the  plate  is  short  enough  compared  to  a  reduced 

wavelength  (Ax«*)  and  far  enough  from  both  velocity  and 
pressure  nodes  that  Pj  and  Uj  can  be  considered  uniform 
over  the  entire  plate.  We  assume  the  gas  has  zero 

viscosity,  so  that  Uj  does  not  depend  on  y  (ue  had  already 
made  this  assumption  by  looking  at  Uj  in  the  x-direction 
only).  We  assume  that  the  plate  has  a  large  enough  heat 

capacity  per  unit  area  that  its  teaperature  does  not  change 
appreciably  at  the  acoustic  frequency.  We  assume  that  the 
plate  has  a  mean  temperature  gradient  in  the  x-direction 
Finally  ue  neglect  the  plate’s  and  gas’  thermal 
conductivity  in  the  x-direction. 
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Applying  these  Assumptions,  we  see  that  the  mean 
fluid  temperature  (TB<x))  is  the  same  as  that  of  the  plate. 
Swift  calculated  the  oscillating  fluid  temperature  using 
the  general  equation  of  heat  transfer.  He  kept  only  first- 
order  terms,  neglected  thermal  conduction  along  x  and 
applied  the  boundary  condition,  ^(0)=®,  imposed  by  the 
plate.  The  resulting  equation  is: 

T .  (y  )  =  t  (T  P/ )p.-(7TM/w)u . 3 C 1-expC- 1  1+i )y/fi„13.  (II-?) 

1  I  9  m  p  i  n  i 

To  interpret  this  equation,  we  look  at  it  in  the  limit  that 
the  gas  is  far  enough  from  the  plate  <y»SK)  to  make 
negligible  thermal  contact  with  the  plate.  This  gives: 

Tx(y)  «C(TJ#P/f|icp)p1-(VTB/u.)u13.  ( 1 1-8) 

The  first  term  in  Equation  1 1-8  is  the  temperature 

oscillation  due  to  the  adiabatic  expansions  and 

compressions  of  the  gas  (see  Equation  II-4).  The  second 
term  comes  from  the  mean  temperature  gradient  in  the  gas. 
As  the  gas  oscillates  along  x  with  displacement  amplitude 
Uj/w,  the  temperature  at  a  given  point  in  space  oscillates 
by  an  amount  even  if  the  temperature  of  a  given 

piece  of  fluid  remains  constant.  The  actual  temperature 
oscillations  are  a  linear  superposition  of  these  two 
effects. 

The  y  dependent  part  of  Equation  1 1-7  is  complex.  It 
approaches  1  for  y»SK,  as  given  previously.  It  approaches 
zero  for  y«SK,  where  the  plate  imposes  the  condition  Tj=0. 
Most  importantly,  for  y *6*.  its  magnitude  is  still  on  the 
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order  of  1,  but  it  has  a  substantial  imaginary  part.  This 
phasa  shift  m  tha  oscillating  tamparatura  of  tha  standing 
wave  at  y * ,  dua  to  tha  tharmal  prasanca  of  tha  plats,  is 
an  important  rasult  bacausa  it  laads  directly  to  the  time 
averaged  heat  flow  in  tha  x-direction.  This  is  tha  same 
phasa  shift  that  was  emphasized  in  tha  description  of  the 
gas  parcel  motion. 

a.  Heat  Flow 

Swift  argues  that  since  we  are  neglecting  ordi¬ 
nary  thermal  conductivity  in  tha  x  direction,  the  only  way 
heat  can  be  transported  along  x  is  by  tha  hydrodynamic 
transport  of  entropy,  carried  by  tha  oscillatory  velocity 

V 

42*Tmf*5Tri=(  1/2)fmcpUCTi3ur  ( 1 1-9 ) 

The  line  above  the  quantity  SjUj  indicates  the  time-average 

of  the  first  order  entropy  and  velocity  product.  The  heat 
flow  is  a  second  order  quantity,  signified  by  the  subscript 
2. 

The  total  heat  flow  Q2  along  the  plate,  in  the  x 
direction,  is  found  by  integrating  q2  over  the  y-z  plane: 

WVy  di-la) 

where  TT  is  the  perimeter  of  the  plate  in  the  y-z  plane. 
Substituting  for  T^  and  performing  the  integration  gives: 

Q2=-<  1/4  )TTSx  (T^JpjUjt  r-1)  ( 11-11) 

where  7T&*  is  the  thermodynamically  active  area  in  a  plane 
perpendicular  to  the  longitudinal  acoustic  motion,  Tmp  is 
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the  heat  parameter  of  the  got.  And  r  is  defined  as  the 
rAtio  of  ActuAl  temperature  gradient  to  tho  critiCAl 
tesporAturo  gradient  (7T/vTcrit).  ^**e  critiCAl  temperature 
gradient  occurs  when  tho  temperature  change  Along  the  plate 
just  eatches  the  temperature  change  due  to  the  Adiabatic 
compression  of  the  gas,  and  no  heat  flows  between  the  gas 
and  the  plate.  This  is  the  boundary  between  the 
refrigerator  and  prime  mover  functions  of  thermoacoustic 
engines. 

Equation  1 1- 11  shows  that  when  r<l,  heat  flows 
up  the  temperature  gradient  from  cold  to  hot  and  work 
(acoustic  power)  is  absorbed,  as  for  a  refrigerator.  When 
r=  1  there  is  no  heat  flow.  When  r> 1 ,  heat  flows  down  the 
temperature  gradient  from  hot  to  cold  and  acoustic  power  is 
produced,  as  for  a  prime  mover. 

Note  that  the  total  heat  flow  is  proportional  to 

the  area  7T&K  ,  and  to  T^P  (which  equals  1  for  ideal  gases). 

It  is  also  proportional  to  the  product  PjUj,  and  so  equals 

zero  if  the  plate  is  at  either  a  pressure  or  velocity  node 

of  the  standing  wave.  The  maximum  value  of  the  product 

occurs  halfway  between  the  nodes.  Finally,  the  heat  flow  is 

proportional  to  the  temperature  gradient  factor  r-l.  For 

vT  >vT  .  r-l  is  greater  than  zero  and  the  heat  flow  is 
a  crit 

toward  the  pressure  node,  while  for  VT<VT _ . . ,  r-l  is  less 

m  cr  1 1 

than  zero  and  the  heat  flow  is  away  from  the  pressure  node. 
If  suitable  heat  exchangers,  at  and  Tg,  are  installed  at 


the  end a  of  tht  plot*  (with  fy*Tc«  Tft4x),  this  hoot  flow 
carries  hoot  froo  on*  exchanger  to  tho  other. 

Tho  hoot  flow  is  ssoll  under  ordinory 

circuoatoncos.  However,  in  closed  resonators  it  is  possiblo 
to  ochiovo  sound  ooplitudos  oony  or dors  of  magnitude  highor 
thon  thoso  of  ordinory  convorsot ion.  Sine*  Q2  is 
proportionol  to  P^*,  ond  smeo  in  procticol  acoustic 
engines  tho  ontiro  cross  section  of  tho  standing  wovo  is 
filled  with  plotos  (spaced  roughly  46*  opart),  very  high 
heat  flows  aud/or  a  largo  T^-Tg  *ay  bo  achieved, 
b.  Work  Flow 

The  work  flow  (i.e.  acoustic  power)  is  given  by 
the  work  per  cycle  tioes  the  rote  at  which  that  work  occurs 
(the  acoustic  frequency  f).  Froo  theroodynaoics,  the 
average  acoustic  power  produced  per  unit  voluoe  is: 

w2=-(w/?>)  ip~f1=-(  l/2)w*plIsCT13.  (11-12) 

The  gas  about  a  thermal  penetration  depth  away 
from  the  plate  "breathes",  because  of  thermal  expansion  and 
contraction,  with  the  right  time  phase  with  respect  to 
oscillating  pressure  to  do  (or  absorb)  net  work.  This  is 
exactly  the  same  gas  that  we  have  seen  is  responsible  for 
the  heat  flow.  Gas  elsewhere  is  ineffective  in  doing  (or 
absorbing)  work.  The  density  oscillations  for  y«S„  and  for 
y»&x  are  in  phase  with  the  pressure  oscillations,  and  hence 
do  (or  absorb)  no  net  work. 
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The  total  acoustic  power  produced,  W^,  is  found 

by  integrating  ovtr  *U  apace,  aa  with  the  boat  float 

w2«(  l/dm*  (f-lXpjV^aMr-lUx/x,  (11-13) 

uhoro  all  the  tins  have  been  defined  previously. 

The  acoustic  power  is  proportional  to  the  voluse 

tf&'&x  of  a  fluid  that  is  about  a  thermal  penetration  depth 

froa  the  plate.  It  is  proportional  to  p^*,  and  so  is 

quadratic  in  the  acoustic  aaplitude  (as  was  the  heat  flow) 

and  vanishes  at  pressure  nodes.  Finally  W2  is  proportional 

to  r-l,  the  sane  tenperature  gradient  factor  as  appeared  in 

Q2  (Equation  11-11).  When  VTm*vTcrlt,  i*-i*e,  and  there  are 

no  tenperature  oscillations  in  the  fluid,  other  than  those 

due  to  adiabatic  processes,  and  no  acoustic  power  is 

absorbed  or  generated.  For  7T*>vTcrlt»  •*-!>•  end  acoustic 

power  is  produced  near  the  plate.  Whether  this  power 

increases  the  amplitude  of  the  standing  wave,  is  radiated 

away  to  infinity,  is  simply  absorbed,  or  flows  through  an 

acoustic-to-electric  transducer  to  generate  electric  power 

depends  on  details  of  the  resonator,  not  on  the  plate 

itself  or  the  standing  wave  near  the  plate.  For  vT  <vT„..., 

■  cr  l  * 

r-l<0  and  acoustic  power  is  absorbed  near  the  plate.  For  a 
tube  without  plates,  which  has  a  diameter  less  than  the 
wavelength,  at  constant  temperature  (VT=8),  this  work  flow 
is  responsible  for  the  ordinary  thermal  attenuation  of  a 
sound  wave  (Kinsler,  et  al.,  1982). 
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c.  Efficiency 

He  coabine  Equations  I I- 11  and  11-13  to  tot  tha 
efficiency  of  the  plate  with  no  viscous  losses  or 
longitudinal  conduction  losses* 

j(*Axpl)/<faa*ul>. 

Since  Uj-u^sintx/x),  and  pt*  Mau#cos(x/*>,  we  get* 
*}*C(*-l)/(Tn*)3(  4x/(*tan(x/*) ) ). 

For  x« *,  this  reduces  to* 

7*CU-l)/(TM)3Ux/x), 

and  ,?*l7carnot^r* 

This  efficiency  approaches  the  Carnot  efficiency 

as  the  power  output  and  heat  transfer  rates  approach  aero. 

He  can  nake  a  stellar  calculation  for  the 
refrigerator  node  of  operation,  where  the  relevant 
efficiency  is  the  coefficient  of  perforaance  (COP).  Given 
that  COP-Qg/Hg,  ue  find  that* 


(11-14) 

( I I- IS) 

(11-16) 

(11-17) 


COP* r( COP 


carnot 


). 


(11-18) 


3.  Short  Engine 

Having  developed  the  slaplified  case  of  a  single 
plate  ue  now  turn  our  attention  to  a  more  realistic  model, 
with  an  entire  stack  of  plates.  Since  Suift  (1989)  has 
incorporated  the  thernoacoustic  theory  of  Rott  ( 1969,  1973, 
1974,  1975  and  1986),  we  will  once  again  follow  his 
developaent  closely.  The  assuaptions  that  the  length  of  the 
stack  is  such  shorter  than  a  wavelength  and  that  the 
temperature  spanned  is  auch  less  than  the  absolute 
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temperature  will  be  retained.  However,  this  more  realistic 
model  includes  viscosity,  longitudinal  thermal  conduction, 
finite  (instead  of  infinite)  plate  heat  capacity  and  a 
stack  of  plates  (instead  of  just  one).  The  goal  in  our 
development  is  to  include  these  effects,  which  produce  many 
complications,  in  such  a  manner  as  to  maintain  a  connection 
to  the  more  simplified,  but  more  easily  understood, 
development  of  the  previous  section. 

Swift  defines  a  co-ordinate  system  where  the  x-axis 
is  along  the  direction  of  acoustic  vibration  (longitudinal 
in  the  stack)  and  the  y-axis  is  perpendicular  to  the  planes 
of  the  plates  (using  a  parallel  plate  geometry)  with  the  y- 
origin  centered  between  two  plates.  This  gives  a  plate 
spacing  of  2y0  and  a  plate  thickness  of  21.  He  begins  with 
an  equation  of  motion  and  a  continuity  equation  for  the 
fluid: 

(<5  v/Jt )  +  (  v  v)v3=-vp+vv*v+(5  +  v/3)7(v  -  v) ,  (11-19) 

and  (<Jf/Jt)  +  V-(?v)=0,  (11-20) 

where  f  is  density,  v  is  velocity,  p  is  pressure,  v  is 
kinematic  viscosity  (Swift  uses  )*  for  this  parameter)  and  £ 
is  second  viscosity.  He  imposes  the  boundary  condition  that 
v=0  at  the  fluid-plate  interface.  He  then  uses  the  equation 
for  heat  flow  in  a  fluid  and  solid: 

fTC<Js/it+v-vs]=  v*()WT)  + (quadratic  velocity  terms),  (11-21) 

and  p  c  <<JT  /Jt)  =  A  v2T  .  (11-22) 

where  is  thermal  conductivity,  s  is  entropy  per  unit 
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mass,  T  is  absolute  temperature  and  the  subscript  s  refers 
to  the  solid.  The  applicable  boundary  conditions  are  that 


T=T  and  AiT/jy=  JT_/iy  at  the  fluid-plate  interface. 

......  S  .  S  S.  . ^  . 

Equations  for  each  of  the  needed  parameters  are  then 
linearized,  assuming  that  they  all  oscillate  at  the  angular 
frequency  u>: 


p=Vpi,x>® 

f=fB+  fl1*'?’® 


iu*t 


iu*t 


.  .  iu> t  ,  a  .  v  iu/t 

v=xUj(x,y)e  +yVj(x,y)e  , 

T=T  (x)+T. (x,y)elwt, 
m  l 


T  =T  (x)+T.  (x,y)e 
s  m  l ,  s 


iu*t 


iu;+ 

and  s=s  (x)+s, (x,y )e  , 

ill  1 


(11-23) 
( 11-24) 
(11-25) 
(II -26 ) 
(1 1-27) 
(11-28) 


where  the  subscript  m  denotes  a  mean  value  and  the 
subscript  1  refers  to  a  first-order  quantity.  Swift  then 
manipulates  these  equations,  by  integrating  with  respect  to 
y,  into  the  form  of  a  wave  equation  for  Pj(x)  and  an 
equation  for  energy  flux  H2  (which  includes  both  heat  and 
work  since  H=Q+H)  along  x.  The  wave  equation  is:  (11-29) 


( v  - 1 )  f  „ 


1  +  - 


l+< 


s 


V2  d 
V - 


1-fy  dPi 
Pm 


dx 


a 


dTm  d*l 


=0 


u> 2  dx 

and  the  equation  for  energy  flux  is: 
[* 

.  Try0  ,  dp 


u/2  (  1  —  <r )  (  1+  <  )  dx  dx 


H2= 


2«f 


Im 


m 


Tty^c 


i  ^  TP(fK  -fv  )  " 
1  .  g  m  * 

P 1  I  ^ 

dx  [  (  l+«_) ( 1+r ) 


0  p 


dTn  dPl  dp 


1 


1 


s 


Im 


2uA 


fm 


(  1—  <r )  dx  dx  dx 


yj 


(f„  -f.,  )(!+<)  (f,,/fJ 


fv+' 


(l+(  )(l+c) 
s 
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(11-30) 


dT 

-T(M+1V~r 

dx  . 

We  have  used  as  Prandtl  number, 

fM=tanhC(  1  +  i  )y0/Srt  ]/[  (  1  +  i  )y0/ 1 ,  (11-31) 

fv  =tanhC  (  1  +  i  )y0/Sv]/C  ( 1  +  i  )y0/Sv  ] ,  (11-32) 

and  i  =  c  SK  tanh  C  (  1  +  i )  y,./ SK  3/ f  c  4  tanh  C  (  1  +  i )  1  /  ] ,  ( 1 1-33 ) 

snip’'  ^  0  n  '  s  s  s  s 

where  S  is  a  specific  penetration  depth  denoted  by  the 
appropriate  subscript,  *  for  thermal  penetration  depth  in 
the  fluid,  s  for  thermal  penetration  depth  in  the  solid  and 
v  for  viscous  penetration  depth  in  the  fluid  (which  is 
roughly  the  distance  that  shear  momentum  diffuses  in  time 
\/*>) .  Swift  then  applies  some  simplifying  assumptions  to 
produce  equations  that  can  be  related  to  the  results  of 
section  2.  A  "boundary  layer"  approximation  (i.e.,  y0*>iK 
and  1»6  )  sets  the  hyperbolic  tangents  to  unity.  Since  in 
practical  engines  the  plate  spacing  is  2  to  4  thermal 
penetration  depths  (or  greater)  and  the  hyperbolic  tangent 
function  is  within  a  few  percent  of  unity  at  y0=24K,  this' 
is  a  reasonable  approximation.  A  "short  stack" 
approximation  (i.e.,  &x*A)  yields  relations  for  p1  and 
that  are  identical  to  those  of  section  2,  with  the 
exception  of  a  (l+l/y0)  term  in  u1  that  accounts  for  the 
velocity  being  higher  inside  the  stack.  This  approximation 
also  gives  us  the  fact  that  p^,  u^  and  all  thermophysical 
properties  are  independent  of  x  within  the  stack.  Finally, 
the  zero  viscosity  assumption  is  made  (although  viscosity 
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will  be  added  back  in  later)  thereby  setting  or  and  fv  equal 
to  zero.  Using  these  approximations  Equations  11-29,  11-30 
and  11-33  can  be  reduced  to: 

1  dp, "I  U-l) 

—  — i  - - — (r-l)p.  ,  (11-34) 

f,  dxj  (i+i)(i+«s)y0 

1  TmP  „  dT 

H  =--¥sK -JL.p  su.s(  r-l)-TT(y  A  +  :  X  )  — 2  ,  (11-35) 

*  4  l+<  1  1  w  s  dx 

s 

Q  C  S 

and  (=-”■-£—  ,  (11-36) 

s  f  c  5 
Is  s  s 

where  all  terms  have  previously  been  defined.  It  can  be 

seen  that  the  first  term  in  Equation  11-35  is  identical  to 

the  hydrodynamic  heat  flux  Qg  derived  for  the  single  plate, 

except  for  the  l/(!+(  )  term.  This  term  accounts  for  the 

finite  heat  capacity  of  the  plate  and  modifies  the  boundary 

condition  for  T1  at  the  fluid-plate  interface  to  a  non-zero 

value,  thereby  reducing  the  heat  flux.  For  the  short  stack 

U0  is  small  ( W/Q0=  r=&T/Tm r«l)  so  heat  flux  and 

2  22  carnot  m 

energy  flux  are  approximately  equal.  The  second  term  in 
Equation  11-35  accounts  for  ordinary  conduction  of  heat 
down  the  temperature  gradient  by  the  fluid  and  solid. 

The  next  step  is  to  derive  an  equation  for  the 
acoustic  power.  Since  net  acoustic  power  cannot  enter  or 
leave  the  engine  in  the  y-direction  the  acoustic  power 
generated  or  absorbed  must  involve  the  difference  in 
average  acoustic  intensity  p . u .  between  the  two  ends  of  the 
stack.  Using  this  and  the  fact  that  dp  ^/dx=-i  w^u  ^ ,  Swift 
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manipulates  Equation  11-34  to  arrive  at  a  relation  for 
acoustic  power: 

1  (*-l)u» 

W,-=-7rsHAx - —  |P il  2  <  r- 1 )  ,  (11-37) 

^4  f  aMl+i  )'  11 

;m  s 

which  is  the  same  as  Equation  11-13,  with  the  exception  of 

the  l/(l+(  )  term  which  was  discussed  above.  From  these 
s 

results  it  is  clear  that,  under  these  simplifying 
assumptions,  the  short  engine  performs  in  a  manner  similar 
to  that  predicted  by  the  single  plate  results. 

We  now  follow  Swift  in  repeating  this  derivation, 
but  include  non-zero  viscosity.  The  inclusion  of  viscosity 
adds  some  complexity  to  the  situation.  Under  these  circum¬ 
stances  Uj  becomes  a  function  of  y  and  the  concept  of 
critical  temperature  gradient  becomes  less  well  defined. 
Since  Tj  now  depends  on  y,  in  an  irregular  fashion,  there 

is  no  natural  definition  of  *T  ...  However,  we  will  use 

cr  it 

(as  did  Swift)  r=7T  /7T  and  ,  =T  Pp.S/P  c  <u.S> , 

m  crit  crit  m  M  >m  p  1 

where  the  superscript  s  defines  the  value  for  the  standing 
wave.  The  addition  of  viscosity  modifies  the  relation  for 
dp^/dx  to: 

s 

dp  ,  iu/o  <  u , '  u>q  <  u ,  > 

— i - i»— 1—  ,  (11-38) 

dx  l-f„  l-fv 

5 

which  is  complex  even  though  <Uj  >  is  real. 

Substituting  into  Equations  11-29  and  11-30  and 
applying  the  boundary  layer  approximation  yields  the 
following  relations: 
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r  & 

m 

l+<r  +  <r+  ff  i 

P  S 
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l+ff5 — 

L  y<2L 
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l+ff5 
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-  <y0A+lA 
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( 1 1 -39  ) 


and 


p  +!sll  _ 


d 

1-f, 

*2i 

Of-l)  shPi 

r 

dx 

-  ^m 

d  x 

( i+i)  (  i+«s  )y0 

(  l+ff'4)  (  l-fv  ) 

—  mt 

. ( 11-40) 


For  the  zero  viscosity  case,  ff=0,  these  relations  reduce  to 
Equations  11-34  and  11-35. 

Work  flow  is  calculated,  using  the  short  stack 
approximation,  from  the  wave  equation.  This  gives: 


s 


w_= -TrsK  ax 
z  4 


(*-l)u»(pi  ) 

fa* ( l+<  ) 
/  m  s 


—  1 


<  l‘+»  )  <  l-sv/y0+6v  V2y0* ) 


1  <u,  >* 

-  -iia,  ax——  - - 

4 


(11-41) 


which  reduces  to  Equation  11-37  when  <r=0. 

From  these  relations  it  is  clear  that  viscosity 
makes  the  situation  much  more  complex.  Since  the  goal  is  an 
expression  for  the  coefficient  of  performance,  it  is 
necessary  to  determine  Qg/Wg.  We  know  that  H2=Q2+W2>  suc^ 
that  C0P=Q2/W2=(H2-W2)/W2=(H2/W2)-1.  In  order  to  maintain 
the  essential  effects  of  viscosity,  while  reducing  the 
complexity  of  the  resulting  equations,  we  will  include 
viscosity  only  to  the  lowest  order  present,  that  is  ff  . 
Additionally  we  will  assume  that  y0ssK,  which  is  the  case 
for  most  practical  engines  as  previously  discussed.  This 
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gives  us  that  $W/SK  **» .  Applying  these  approximations 


to  Equations  I I —40  and  11-41  yields: 

1  T  Ap .s<u .S>  dT 

H  =~ ffi,  -i-4 - L_-<r-l)-ir(y  *  +  lx  )— B  , 

*  4  <l+<  )(l-ff*)  w  s  dx 

s 


( 11-42) 


.  l  (y-i)w(p  s)i  l 

W0=-TTSHiix - i - (  r-l)--TTSvixu;P  <tt,S>i  ,  (11-43) 

2  4  eaMl+«  )  4  ®  1 

>m  s 

Now  we  can  take  the  ratio  of  these  two  relations  and 
get  a  useable,  understandable  relation  for  the  coefficient 
of  performance,  which  is  what  we  are  trying  to  improve 
experimentally.  The  result  is: 


- ( 1/4 ) TTSkT  pp1s<u1s>(  r-l) 
m  i  l 

(  1+7”  )  (  l-<r*> 

S 


-TT(y  a+1A  )— a 
v  s  dx 


-1  . (11-44) 


( l/4)7r$KAx(*-l)u><Pj  )  *  (  r-1)  1 

-■ - - — - — 7T&V  Ax<4>P  <u,S>  2 

f  aMl+<  )  4  tt  1 

1  m  s 

The  above  relations  are  approximate  expressions,  which 
should  apply  to  realistic  situations,  with  the  only 
remaining  .  assumptions  being  1)  the  short  stack 
approximation,  2)  the  boundary  layer  approximation,  and  3) 
a  small  Prandtl  number  approximation.  We  have  previously 
explained  that  the  boundary  layer  approximation  is 
reasonable  for  the  plate  separations  of  our  stacks.  The 
short  stack  approximation  is  also  not  particularly 
restrictive  for  realistic  engines.  With  a  Prandtl  number  of 
approximately  0.027  for  helium-xenon  we  have  a  relatively 


small  Prandtl  number,  so  this  approximation  is  also  very 
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reasonable.  In  his  comparison  to  theory,  Hofler  (1986) 
shows  that  in  relatively  "long"  acoustic  engines  Equations 
11-29  and  11-39  predict  the  performance  within  29  to  39 
percent . 

4.  Design  Considerations 

The  purpose  of  this  section  is  to  discuss  the 
primary  modification  that  we  have  made  based  on  theory, 
that  is,  the  reduction  of  the  Prandtl  number  of  the  system. 
To  evaluate  how  this  affects  performance  we  will  first  look 
at  Equation  11-44.  We  know,  from  experience,  that  the 
second  term  in  the  numerator  is  small  compared  to  the  first 
term.  We  also  know  that  *  is,  in  general,  small  compared 
to  unity,  so  we  will  set  *  =0  for  evaluation  purposes.  The 
product  Tn?  is  known  to  equal  unity  for  ideal  gases. 
Applying  these  conditions,  with  Swift’s  Sv  =  SH(r'  and  y0=s6# 
we  can  reduce  Equation  11-44  to?  (11-45) 

"{C-H(Ty0.p1S<u1S>(  r-1)  ]/(  l-ir’Sj-CTUy^+lA  )dT  /dx3  1 
COP=  - 2_J - i_ - - m  ■  -1  . 

UH^axPf-DpjM  r-l)3/(?  a*  >-CHiry0* 

£ 

T‘  is  clear  from  this  that  as  <r  is  reduced,  1  —  a-  is 
increased  resulting  in  an  overall  decrease  in  the  numerator 
of  our  relation,  and  »  is  reduced  resulting  in  a  decrease 
in  the  denominator  of  our  relation.  These  two  effects,  when 
evaluated  fully,  show  that  COP  will  decrease  with  reduced 
Prandt’  number.  However,  we  know  based  on  our  experiments 
that  this  cannot  be  true  for  our  system.  This  effect  that 
is  apparently  the  opposite  of  what  we  expect,  and  know  to 
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be  true,  is  therefore  dua  to  some  (or  one)  of  the 
approximations  that  we  have  made  being  invalid  for  our 
system.  As  ue  are  about  to  show,  the  relations  given  in 
Equations  H-39  and  11-41  yield  the  proper  result, 
indicating  that  the  inclusion  of  viscosity  only  to  the 
lowest  order  does  not  provide  adequately  for  the  effects  of 
viscosity  in  our  system. 

Looking  at  Equations  11-39  and  11-41  we  see  that 

they  are  fairly  complex.  We  will  use  the  simplifications 

that  «  =0,  T  Pal  and  f  a*=  3fp_ ;  and  ue  will  consider  the 
s  m  'mm 

thermal  conduction  term  in  the  heat  flow  to  be  negligible. 

u  o 

Ue  also  know  that  &„*&*/* ' ,  and  we  will  use  Pj  ap^coskx  and 

<u1s>=a(  1/D (p^/p^Jsinkx  to  further  simplify  our  relations. 

This  allows  us  to  re-write  Equations  11-39  and  11-41  ass 

4H0Yp  acoskxsinkx 

^  m _ 


1TSKpfti  <l+*)Cl-Sv/yft+J4($v/y0)*] 


0 

(11-46)  and 
4U  *p 

— — -=  ( V  - 1 ) 
P0jTr6« 


l+ir^+ir 


l+<r 


-<l+^-$v/y0) 


uixcos 2  kx 


(l+<r*)(l-Sv/y0+b($v/yft) 


-1 


0‘ 


K 

^Axff  sinJkx 


c  1_Su/y0+*(  sv/y0)  n 


(11-47) 


Now  we  can  combine  Equations  11-46  and  11-47  to  get  a 
relation  for  COP  as  follows? 
coskxs inkx 


COP=- 


kx (  l+<r)  (  l-5v/y0+b<  iy/yft)  2  ) 


- -1  . 


(11-48) 


r  0 
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where  the  brackets  < - >  represent  the  following  relations 


and  we  have  used  a/w Ax«a/wxs  1/kx.  It  is  not  clear  from  this 
relation  how  COP  will  vary  with  Prandtl  number.  However,  if 
we  define  a  set  of  variables  and  perform  a  numerical 
analysis  we  can  see  the  desired  effect.  He  will  define  a 

U 

quantity,  6^= ( SK )  =constant ,  as  the  geometric  mean  of  the 

u  u 

penetration  depths,  such  that  =  and  S^/y^) *  . 

This  roughly  corresponds  to  our  experiment.  Setting 

reasonable  values  for  the  four  parameters  in  our  relation: 

kx=0.2,  r=0.7,  S^/y^aC.5  and  Y=1.667,  we  can  see  how  COP 

varies  with  Prandtl  number.  Figure  II-2  is  a  plot  showing 

hou  heat  (Hg),  work  (Hg)  and  COP  vary  wi  tx  Prandtl  number, 

using  Equations  11-46,  11-47  and  11-48.  It  is  clear  that  as 

Prandtl  number  is  reduced  COP  increases.  It  should  be  noted 

that  as  Prandtl  number  becomes  very  small  we  actually  see 

COP  decrease  with  reduced  Prandtl  number,  as  we  estimated 

using  the  approximation  with  only  the  lowest  order 

viscosity  terms.  Since  Swift’s  engine  (1989)  with  liquid 

sodium  has  a  Prandtl  number  of  approximately  0.004,  his 

approximation  is  in  fact  valid  for  his  system,  although  it 

is  not  for  ours,  as  we  discussed  above. 
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b.  BINARY  GAS  HIXTUURE  THEORY 
1.  Sound  Speed  Theory 

One  of  the  parameters  of  importance,  when  consi¬ 
dering  refrigerator  performance,  is  the  velocity  of  sound 
in  the  fluid  being  used.  The  purpose  of  this  section  is  to 
show  how  this  value  is  calculated  for  binary  gas  aixtures. 

In  general,  the  sound  speed  *a”  in  a  gas  is  given  by 
the  relation?  a*(YR^T/a)  ,  where  t*cp/cy  is  the  ratio  of 
isobanc  to  isochoric  specific  heats,  Rq  is  the  universal 
gas  constant,  T  is  the  absolute  temperature  and  a  is  the 
nolecular  weight  of  the  gas.  This  relation  holds  for  a 
binary  Mixture  of  gases  as  well.  However,  Y and  m  may  be 
different  for  a  mixture  and  must  be  calculated.  From  Huang 
(1976)  we  get  a  relation  for  the  specific  heat  capacity  of 
a  binary  Mixture;  (cp)  t  j*CxiMicpi+XjM jcpj  l/Cx^+x^Mj  1 , 
where  xt  is  the  mole  fraction  of  the  i  species,  Hi  is  the 
molecular  weight  of  the  ith  species  and  cpi  is  the  isobanc 


specific  heat  capacity  of  the  i 


th 


species.  A  similar 


relation  holds  for  cy  as  well.  The  equivalent  molecular 

weight  m,  of  the  mixture  is  given  by: 

Combining  the  above  relations  we  get  a  final  equa- 

ion  for  the  speed  of  sound  of  a  binary  mixture,  as  follows: 

a={CRftT(x,  H.  c_.  ,+x  H.c  )  1/C  (x.H.+x.H. )  (x.H.c  +x  ,H.c„  .  )  ]>^. 

®  i  i  pi  j  j  pj  i  i  j  j  i  i  vi  j  j  vj 

Since  Y  is  approximately  5/3  for  any  ideal  monatomic 
gas  it  is  clear  that  any  variation  due  to  c^  and  cv  will  be 
small.  Therefore,  the  primary  effect  on  sound  speed,  of 
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mixing  two  gases  is  that  it  is  inversely  proportional  to 
the  square  root  of  the  total  mass.  At  3*6  degrees  Kelvin 
the  sound  speeds  of  the  mixtures  used  in  our  experiment  are 
as  follows: 

pure  helium  asie®8  m/s 

helium-argon  ( 19X  Ar)  a  -614  m/s 

helium-xenon  ( 12X  Xe)  a»46Q  m/s 

2.  Prandtl  Humber  Theory 
a.  General  Discussion 

Before  going  into  the  details  of  how  Prandtl 
number  varies  for  binary  gas  mixtures,  it  is  important  to 
understand  the  general  concept  of  Prandtl  number,  and  to 
discuss  the  expected  variation  based  on  this  concept. 
Prandtl  number  is  a  dimensionless  ratio  that  relates 
transport  properties  of  fluids.  In  Fluid  Mechanics  it  is 
given  the  symbol  "Pr*  and  in  Thermodynamics  we  use  *»". 
Prandtl  number  can  be  represented  in  many  different  ways. 
For  calculat ional  purposes  we  define  Prandtl  number  as 
f=l?Cp/A,  where  y  is  dynamic  viscosity,  cp  is  the  isdbaric 
specific  heat  capacity  and  A  is  the  thermal  conductivity  of 
the  fluid.  This  relation  is  used  because  these  parameters 
are  well  tabulated  and  more  readily  calculated  than  others 
that  are  available.  For  an  understanding  of  the  physics  of 
Prandtl  number  and  how  it  relates  to  the  refrigeration 
process  we  define  Prandtl  number  as  <t=v/k,  where  v  =  ^/f  is 
the  kinematic  viscosity  and  K  =  Vfcp  is  the  thermal 
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diffuaivity.  Thia  definition  ia  uaed  because  kineaatic 
viacoaity  represents  aoaentua  diffusion  through  the  fluid 
due  to  a  gradient  in  the  velocity,  and  theraal  diffuaivity 
represents  the  diffusion  of  theraal  energy  due  to 

teaperature  gradients  in  the  fluid.  Thia  aaya  that  Frandtl 
nuaber  is  a  ratio  of  aoaentua  diffusion  to  theraal 
diffusion.  A  third  nay  of  describing  Frandtl  nuaber  is 
<r»(  sv/6„  )  * ,  where  is  the  viscous  penetration  depth  and 

is  the  theraal  penetration  depth.  The  capacity  for 
diffusion  of  a  given  quantity  is  directly  related  to  the 
aagnitude  of  the  associated  penetration  depth.  Since 
&v  ■  t 2v/u» )  is  a  function  of  kineaatic  viscosity,  and 

5h  =  (2K/u>)*  is  a  function  of  theraal  diffusivity,  we  will 
restrict  the  discussion  to  v  and  k  effects. 

Hirschfelder  et  al.,  (1954,  p. 14)  derive  a 

relation,  using  siaple  kinetic  theory,  for  the  theraal 
conductivity  of  a  lonatoaic  gas,  **158^7/4H,  where  is 
the  universal  gas  constant  and  N  is  the  aolecular  weight 
ofthe  gas.  Froa  this  relation  and  the  knowledge  that  c^-> 
cy=F  and  Y=Cp/cv,  where  R=F^/H,  it  is  easily  shown  that 
,=l?cp/A=4),/ ( 15(V-1) )  for  a  aonatoaic  gas.  It  is  also  known 
that  r  =5/3  for  any  aonatoaic  gas,  which  results  in  <r=2/3 
for  all  aonatoaic  gases.  Calculations  using  rigorous 
kinetic  theory  also  show  that  Prandtl  nuaber  is  approxi- 
aately  2/3  for  aonatoaic  gases.  This  theory  will  be 

discussed  in  the  next  section. 
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When  two  monatomic  gases  are  combined  the 


molecular  interactions  produce  changes  in  the  physical 
properties,  resulting  in  Frandtl  number  deviating  from  the 
2/3  value.  For  a  mix  of  a  heavy  gas  and  a  light  gas,  as  we 
have  used,  the  Prandtl  number  drops  rapidly  as  the 
percentage  of  heavy  gas  increases  from  2ero.  This  is  due 
primarily  to  a  very  rapid  rise  in  the  molecular  density 
(  =grams/mole)  of  the  gas  as  the  concentration  of  the 
heavier  gas  increases  from  zero  to  some  small  percentage. 
As  an  example,  a  mix  of  two  percent  xenon  with  98%  helium 
produces  a  64%  increase  in  this  molecular  density.  Since 
the  kinematic  viscosity,  V  =  ,  is  inversely  proportional 
to  this  molecular  density,  while  dynamic  viscosity  is 
relatively  invariant,  there  is  a  rapid  drop  in  kinematic 
viscosity  and  therefore  momentum  diffusion.  Thus,  we  see 
that  the  mi.  ^entum  diffusion  process  is  dominated  by  the 
presence  of  the  heavier  gas  molecules.  The  thermal 
diffusivity,  K=A/fCp,  on  the  other  hand  is  inversely 
proportional  to  the  product  fcp,  which  is  essentially 
constant  with  changes  in  concentration,  and  directly 
proportional  to  the  thermal  conductivity,  thereby  changing 
strictly  due  to  variations  in  thermal  conductivity.  Reif 

K 

(1965)  represents  thermal  conductivity  byA=c(kT/m)  ,  where 
c  is  a  heat  capacity  per  molecule  (which  is  essentially  a 
constant),  k  is  the  Boltzmann  constant,  T  is  the  absolute 
temperature  and  m  is  the  molecular  weight.  It  is  clear  from 
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this  relation  that  the  thermal  conductivity,  and  therefore 
thermal  diffusivity  and  thermal  energy  diffusion,  will  be 
much  greater  for  a  lighter  gas,  and  their  values  will  be 
dominated  by  the  concentration  of  the  lighter  gas  in  the 
mixture.  Combining  these  two  effects  it  is  clear  that  both 
momentum  and  thermal  energy  diffusion  rates  decrease  with 
increasing  concentration  of  the  heavier  gas.  However,  the 
rate  of  change  in  momentum  diffusion  is  much  more  dramatic 
at  the  lower  concentrations  and.  drops  off  rapidly  as  the 
concentration  of  the  heavier  gas  increases,  while  the  rate 
of  change  of  thermal  energy  diffusion  is  much  more  gradual 
and  continuous  over  the  concentration  range.  This  results 
in  the  ratio  of  the  two,  which  is  effectively  the  Prandtl 
number,  dropping  rapidly  from  the  2/3  value  for  the  pure 
light  gas  as  the  heavier  gas  is  initially  introduced, 
reaching  a  minimum  as  the  momentum  diffusion  steadies  out 
and  rising  back  to  the  2/3  value  as  energy  diffusion 
continues  to  decrease  while  the  heavier  gas  concentration 
approaches  10054. 

b.  Kinetic  Theory  and  Quantitative  Analysis 

The  purpose  of  this  section  is  to  discuss  the 
theory,  with  its  assumptions  and  limitations,  and  the 
equations  used  to  calculate  the  Prandtl  number  of  a  binary 
gas  mixture.  Much  of  the  information  presented  here  is 
reproduced  in  Appendix  A  in  order  to  allow  the  removal  and 
use  of  Appendix  A  without  reference  back  to  the  main  body 
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of  the  text.  Hirschfelder ,  et  el.  (1954)  introduce  a  basic 
kinetic  theory  in  Chapter  one,  develop  a  rigorous  kinetic 
theory  in  Chapter  seven  and  present  equations  for  the 
calculation  of  thermal  conductivity  and  dynamic  viscosity 
for  binary  gas  mixtures  in  Chapter  eight  of  their  text. 
Using  these  relations  with  another  for  the  isobaric 
specific  heat  capacity  of  the  mixture  we  can  determine 
Prandtl  number,  <r=rjc^/A..  Rather  than  reproduce  the  work  of 
Hirschfelder,  et  al.  (1954)  we  will,  present  a  discussion  of 
the  methods  they  used  in  deriving  their  equations, 
including  the  assumptions  and  limitations  that  apply.  Then 
we  will  present  all  of  the  equations  necessary  for  the 
calculation  of  Prandtl  number  for  a  binary  gas  mixture,  and 
a  sample  curve  of  Prandtl  number  as  a  function  of 
concentration  that  was  created  using  the  output  data  of  the 
program  presented  in  Appendix  A,  which  uses  these 
equations. 

The  development  of  the  kinetic  thoery  of  gases 
is  based  on  knowledge  of  the  distribution  function  that 
represents  the  number  of  molecules  of  a  specific  species, 
in  a  unit  volume  element,  about  a  point  in  space,  with 
velocities  in  a  unit  range  about  a  specified  velocity,  at  a 
given  instant  in  time.  This  distribution  function  is 
represented  by  f^r.v^t).  Hirschfelder,  et  al.  derive  the 
Boltzmann  equation,  in  terms  of  this  distribution  function, 
which  is  the  basis  for  discussion  of  transport  properties 
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of  gases.  The  Boltzmann  equation  is  given  as: 

’(fi,fj'-fifj)gij1>dbd(dVj(n-49i 


11 u 

'  If  \ 
u  1 

i  V  i{i 

Jt  \ 

, 1  iry 

*  *A  ijvi 

where 
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superscr 

8iJ 

is  the  pre-co 
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the  impact  parameter,  mA  is  the  mass  of  the  it^1  species,  X 
is  an  external  force  on  the  molecule,  V  is  the  velocity  of 
the  molecule  and  <  is  a  rotational  impact  parameter.  The 
hydrodynamic  equations  of  change  (conservation  of  energy, 
conservation  of  momentum,  etc.  )  are  derived  from  the 
Boltzmann  equation. 

Hirschfelder ,  et  al.,  use  a  perturbation  method 
of  successive  approximations,  which  they  attribute  to 
Enskog,  to  obtain  solutions  to  the  Boltzmann  equation. 
Using  this  method  they  derive  a  form  for  the  perturbation 
function  *,  that  depends  on  space  and  time  through  the 
quantities,  number  density,  mass  average  velocity  and 
temperature,  and  their  spatial  derivatives: 


i  s-U^lnT/arJ-tB^V^D+nE,  (C  'J'd,).  (11-50) 


i  i  '  ' "i  ’0' 

The  coefficients  Bi 


.  ...  (J>u 
J  i  J 


and  Ci<^)  each  represent  rather 


involved  integral  equations  which  must  be  solved.  They  are 
solved  by  a  series  approximation  method,  using  the  infinite 
Sonine  polynomial  series.  The  Sonine  polynomials  are 


defined  as 
„  (m) 


n 


(x)*E  C (-1 ) J (m+n) !xJ 1/C (n+j ) ! (m-j ) ! j ! . 


(11-51) 


Linear  combinations  of  a  finite  number  of  these  polynomials 
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are  shown  to  be  the  solutions  to  the  integral  equations  for 

At,  Bi  and  C^*.  However,  the  transport  coefficients 

(which  is  what  we  are  after)  are  expressed  in  terms  of  the 

Sonine  polynomial  expansion  coefficients.  These  Sonine 

polynomial  expansion  coefficients  involve  complicated 

combinations  of  integrals,  called  "bracket"  integrals. 

Chapman  and  Cowling  (1939)  express  these  bracket  integrals 

as  linear  combinations  of  a  set  of  integrals: 

ftiJ  ( 1,s)  =  (2fkT/HiJ  )  *  (e"yi  j  * )  (2s+3)  ( l— cos 1  (X)  )bdbdKiJ  , 

Jo  Jo 

(11-52),  where  )»  is  the  reduced  mass  of  the  colliding 

molecules  i  and  j,  is  the  reduced  initial  relative 

speed  of  the  colliding  molecules,  b  is  the  impact  parameter 

and  X  is  the  angle  of  deflection  of  the  molecules. 

(1  s ) 

Hirschfelder,  et  al.,  have  tabulated  values  of  these  ft  ’ 

integrals  for  frequently  encountered  bracket  integrals  in 

Table  I-M  of  their  text.  Hirschfelder,  et  al.,  go  on  to 

(1  s ) 

define  the  transport  coefficients  in  terms  of  these  ft 
integrals.  It  is  these  relations  that  we  have  used  in  our 
determination  of  Prandtl  numbers  for  binary  gas  mixtures. 

Before  presenting  these  relations  it  is  impor¬ 
tant  to  understand  the  assumptions  used  in  this  derivation 
and  the  limitations  they  impose,  if  any.  Since  the  usual 
definitions  of  transport  coefficients  only  apply  under 
conditions  only  slightly  displaced  from  equilibrium,  our 
results  are  limited  by  this  restriction. 
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Only  binary  collisions  are  considered,  which 
Halts  us  to  densities  low  enough  that  three-body 
collisions  are  unimportant.  The  theory  predicts  that 
viscosity  is  independent  of  pressure  at  constant 
temperature.  Hirschfelder ,  et  al. ,  state  that  experiments 
on  nitrogen  have  shown  that  only  a  4X  deviation  in 
viscosity  occurs  as  pressure  is  varied  from  1  to  60 
atmospheres.  Our  operating  condition  of  10  atmospheres  is 
well  within  this  range  and  should  therefore  closely  conform 
to  the  theory  in  this  respect. 

The  theory  limits  itself  to  Classical  Mechanics 
and  thereby  precludes  low  temperature  situations  where 
quantum  effects  become  significant.  These  quantum  effects 
are  less  than  IX  for  helium  above  200  Kelvin  (and  even 
smaller  for  heavier  isotopes).  The  coldest  temperature  we 
have  achieved  to  date  has  been  193  Kelvin,  under  a  no  load 
condition,  so  our  normal  operations  are  all  above  200 
Kelvin  and  we  are  not  concerned  with  quantum  effects. 

The  Chapman-Enskog  (1939)  theory  uses  a  first- 
order  approximation  in  their  perturbation  method  solution 
to  the  Boltzmann  equation.  This  is  only  valid  when  the 
gradients  of  the  physical  quantities  (density,  velocity, 
temperature,  etc.)  are  small.  This  means  the  changes  in 
these  properties  over  one  mean  free  path  must  be  small  with 
respect  to  unity.  Based  on  the  mean  free  paths  for  gases  at 
pressures  above  one  atmosphere,  this  condition  is  satisfied 
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except  under  conditions  of  extreme  gradients,  such  as  those 
present  in  shock  waves.  Since  the  observed  pressure 
waveforms  in  our  resonators  are  sinusoidal,  we  do  not  have 
any  shock  waves  and  we  are  not  restricted  by  this 
assumption. 

The  theory  further  assumes  that  the  dimensions 
of  the  containing  vessel  and  any  obstacles  within  it  are 
large  compared  to  the  mean  free  path,  which  provides  no 
limitations  for  our  gases  and  pressures.  Although  the 
'Jhapman-Enskog  theory  applies  strictly  to  monatomic  gases 
the  results  have  been  shown  to  be  good  even  for  polyatomic 
gases,  provided  that  the  molecules  are  not  too  non- 
spherical.  Finally,  the  development  of  the  theory  results 
in  a  set  of  equations  for  transport  coefficients  which  rely 
upon  the  previously  discussed  integrals.  These 
integrals  are  dependent  upon  the  potential  function  that 
represents  the  molecular  interactions.  Throughout  the 
derivation  by  Hirschf elder ,  et  al.,  and  in  all  of  the 
equations  used  for  our  calculations,  the  Lennard-Jones  (6- 
12)  potential  is  assumed. 

Now  that  we  have  presented  the  theory  and  its 
limitations  we  will  present  the  equations  necessary  for  the 
calculation  of  dynamic  viscosity,  thermal  conductivity  and 
isobaric  specific  heat  capacity  of  a  binary  mixture  of 
monatomic  gases.  From  these,  as  we  have  already  discussed, 
we  get  Prandtl  number,  where  <r=^c  /*. 
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The  dynamic  viscosity  for  a  binary  gas  mixture 


is  given  by: 


where 
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In  these  equations  x^  and  x2  represent  the  mole  fractions 
of  species  1  and  2,  Mj  and  »2  represent  the  molecular 
weights  of  species  l  and  2  and  A124=ft(2*2) 4/ft( *’ 4,  where 
the  superscript  ♦  refers  to  a  value  that  takes  into  account 
the  deviation  of  a  particular  molecular  model  (in  our  case 
the  Lennard-Jones  (6-12)  potential)  from  the  idealized 
rigid  sphere  model  of  molecular  interactions.  The  values 


for 

substance : 


and  C?23l  are  8iven  ky  'the  equation  for  a  pure 


CT  3  j  x  107=266 ,  93  ( HT )  */ C  <r  *  ft (  2 1  2  *  4  ( T* ) , 
where  ^  is  dynamic  viscosity  in  gm/cm-sec,  T  is  temperature 
in  ®K,  T*=kT/(  is  reduced  temperature,  M  is  molecular 

weight,  9  is  collision  diameter  in  angstroms  and  </k  is 
the  potential  parameter  in  ®K,  The  remaining  term  C*f123i’ 
represents  an  artificial  quantity  for  a  hypothetical  pure 
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substance  with  Molecules  of  Molecular  weight  2M1M2/(M1+M2) , 
that  interact  according  to  a  potential  curve  specified  by 
interaction  parameters  *12  and  <12,  These  teres  are 
described  in  teras  of  the  Lennard-Jones  (6-12)  potential, 
where  <t  is  the  distance  of  closest  approach  (in  angstroas) 
of  two  Molecules  which  collide  with  zero  initial  relative 
kinetic  energy,  and  i  is  the  maximum  energy  of  attraction 
of  the  two  Molecules.  C^123j  is  given  by: 
C712]1xl07*266.93C2M1M2T/(M14-H2)]VC(r12tft<2’2)4(T12*), 
where  the  subscript  12  on  individual  teres  refers  to  a 


combined  value  for  species  1  and  2.  The  values  for  <r12  and 
are  given  by  the  eepirical  combining  laws  t rt  +  <r0) 


b  * 

and  €12=(<1^2)  *  respectively,  and  TJ2 


12  12 
is  given  by 


T,2  =kT/(12. 


The  thermal  conductivity  for  a  binary  mixture 


is  given  by: 
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The  definitions  fop  all  terns  are  similar  to  those  given 
for  the  dynamic  viscosity  relations. 

The  remaining  quantity  necessary  for  the  calcu¬ 
lation  of  Prandtl  number  is  the  isobaric  specific  heat 
capacity  of  the  binary  mixture.  Huang  (1976)  gives  a 
relation  for  this  quantity*  cp*ixsEilicpi/,,S!SxiMicpi/SxiMi ' 


where  x. 


is  the  mole  fraction  of  species  i,  H  is  the 


and  c  .  is  the  isobaric 
Pi 


molecular  weight  of  species  i 
specific  heat  capacity  of  species  i.  The  required  units  of 
this  quantity,  for  compatibility  with  the  previous 
equations  are  cal/gm-*K. 

All  of  these  equations  were  formed  into  the 
Fortran  computer  program  presented  in  Appendix  A.  Outputs 
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have  been  created  (from  this  program)  for  temperatures  of 
200'K  and  3$0*K,  and  for  Mixtures  of  helium-xenon  and 
hel iua-argon,  for  all  concentrations  froa  1GGX  heliua  to  Ox 
helium.  A  composite  plot  showing  all  of  these  cases  is 
provided  as  Figure  II-3.  It  is  clear  from  these  curves  that 
the  Prandtl  number  variation  with  concentration  is  just  as 
predicted  by  the  qualitative  analysis,  given  in  the  first 
part  of  this  section. 

3.  Selection  of  a  Gas  Mixture 

Now  that  we  have  shown  how  the  sound  speeds  and 
Prandtl  numbers  of  binary  gas  mixtures  vary  with  concentra¬ 
tion,  it  is  necessary  to  consider  how  each  of  these 
parameters  affect  performance.  We  know  from  the 
thermoacoustic  theory,  presented  in  section  A  of  this 
Chapter,  that  power  density  will  decrease  with  reduced 
sound  speed  (due  to  the  increased  mass),  and  that  COP  will 
increase  with  reduced  Prandtl  number.  However,  as  Prandtl 
number  decreases,  due  to  the  addition  of  a  heavy  gas,  the 
overall  mass  of  the  gas  increases,  and  sound  speed  is 
reduced.  It  is  therefore  necessary  that  the  mixture  have  as 
low  a  Prandtl  number  as  possible  while  the  mass  increase  is 
kept  to  a  minimum  (and  therefore  the  reduction  in  sound 
speed  is  minimized).  Looking  at  Figure  1 1 —3  it  is  clear 
that  we  need  a  mixture  that  is  at  a  point  on  the  very 
steeply  sloped  section  of  the  curve,  as  the  concentration 
begins  to  deviate  from  pure  heliua.  Our  selection  of  12. 5X 
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xenon  and  18. Sx  argon  in  our  helium-xenon  and  helium-argon 
Mixes,  respectively,  was  Made  baaed  on  this  requirement. 
With  the  xenon  we  achieved  approximately  9®x  of  the  maximum 
possible  Prandtl  nuaber  reduction,  while  getting  only  40X 
of  the  possible  Hass  increase.  With  argon  we  once  again 
achieved  about  90X  of  the  MaxiMUM  possible  Prandtl  number 
reduction  but  got  approxiMately  6®x  of  the  possible  mass 
increase.  The  results  of  these  selections  will  be  discussed 
in  Chapter  V. 

C.  STACK  DESIGN  THEORY 

The  purpose  of  this  section  is  to  explain  the  diff¬ 
erences  between  the  two  stacks  built  for  this  project  and 
that  used  iu  previous  experiments,  including  the  basis  for 
any  and  all  Modifications.  To  provide  an  orderly  flow  of 
information,  the  differences  between  the  stacks  will  be 
fully  described  and  then  an  explanation  of  each  of  the  new 
features  will  be  given.  Details  of  stack  construction  are 
provided  in  Chapter  III. 

Ue  will  discuss  three  stacks  in  this  section:  the 
original  stack,  built  by  Dr.  Tom  Hofler,  will  be  referred 
to  as  the  "old"  stack,  the  first  stack  built  for  this 
experiment  will  be  referred  to  as  the  "new"  stack,  and  the 
second  stack  built  for  this  experiment  will  be  referred  to 
as  the  "dual"  stack.  The  actual  designs  for  all  of  the 
stacks  are  based  on  the  work  and  ideas  of  Dr.  Hofler.  All 
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of  the  stacks  have  the  same  outer  diameter  of  1.596  inches 
and  the  same  size  center  post  with  a  9.25  inch  diameter. 
Additionally  all  three  have  the  same  plate  thickness  of 
9.993  inches,  and  were  built  using  a  loom,  discussed  in 
Chapter  HI,  with  a  line  spacing  of  9.2  inches.  The  primary 
differences  between  the  stacks  are  in  the  areas  of  average 
plate  spacing,  overall  stack  length,  and  final  arrangement 
of  the  fishing  line  used  to  provide  the  plate  spacing. 

The  old  stack  has  an  overall  length  of  3.989  inches  with 
an  average  plate  spacing  of  9.9149  inches,  with  no 
modifications  to  the  fishing  line  arrangement.  The  new 
stack  has  an  overall  length  of  3.966  inches  with  an  average 
plate  spacing  of  9.9133  inches.  Additionally,  it  has  a 
special  pattern  cut  into  the  fishing  line,  as  shown  in 
Figure  1 1-4.  The  dual  stack  has  an  overall  length  of  3.968 
inches  with  an  average  plate  spacing  of  9.916  inches  at  the 
hot  end  and  9.935  inches  at  the  cold  end.  There  is  no 
pattern  cut  into  the  fishing  line  of  the  dual  stack, 
however,  there  is  a  gap  of  approximately  9,75  inches 
between  the  rows  of  thin  lines  and  thick  lines,  as  shown  in 
Figure  II-5.  Additionally,  on  the  layer  of  plastic  with 
both  sizes  of  fishing  line,  every  other  piece  of  the  thick 
fishing  line  has  been  removed,  starting  at  a  point  18 
inches  from  the  end  of  the  plastic  sheet  that  is  at  the 
center  of  the  roll.  The  only  remaining  difference  between 
the  stacks  is  that  the  old  and  new  stacks  have  center  posts 
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of  a  cloth  phenolic  type  rod  and  the  dual  stack  has  a  nylon 
plastic  type  rod.  This  last  difference  is  strictly  due  to 
availability  of  material  at  the  time  of  construction  and  is 
not  considered  to  be  of  any  significance  in  relation  to 
refrigerator  or  stack  performance. 

Each  of  the  modifications,  relative  to  the  old  stack, 
will  now  be  explained,  beginning  with  those  for  the  new 
stack.  The  changes  in  this  stack  are  primarily  intuitive  in 
nature,  based  on  experiments  by  Hofler  and  Wheatley,  and 
therefore  do  not  lend  themselves  to  quantitative  analysis. 
The  change  in  overall  length  was  selected  to  provide  a 
space  between  the  end  of  the  stack  and  primarily  the  cold 
heat  exchanger.  This  space  is  desired  to  preclude  any  flow 
restrictions  that  could  occur  with  the  stack  plates  resting 
again.' t  the  plates  of  the  heat  exchanger  surface.  A  piece 
of  fishing  line  (approximately  0.020  inches  in  diameter) 
was  glued  to  the  end  of  the  stack,  to  fill  the  void  left  by 
the  reduced  length,  thereby  preventing  the  stack  from 
moving  around  in  the  resonator  relative  to  the  heat 
exchangers.  The  pattern  of  Figure  1 1-4  was  selected  for  two 
reasons.  Since  the  amount  of  heat  transport  is  related  to 
the  total  surface  area  of  the  plates,  the  removal  of  some 
of  the  fishing  line  would  increase  this  surface  area  and 
thereby  increase  the  heat  transport.  Also,  since  it  is 
virtually  impossible  to  have  precisely  uniform  spacing  of 
the  fishing  line  and  plates,  there  are  slight  differences 
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in  the  longitudinal  flow  channels  in  the  stack.  This  in 
turn  can  produce  flow  variations  and  pressure  differences 
from  one  channel  to  the  next.  Due  to  there  being  no 
continuous  pieces  of  fishing  line  along  the  full  length  of 
the  stack,  with  this  pattern,  the  channels  are  all 
interconnected  and  the  pressure  and  flow  differences,  if 
any,  will  be  minimized.  The  overall  result  of  this  pattern 
was  a  reduction  of  approximately  50%  in  the  amount  of 
fishing  line  present  in  the  stack.  This  reduction  should 
also  have  the  effect  of  reducing  the  overall  obstruction  to 
flow  in  the  resonator.  It  should  also  be  noted  that  this 
pattern  may  have  had  a  significant  effect  on  the  degree  of 
crush  experienced  during  the  rolling  process,  resulting  in 
the  smallest  average  plate  spacing  of  any  of  the  stacks. 
This  plate  spacing  is  very  important  in  the  thermoacoustic 
performance  of  the  refrigerator. 

The  changes  in  the  dual  stack  are  somewhat  more  complex 
and  theoretical  in  nature.  Within  the  accuaracy  of 
construction  methods,  the  overall  length  of  the  dual  stack 
is  the  same  as  that  of  the  new  stack,  for  reasons  explained 
above.  The  removal  of  every  other  piece  of  the  thick 
fishing  line,  after  the  first  18  inches,  is  based  on 
reducing  the  overall  amount  of  fishing  line  in  the  stack, 
as  discussed  for  the  new  stack.  All  lines  were  left  intact 
for  the  first  18  inches  to  provide  a  more  uniform  plate 
spacing  toward  the  center  of  the  roll,  where  the  diameter 
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of  the  roll  is  not  very  large  with  respect  to  the  fishing 
line  spacing.  The  0.75  inch  gap  between  the  rows  of  thin 
and  thick  strings  was  used  for  two  reasons.  A  certain 
minimal  gap  is  necessary  to  prevent  the  inner  plastic 
layer,  on  the  hot  end,  from  butting  up  against  the  thick 
strings  when  the  ends  of  the  stack  are  pushed  together.  The 
extra  distance,  out  to  0.75  inches,  was  used  to  give 
approximately  the  same  total  length  of  string  as  we  had  in 
the  new  stack. 

Understanding  the  reasons  for  the  dual  nature  of  the 
stack  requires  some  knowledge  of  the  basic  characteristics 
of  large  temperature  span  refrigerators,  as  well  as  the 
basic  heat  transfer  process  that  occurs  in  the  stack.  The 
heat  flow  at  the  two  ends  of  the  stack  is  not  the  same, 
generally  being  significantly  higher  at  the  hot  end,  where 
heat  is  rejected.  This  difference  is  functionally  dependent 
on  efficiency  and  increases  with  reduced  efficiency.  To 
properly  explain  this  effect  some  terms  must  be  introduced. 
Although  the  following  relations  are  simplified  they  are 
adequate  for  the  explanation  of  this  effect.  The  heat 
transfer  capacity  q  is  given  by: 

q(x)=Q(x)/<  r(x)-l  )=-TT6Kp(x)u(x) ,  (11-52) 

and  the  local  heat  flow  Q  is  given  by: 

Q(x)=-TT5Kp(x)u(x)  (  r(x)-l),  (11-53) 

where  r(x)=vT/vT  vT  _  is  the  iuXiaum  temperature 

cnt  cnt 

gradient  in  the  x  direction,  TT  is  the  surface  area  per  unit 
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length  of  the  stack,  p(x)  is  the  acoustic  pressure 
distribution,  u(x)  is  the  acoustic  velocity  distribution, 
is  the  theraal  penetration  depth  (which  according  to 
Swift  (1989),  is  the  distance  that  heat  can  diffuse  through 
the  fluid  during  the  time  IAj),  and  x  is  the  longitudinal 
distance  along  the  stack.  In  a  standard  stack  u(x)  and  r(x)- 
1  are  the  only  paraaeters  that  vary  significantly  in  the 
longitudinal  direction.  The  quantity  r(x)-l  represents  the 
temperature  gradient  perpendicular  to  the  x  direction  (i.e. 
the  transverse  direction).  This  quantity  will  vary  greatly, 
depending  on  heat  load,  and  aust  be  large  in  order  to 
transfer  large  aaounts  of  heat.  In  a  large  teaperature  span 
refrigerator  |r(x)-l  is  very  small  at  the  cold  end  and 
large  at  the  hot  end.  The  Q(x)  distribution  is  dominated  by 
|r(x)-l|.  Conversely,  u(x)  is  generally  small  at  the  hot 
end  of  the  stack  and  large  at  the  cold  end.  It  is  u(x)  that 
dominates  the  q(x)  distribution.  These  effects  give  a  low 
heat  flow  at  the  cold  end,  where  the  heat  transfer  capacity 
is  highest,  and  a  high  heat  flow  at  the  hot  end,  where  the 
heat  transfer  capacity  is  lowest. 

The  object  of  the  dual  stack  design  is  to  reduce  this 
conflict  by  changing  the  distribution  of  q(x),  making  it 
more  uniform  and  similar  to  that  of  Q(x).  Since  r(x)-l  must 
be  geater  than  zero  to  transfer  heat,  and  is  very  dependent 
on  heat  load,  it  was  determined  that  q(x)  was  the  quantity 
that  should  be  modified.  To  accomplish  this,  the  acoustic 


velocity  diatibution  had  to  be  countered  by  soee  other 
effect.  The  only  parameter  affected  enough  by  stack  design 
to  accomplish  this  is  TT.  By  varying  the  amount  of  surface 
area  along  the  stack  length  the  q(x)  distribution  can  be 
altered.  By  having  half  the  number  of  plates  at  the  cold 
end,  we  cut  the  surface  area  by  a  factor  of  two, 
essentially  countering  the  u(x)  distribution,  and  making 
the  capacity  for  heat  transfer  correspond  more  closely  to 
the  actual  heat  flow  conditions.  Ideally,  the  surface  area 
should  vary  continuously  in  the  longitudinal  direction  to 
precisely  counter  the  velocity  distribution  and  produce  an 
exactly  uniform  q(x)  distribution,  thereby  giving  optimum 
results.  This  situation  would  be  much  more  difficult  to 
construct  and  methods  that  have  been  discussed  to  date  have 
other  drawbacks  that  make  them  less  acceptable  than  the 
current  design.  Although  this  is  the  major  reason  for  an 
expected  increase  in  the  COP  with  the  dual  stack  design 
there  are  two  other  minor  effects  which  should  also  help. 
An  area  not  covered  by  the  simplified  equations  above  is 
that  of  viscous  friction  as  the  gas  moves  along  the  plates. 
The  energy  loss  due  to  this  friction  is  proportional  to 
both  the  square  of  the  acoustic  velocity  and  the  total 
surface  area  exposed  to  the  gas.  Since  the  surface  area  has 
been  cut  in  half  at  the  cold  end,  where  the  velocity  is 
highest,  this  source  of  energy  loss  should  be  significantly 
reduced,  thereby  increasing  the  COP.  The  second  area  of 
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interest  relates  to  the  uniformity  of  the  q(x> 
distribution.  Since  q(x>  depends  on  u(x),  anything  that 
■akes  u(x)  lore  umfora  will  also  sake  q(x )  lore  uniform, 
which  in  turn  tends  to  improve  the  COP.  The  velocity 
distribution  will  depend  somewhat  on  the  degree  of  blockage 
in  the  flow  path.  The  plates  have  an  "effective"  size  that 
is  larger  than  their  thickness,  due  to  viscous  effects.  By 
reducing  the  number  of  plates  the  degree  of  blockage  is 
reduced  and  the  velocity  distlbution  is  more  uniform, 
thereby  tending  to  increase  the  COP. 
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III.  STACK  CONSTRUCTION 


The  theraoacoust ic  stack,  necessary  for  the  transport  of 
heat  as  discussed  m  chapter  il,  consists  of  a  continous 
roll  of  very  thin  plastic,  the  layers  of  which  are 
separated  by  spacers  attached  to  the  plastic  surface.  This 
geonetry  was  selected  instead  of  a  parallel  plate 
arrangeaent  based  on  both  ease  of  construction  and 
structural  rigidity  of  the  finished  product.  When  the 
diaensions  of  a  stack  are  considered,  plate  spacing  of 
9.915  inches  (or  9.933  inches  depending  on  which  stack), 
plate  thickness  of  9.993  inches,  overall  dianeter  of  1.596 
inches,  and  overall  length  of  about  3.966  inches,  it  is 
easily  seen  that  the  parallel  plate  geoaetry  would  be  auch 
aore  difficult  to  construct  and  extrenely  fragile.  The 
plastic  used  in  this  expenaent  was  a  polyester  Hylar  TH 
type.  There  was  no  particular  reason  for  this  selection 
other  than  ease  of  handling.  The  plate  spacers  consist  of 
strips  of  high  grade  aonofilaaent  nylon  fishing  line  which 
are  affixed  to  the  plastic  surface  using  a  spray  adhesive. 

Prior  to  construction,  decisions  aust  be  nade  concerning 
overall  length  of  the  stack,  which  deteriines  the  width  of 
the  plastic  sheet,  required  dianeter,  which  deternines  the 
overall  length  of  the  plastic  sheet  and,  approxiaate  plate 
spacing,  which  deternines  the  size  of  spacers  to  be 
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Attached  to  the  plastic  surface.  For  this  project  two 
stacks  sere  built.  The  first  had  an  overall  length  of  3.e66 
inches,  uith  a  diameter  of  1.596  inches  and  9.915  inch 
spacers.  The  second  was  a  dual  stack  uith  the  sane 
diaseter,  an  overall  length  of  3.968  inches  and  spacers  of 
9.915  inches,  uith  tuo  layers  of  plastic  along  half  of  its 
length,  and  spacers  of  9.933  inches  along  the  other  half. 
This  provided  a  variable  plate  spacing  for  reasons 
discussed  in  chapter  II.  The  overall  length  of  plastic 
sheet  needed  to  make  a  stack  can  be  determined  by 
mathematical  analysis  of  the  spiral  uound  roll  once  the 
dimensions  are  specified.  The  resulting  equation  for  length 
is  L=(  T/d)  (r  where  d  is  the  sun  of  the  spacer  size 
and  plate  thickness,  r j  is  the  radius  of  the  stack  (and 
also  the  radius  of  the  resonator  bore),  and  r^  is  the 
radius  of  the  center  post,  on  which  the  plastic  is  rolled. 
The  center  posts  that  have  been  used  are  9.25  inch  diameter 
•and  are  either  cloth  phenolic  or  a  nylon  type  plastic  rods. 
Usually  a  safety  margin  of  twelve  to  eighteen  inches  is 
added  to  the  length  calculated  by  the  above  equation.  This 
is  due  to  the  fact  that  a  precise  fit  is  desired  in  the 
resonator  bore  and  it  is  easier  to  trim  small  pieces  from 
an  oversized  roll  than  to  add  pieces  to  one  that  is 
undersized.  The  width  of  the  plastic  sheet  is  cut  very 
precisely  using  a  sharp  paper  cutter,  uith  a  guide  that  has 
been  carefully  positioned  to  the  desired  width.  This  method 
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has  provided  results  with  an  accuracy  of  9.902  inches.  The 
overall  length  of  plastic  sheet  used  for  the  two  stacks  in 
this  experiaent  was  129  inches  for  the  first  stack  and  two 
pieces  of  60  inches  each  for  the  dual  stack. 

Once  the  plastic  sheet  has  been  cut  to  size  and  fishing 
line  of  the  required  size  has  been  acquired,  construction 
can  begin.  The  spray  adhesive  used  in  this  experiaent  was 
3H  brand  Super  77  spray  adhesive.  Since  attaching 
individual  pieces  of  fishing  line  would  be  very  tiae 
consuaing  and  result  in  a  less  than  unifora  spacing  an 
alignaent  device  is  used  to  assist  in  this  process.  An 
aluainua  looa  has  been  constructed  to  provide  uni  fora 
parallel  spacer  positioning  on  the  plastic  surface.  A  long 
piece  of  fishing  line  is  wound  around  the  hooks  of  the 
looa,  positioned  in  grooves  to  provide  parallell isa,  and 
fixed  at  both  ends  with  screws  to  hold  it  in  place.  The 
looa  has  an  open  front  and  back  to  provide  access  to  both 
sides  of  the  fishing  line.  Using  a  spray  shield  to  protect 
the  looa,  spray  adhesive  is  applied  to  the  front  side  of 
the  fishing  line.  The  looa  is  then  turned  over  and  aligned 
over  the  plastic  to  provide  parallel,  unifora  positioning 
of  the  strips  of  fishing  line.  Once  positioned  the  looa  is 
set  down  and  a  saooth  aluainua  block  is  used  to  press  the 
fishing  line  onto  the  surface.  The  working  surface  and 
press  used  in  this  experiaent  were  saoothly  polished 
aluainua  blocks,  which  were  necessary  to  provide  a 
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relatively  uniform  compression  of  the  fishing  line  when 
pressing.  A  teflon  sheet  urns  pieced  under  the  plestic  end 
over  the  working  surfece  to  prevent  the  strings,  which 
extend  beyond  the  plestic,  from  sticking  to  the  working 
surfece.  A  sketch  of  the  loom,  with  fishing  line  etteched 
is  shown  positioned  over  the  plastic  sheet  in  Figure  II 1-1, 
The  distance  between  strips  of  fishing  line  on  the  loom 
used  m  this  experiment  was  9.20  inches.  This  spacing  was 
selected  based  on  experiments  conducted  by  Dr.  Tom  Hofler 
at  the  Naval  Postgraduate  School.  The  loom  covers  a  length 
of  six  inches  and  therefore  requires  twenty  repetitions  of 
this  procedure  to  complete  a  stack.  The  dual  stack  required 
sore  repititions  due  to  there  being  two  different  string 
sizes  on  one  plastic  sheeet  and  a  second  layer  of  plastic 
of  half  the  width,  as  shown  in  Figure  I I 1-2.  Once  the 
fishing  line  is  pressed  down  and  allowed  a  minimal  drying 
time  of  one  to  two  minutes  a  scalpel  or  razor  blade  is  used 
to  cut  the  fishing  line  off  at  the  edge  of  the  plastic.  It 
should  be  noted  that  if  the  line  is  wound  on  the  loom  with 
too  much  tension  the  plastic  will  tend  to  curl  when  the 
line  is  cut.  This  is  obviously  an  undesireable  situation  so 
care  must  be  taken  when  winding  the  line  on  the  loom  to 
provide  just  enough  tension  to  maintain  the  strips  in  the 
grooves.  After  all  of  the  fishing  line  is  attached  to  the 
plastic  some  trimming  must  be  done  prior  to  rolling.  The 
fishing  line  is  carefully  cut  back  approximately  one 
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sixteenth  of  an  inch  from  each  edge,  to  prevent  potential 
flow  problems  at  the  heat  exchanger  interfaces. 
Additionally,  depending  on  the  stack  design,  a  pattern  is 
cut  in  the  strings  or  some  strings  are  removed,  for  reasons 
explained  in  the  stack  theory  of  chapter  II,  Great  care 
should  be  exercised  in  the  triaaing  of  any  patterns  to 
avoid  Baking  any  large  cuts  in  the  plastic  surface  as  this 
could  produce  significant  flow  variations  during  operation 
of  the  unit. 

After  completing  this  process  the  stack  is  ready  to 
roll.  The  center  post  is  carefully  taped,  using  scotch 
tape,  to  the  end  of  the  plastic  3heet.  The  plastic  is 
carefully  rolled  up,  using  just  enough  tension  to  maintain 
the  shape  and  provide  uniforn  layer  spacing,  without 
getting  tor.  much  crush,  which  would  result  in  a  layer 
separation  much  smaller  than  the  string  diameter.  This  is 
in  general  the  most  difficult  part  of  the  stack 
construction  process.  The  uniform  layer  spacing  is  checked 
visually  by  holding  the  stack  up  to  a  light  or  using  some 
sort  of  magnifying  device.  The  spacing  can  also  be  checked 
quantitatively  based  on  the  known  innner  post  diameter, 
the  measured  outer  diameter  and  a  count  of  the  number  of 
layers  along  several  radii.  This  will  provide  an  average 
layer  spacing  which  will  in  general  be  somewhat  smaller 
than  the  sum  of  the  fishing  line  diameter  and  the  plastic 
thickness.  This  is  due  to  a  certain  degree  of  crush  or 
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compression  during  the  pressing  and  rolling  processes. 
After  rolling,  even  with  a  great  deal  of  care,  there  wi 11 
he  some  degree  of  spreading,  or  misalignment,  at  the  ends 
of  the  cylinder.  Using  tuo  smooth  surfaces  the  ends  of  the 
stack  can  be  pressed  and  the  layers  will  slide  into 
position  to  give  flat  ends  and  an  overall  length  as  precise 
as  the  width  of  the  plastic  sheet. 

Drawings  of  the  two  ends  of  the  dual  stack  are  shown  in 
Figures  1 1 1—3  and  1 1 1-4.  The  other  stack  has  the  same 
appearance  as  the  dense  or  hot  end  of  Figure  1 1 1-3. 
Polaroid  photographs  of  the  actual  dual  stack,  showing  a 
radius  from  the  center  post  out,  at  a  6.4x  magnification, 
are  provided  as  Figures  1 11-5  and  1 1 1-6.  The  finished  roll 

is  then  trimmed  down  in  diameter,  by  cutting  small  pieces 

'  |!  ,  ,1  / 

off  the  length  of  the  plastic  strip,  to  fit  snugly  aii  the 

resonator  bore.  The  final  product  is  taped  along  the  seam, 

•  *• 

with  scotch  tape,  to  hold  the  roll  in  the  cylindrical 


shape. 


v.y- 


Although  some  of  the  differences  between  the  two  stacks 
that  were  built  for  this  experiment  are  mentioned  in  here, 
the  details  of  and  reasons  for  the  differences  are 
presented  in  the  stack  theory  of  chapter  II. 


ure  MI-3  Looking  downvpn  to  -  Figure  I  I  I  "4  Looking  Down  on 

Hot  End  of  the  Dual  Stack  --  the  Cold  End  of  the  Dual  Stack 


Figure  III — S.  Photograph  of  Dual  Stack, 
End,  at  6.4x  Magnification 
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Figure  1 1 1-6.  Photograph  of  Dual  Stack,  Hot 
End,  at  6.4x  Magnification 
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IV.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

In  this  Chapter  we  will  present  a  description  of  our 
experimental  apparatus,  including  the  electronics  used,  and 
a  discussion  of  our  experimental  procedure.  We  have  used 
the  prototype  thermoacoustic  refrigerator,  and  associated 
electronics  package,  that  was  built  by  Hofler  (1986,1988). 
We  have  paraphrased  sections  of  his  discussion  of  the 
system  as  necessary  to  provide  a  basic  understanding  of  the 
apparatus  and  its  operation. 

A.  DESCRIPTION  OF  APPARATUS 

In  order  to  properly  describe  the  apparatus  it  is 
necessary  to  break  it  down  into  two  major  areas;  that  of 
the  high  intensity  driver,  with  its  instrumentation  for 
acoustic  power  measurements  and  that  of  the  refrigeration 
components  (resonator,  stack,  heat  exchanger,  etc.). 

1 .  The  Driver  and  its  Instrumentation 

The  driver  consists  of  a  modified  commercial 
loudspeaker,  shown  in  Figure  I V  —  1 .  The  modifications 
consist  of  removal  of  the  original  plastic  housing  and 
fabric  dome  and  attachment  of  a  thin-wall  aluminum  cone- 
shaped  piston  and  a  surround.  Figure  I V—  1  shows  the  unit 
with  a  voice  coil  and  former  1,  iron  pole  pieces  2,  magnet 
3,  aluminum  piston  4,  surround  5,  aluminum  pressure  vessel 


Figure  IV- 1.  Scheaatic  of  Driver  Apparatus  Showing 
Acoustic  Hotor,  Transducers,  and  a  Portion  of  a 
Resonator.  Numbered  Parts  are  Identified  in  the  Text. 
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6  and  7,  aluminum  plate  8,  cooling  water  tubes  9,  annular 
plate  (where  pressure  transducer  is  mounted)  1®,  dynamic 
pressure  transducer  11  and  a  post  for  the  mounting  of  an 
accelerometer  12.  The  unit  has  a  two  inch  diameter  voice 
coil  and  an  effective  piston  driving  area  of  ®. 93  square 
inches.  Additionally,  it  has  a  high  power  capability  that 
is  made  possible  by  suspending  ferrofluid  in  the  magnet 
gap,  to  conduct  heat  from  the  voice  coil.  A  capillary 
filled  hole  is  located  adjacent  to  the  surround  to  break 
the  gas  seal  provided  by  the  surround.  The  flow  impedance 
of  the  capillary  was  chosen  to  be  as  low  as  possible 
without  affecting  the  acoustic  power  measurements.  The 
assembly  is  housed  in  the  aluminum  pressure  vessel  6  and  7, 
and  the  resonator  section  is  bolted  to  6.  The  heat  removed 
by  the  water  cooling  tubes  9,  is  conducted  through  the 
aluminum  housing  6  and  plate  8. 

The  dynamic  pressure  transducer  11,  is  physically 
attached  to  the  structure  of  the  driver  and  can  therefore 
be  calibrated  along  with  other  aspects  of  the  driver, 
independent  of  the  resonator.  The  transducer  itself  was 
designed  and  built  by  Dr.  Hofler.  It  consists  of  a  thin  y- 
cut  quarts  crystal  disk  with  evaporated  electrodes  and  a 
miniature  FET  impedance  converter  circuit  that  buffers  the 
output  signal.  The  dynamic  pressure  transducer  provides  a 
low  noise,  low  distortion,  phase  accurate  signal  which  is 
calibrated  to  determine  the  volume  velocity  in  a  small 


82 


calibration  cavity  that  can  be  attached  to  the  driver. 
This,  in  conjunction  with  the  integrated  accelerometer 
signal  (which  is  then  a  velocity  signal)  and  the  phase 
difference  between  the  two,  provides  for  an  accurate 
measurement  of  the  acoustic  power  provided  to  the 
resonator.  HoflerVs  calibration  of  the  pressure  transducer 
was  accomplished  by  comparison  to  a  “standard"  transducer 
that  had  previously  been  calibrated  to  an  accuaracy  of  ±  1% , 
by  Swift,  et .  al.  (1382),  resulting  in  a  sensitivity  of 
1.292  (tlX)  volts/bar.  The  accuracy  of  the  acoustic  power 
measurement  is  ±3%,  as  shown  by  Hofler  (1988). 

As  previously  noted  in  the  discussion  of  Figure  IV- 
1,  the  accelerometer  is  mounted  on  the  post  12,  which  is 
attached  to  the  back  side  of  the  piston  4.  The  post  is 
necessary  to  thermally  isolate  the  accelerometer  from  the 
piston  (which  becomes  warm  during  operation)  and  to 
position  the  accelerometer  where  its  cable  has  clearance 
from  the  piston  .  wall.  The  accelerometer  is  a  shear 
piezoceramic  type.  The  piston’s  velocity  signal  (which  is 
the  integrated  output  of  the  accelerometer)  is  calibrated 

with  the  dynamic  pressure  transducer  signal  using  a  small 

• 

volume  of  known  acoustic  impedance,  to  get  an  accurate 
value  for  volume  velocity.  This  calibration  may  change  as 
the  properties  of  the  suround  change  with  use  and  must 
therefore  be  repeated  at  periodic  intervals,  to  maintain 
the  accuracy  of  results.  Repeated  calibrations,  over  a  two 


83 


year  period,  have  shown  that  the  sensitivity  of  the 
pressure  transducer  has  been  stable  within  about  0. IX, 
while  the  volume  velocity  calibration  value  has  shown  an 
overall  decrease  of  about  ?X.  It  should  be  noted  that  the 
voluae  velocity  value  dopped  just  over  4X  during  the  first 
four  months  of  driver  use,  and  dropped  only  3x  wore  over 
the  next  18  months. 

2.  Refrigeration  Components 

The  main  components  of  the  refrigeration  section  are 
shown  in  Figure  IV-2.  The  resonator  (18,  20  and  21  in  the 
figure)  is  constructed  of  several  parts.  Where  temperature 
uniformity  is  desired  heavy-gauge  copper  is  used.  In 
section  18,  containing  the  stack  17,  there  is  a  large 
temperature  gradient  and  an  epoxy-fiberglass  substance  is 
used  to  provide  high  strength  with  very  low  thermal 
conductance.  This  epoxy-fiberglass  section,  by  itself  would 
be  very  permeable  to  gas  diffusion  if  left  alone.  To 
correct  this  situation.  Dr.  Hofler  evaporated  a  thin  film 
of  indium  over  the  outside  of  the  fiberglass  section,  and 
then  wrapped  more  fiberglass  over  the  metal  film. 

The  primary  refrigeration  component  is  the  stack  (17 
in  the  Figure  IV-2)  which  is  constructed  of  a  long  sheex  of 
plastic  that  is  spiral  wound  on  a  0.25  inch  diameter 
plastic  or  cloth  phenolic  rod  16.  The  layers  of  plastic  are 
separated  by  monofilament  fishing  line  spacers.  The  details 
of  stack  construction  were  presented  in  Chapter  III. 
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The  other  components  of  importance  in  the 
refrigeration  section  are  the  hot  and  cold  heat  exchangers 
15  and  22  in  Figuure  IV-2,  respectively.  The  hot  heat 
exchanger  provides  for  removal  of  heat  from  the  unit,  to 
prevent  overall  heating  of  the  stack  during  operation.  The 
cold  heat  exchanger  provides  good  thermal  contact  uith  the 
applied  heat  load,  from  the  heater  coil  19.  Dr.  Holler 
constructed  the  heat  exchangers  using  electroplating  and  a 
chemically  removable  aluminum  form,  shown  in  Figure  IV-3 
prior  to  removal  of  the  aluminum.  The  resulting  heat 
exchangers  had  parallel  copper  strips  with  the  planes  of 
the  strips  parallel  to  the  tube  axis.  The  details  of  this 
construction  are  provided  in  the  doctoral  dissertation  of 
Hof ler  (1986). 

The  remainder  of  the  physical  apparatus  consists  of 

a  vacuum  system,  that  attaches  below  the  driver  housing, 

thereby  enclosing  the  entire  refrigeration  section,  and  a 

section  of  super-insulation  wrapped  around  the  resonator. 

All  of  this  is  to  provide  thermal  insulation,  thereby 

reducing  nuisance  heat  loads,  A  layer  of  super-insulation 

consists  of  a  layer  of  glass  micro-fiber  web  followed  by  a 

layer  of  aluminized  polyester  film.  The  entire  resonator 

section  is  wrapped  with  10  to  12  of  these  layers.  The 

vacuum  system  uses  an  ion  gauge  to  monitor  the  vacuum, 

-4  -6 

which  varies  between  3x10  and  5x10  torr  during 
operation,  depending  on  operating  conditions  and  gas  types. 


B.  DESCRIPTION  OF  ELECTRONICS 


The  electronics  portion  of  our  experimental  apparatus 
consists  of  two  primary  sections*  These  are  the  control  and 
protective  section  and  the  data  monitoring  section.  The 
control  and  protective  systems  will  be  discussed  here  and 
the  data  monitoring  aspects  will  be  covered  in  the  nex\ 
section,  where  ue  discuss  the  experimental  procedure. 

The  control  and  protective  electronics  system  provides 
feedback  control  of  the  drive  frequency  and  amplitude, 
protection  against  burn-up  of  the  driver  voice  coil,  and 
signal  conditioning  to  increase  the  precision  of  acoustic 
power  measurements.  The  protective  circuit  is  particularly 
important  in  our  situation  because  the  required  operating 
conditions  are  very  close  to  the  maximum  limits  of  the 
driver. 


Figure 

I V-4 

is  a 

block  diagram  of  the  control  circuits 

for  drive 

frequency 

and  amplitude. 

This  figure 

iS 

reproduced, 

with 

minor 

equipment  updates, 

from  the  doctoral 

dissertat  ion 

of 

Tom 

Hof ler  (1986).  The 

explanation  of 

its 

operation , 

for 

the 

most  part,  is 

quoted  here 

for 

comp  1 eteness : 


"The  frequency  control  is  based  on  a  lock-in  amplifier 
that  is  used  as  a  phase  null  detector  uith  an  analog 
output.  The  velocity  signal  is  fed  into  the  lock-in 
preamp,  and  the  pressure  signal  is  used  as  a  reference 
channel  input.  The  level  of  this  reference  signal  is 
adjusted  to  bj  IV  rms,  in  order  to  minimize  lock-in 
phase  errors.  The  Y-channel  output  voltage  of  the  lock-m 
is  proportional  to  the  measured  phase  deviation  and  is 
fed  into  a  simple  integrating  controller.  The  controller 
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output  is  connected  to  the  voltage-controlled-frequency 
(VCF)  input  of  the  signal  generator,  and  the  loop  is 
closed  by  the  driver-resonator  system.  The  phase, 
including  calibration  corrections,  is  set  on  the  lock-in 
so  that  the  desired  Y-channel  output  voltage  is  zero  when 
the  pressure  and  velocity  are  in  phase.  The  controller 
has  a  set  point  voltage  fixed  at  zero,  so  the  frequency 
is  controlled  to  maintain  the  resonance  condition  of  the 
pressure  and  velocity  being  in  phase.  The  pressure 
amplitude  control  group  begins  by  adjusting  the  amplitude 
of  the  pressure  signal  to  a  pre-defined  level  with  an 
adjustable  gain  stage.  This  ensures  that  the  loop  gain 
and  stability  margin  are  constant  for  different  drive 
levels.  After  the  pressure  signal  is  adjusted  in 
amplitude,  it  is  rectified  and  filtered,  and  this  DC 
signal  is  used  by  another  integrating  controller  to 
produce  a  control  voltage  for  an  adjustable  gain  control 
circuit  (AGO.  The  output  of  the  signal  generator  is  then 
passed  through  the  AGC  and  into  the  power  amplifier 
input.  The  loop  is  again  closed  by  the  acoustic  system, 
with  its  pressure  amplitude  thereby  being  regulated. 

Integrating  the  accelerometer  signal  is  convenient 
since  it  is  the  velocity  that  is  the  important  quantity 
in  a  power  measurement...." 

To  finish  the  explanation  we  will  need  to  discuss  the 
velocity  signal  measurement  since  the  method  used  by  Hofler 
was  not  used  in  our  experiment.  In  our  experiment  ue 
measured  the  velocity  amplitude  using  the  lock-in  directly. 
It  was  determined  that  this  would  not  introduce  any  errors 
into  the  power  measurements,  as  previously  thought.  Unlike 
the  previous  apparatus,  our  lock-in  amplifier  uses  true 
sine  wave  conversion,  resulting  in  complete  rejection  of 
the  harmonic  distortion  in  the  accelerometer  signal.  Using 
the  measured  dynamic  pressure  and  velocity,  along  with  the 
sensitivities  of  the  transducer  and  accelerometer,  and  the 
volume  velocity  calibration  factor  we  can  calculate  an 
accurate  acoustic  power,  as  long  as  the  frequency  control 
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circuit  is  maintaining  the  phase  difference  between  the  two 
signals  at  zero. 

The  power  protection  circuit  monitors  the  power  amp¬ 
lifier  current  and  will  trigger  a  crowbar  circuit  to  blow  a 
fuse  if  either  of  two  setpoints  are  exceeded.  If  a  set  peak 
current  value  is  exceeded  the  fuse  is  blown,  to  protect 
against  large  fast  transients.  If  a  set  rms  current  value 
is  exceeded  the  fuse  will  also  be  blown,  to  protect  against 
currents  below  the  peak  limit  but  high  enough  to  burn  out 
the  voice  coil. 

C.  EXPERIMENTAL  PROCEDURE 

In  this  section  we  will  discuss  three  basic  areas  of  our 
experimental  procedure.  First  is  the  setup  of  the  apparatus 
for  a  specific  data  run.  Then  we  will  cover  the  procedure 


for  starting 

up 

and  operating 

the 

refrigerator . 

Ue 

Will 

conclude 

with 

a 

discussion 

of 

the 

actual  data 

that  is 

recorded. 

In 

the 

operating 

section 

we  will  cover 

two 

different  methods  since  the  procedure  was  changed  after  the 
second  data  run,  for  reasons  that  will  be  explained  later. 

Before  discussing  the  actual  setup  it.  should  be  noted 
that  all  temperatures  that  we  recorded  were  measured  using 
commercial  type-E  thermocouple  wire  with  an  accuracy  of 
i  1 . 7  *C  or  i  1%  of  the  reading,  whichever  is  greater.  All  of 
the  thermocouple  wires  are  fed  through  a  Hewlett  Packard 
HP3421A  data  acquisition  control  unit  to  a  computer,  for 


continuous  display  on  the  computer  monitor  during 
operation . 

1 .  Setup  of  Apparatus 

The  setup  of  the  system  for  a  particular  data  run 
consists  of  filling  -he  system  with  the  desired  gas, 
pressurizing  the  system  up  to  the  desired  pressure  of  10 
bar  and  starting  up  tne  external  vacuum  system. 

The  resonator  is  filled  through  a  connection  at  the 
top  of  the  driver  housing,  however  the  controlling  factor 
is  the  capillary  near  the  driver  surround  that  was 
previously  discussed.  The  pressurization  rate  is  monitored 
on  a  differential  pressure  gauge  since  the  surround  is 
capable  of  sustaining  no  more  than  0.25  bar  of  differential 
pressure  without  permanent  damage.  The  system  is  filled 
using  a  process  of  evacuation  and  purging  that  is  repeated 
twice  before  the  system  is  brought  up  to  operating 
pressure.  This  procedure  is  designed  to  remove  any  gas 
other  than  the  desired  gas  (particularly  air  if  the  system 
has  been  opened)  prior  to  system  operation.  After  the 
second  evacuation  the  system  is  pressurized  to  10  bar.  The 
evacuation  and  purging  process  consumes  the  better  part  of 
a  day  due  to  the  low  flow  rates  that  are  possible  within 
the  allowed  differential  pressure  range. 

To  save  time,  the  external  vacuum  system  can  be 
started  in  parallel  with  the  above  procedure.  Since  the 
diffusion  pump  of  the  vacuum  system  creates  a  great  deal  of 
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heat,  an  external  cooling  supply  is  used  to  keep  it  at  a 
reasonable  temperature.  Therefore,  to  start  up  the  vacuum 
system  the  cooling  system  is  lit  off,  the  main  vacuum  pump 
sucticn  is  closed,  the  vacuum  pump  is  started  and  the  pump 
suction  is  slowly  opened,  to  prevent  rapid 
depressurization  and  damage  to  the  previously  discussed 
super-insulation.  Then  the  diffusion  pump  is  turned  on. 

Once  the  system  is  up  to  pressure  and  the  vacuum  is  below 

-3 

1x10  torr  it  is  ready  to  start  up. 

2.  System  Startup  and  Operation 

To  start  the  refrigerator  the  two  controllers 
(frequency  and  pressure  amplitude)  are  turned  off  and  the 
signal  generator  is  set  up  for  a  sine  wave  at  the 
approximate  resonant  frequency  for  the  resonator  and  gas 
type  being  used.  The  drive  amplitude  is  set  to  zero  and  the 
power  is  turned  on.  The  amplitude  is  brought  up  slowly 
until  a  drive  current  of  about  0.5  amps  is  achieved.  At 
this  point  the  frequency  is  manually  adjusted  to  resonance 
(i.e.  zero  quadrature  output)  and  the  frequency  control 
circuit  is  turned  on.  Then  the  amplitude  is  slowly 
increased  to  the  desired  value  (if  possible)  and  the 
reference  channel  is  adjusted  to  precisely  1.00  volts.  Then 
the  pressure  amplitude  control  circuit  is  turned  on. 
Depending  on  the  desired  pressure  amplitude  and  gas  type  it 
is  not  always  possible  to  reach  the  desired  operating  point 
at  initial  startup.  It  is  therefore  necessary  to  select 
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some  lower  amplitude,  until  the  engine  cools  down  and  the 
desired  amplitude  can  be  reached.  The  limiting  factor  is 
generally  the  piston's  displacement  amplitude.  If  the 
maximum  displacement  amplitude  is  exceeded  significant 
damage  to  the  surround  occurs.  Once  the  desired  amplitude 
is  achieved  the  refrigerator  runs  unattended  until  steady 
state  is  reached.  In  this  case,  steady  state  is  defined  as 
less  than  a  one  to  two  percent  deviation  in  the  cold  end 
temperature  over  a  90  minute  period.  This  is  determined  by 
evaluating  the  continuous  computer  monitor  display  of 
temperatures,  which  shows  approximately  90  minutes  of  data 
history  when  set  up  for  display  at  four  minute  intervals. 
This  is  the  basic  procedure  for  startup  and  setup  for  a  no 
(applied)  heat  load  data  point.  In  general  equilibration 
takes  from  8  to  24  hours  depending  on  the  pressure 
amplitude  setting. 


We 

have 

used 

two  different  methods 

for  obtaining 

data 

with 

heat 

loads 

applied  to  the  cold  end.  The  first 

method 

is 

that 

used 

by  Hofler  (1986),  where 

a  designated 

heat 

load 

is 

appl ied 

to  the  heater,  at  the 

cold  end,  and 

equ  i  1  i 

-brat i on 

is  based  on  temperature  stab 

ilisation,  as 

d i scussed 

above 

for  the 

no  load  case.  This  heat  load  is  set 

by  applying  a 

fixed 

DC  current  and  voltage 

to  the  heater 

c  i  1 . 

It 

was 

determined  after  two  data  runs 

that  a  better 

method 

might  be 

available.  This  better  method, 

designed  and 

set  up  by  Dr.  Hofler,  invloves  a  temperature  feedback 
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control  network  'Mch  allows  setting  a  desired  steady  state 
cold  end  temperature.  The  controller  then  establishes  the 
current  flow  to  the  heater  that  is  necessary  to  achieve  and 
maintain  that  temperature.  The  equilibration  point  is  no 
longer  based  on  temperature  but  is  instead  based  on  the 
heater  current  deviation  being  less  than  approximately 
:0.5mA  over  a  90  minute  period.  This  corresponds  to  being 
within  better  than  tlX  over  the  90  minute  period,  for  the 
currents  ue  have  measured.  The  control  network  involves  a 
rather  complex  Fortran  computer  program,  written  by  Dr. 
Hofler,  that  reads  the  heater  thermocouple,  at  an  operator 
selected  time  interval,  and  then  performs  a  crude  time 
integral  of  the  difference  between  the  monitored 
temperature  and  the  selected  set  point  temperature.  The 
program’s  integrator  output  variable  is  used  to  numerically 
control  the  heater  current  from  the  HP6633A  system  D.C. 
power  supply.  The  program  allows  the  operator  to  select  the 
control  time  interval  (usually  5  seconds),  the  periodicity 
of  display  (usually  4  minutes),  the  desired  temperature,  a 
proportional  constant  in  volts/deg  (usually  0.0),  and  an 
integration  constant  in  vo 1 ts/deg-sec  (usually  0.01).  The 
latter  two  of  these  control  the  stabilization  time  and  the 
degree  of  stabilization.  The  usual  values  that  were  given 
are  those  that  have  shown  the  best  results  during  our 
exper i ments . 
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This  temperature  feedback  control  method  was 
determined  to  be  better  l'or  the  STAR  project  for  two  basic 
reasons.  First,  it  provided  a  more  rapid  equilibration 
because  the  refrigerator  was  loaded  more  heavily  at  first, 
to  drive  the  temperature  of  the  cold  end  to  the  desired 
temperature,  and  then  the  load  was  let  off  only  as 
necessary  to  maintain  that  temperature.  This  reduced  the 
equilibration  time  to  6  to  8  hours  in  most  cases.  Secondly, 
this  method  conserves  po-;er  due  to  the  reduced  equil¬ 
ibration  time.  This  is  considered  to  be  vital  to  the  STAR 
project  due  to  its  reliance  on  batteries  for  power,  which 
gives  it  a  limited  operating  lifetime. 

3.  Data  Recording 

Onco  steady  state  has  been  established  a  data  set  is 
recorded  in  a  notebook.  A  data  set  consists  of  nine 
readings,  as  well  as  date  and  time.  First,  the  temperatures 
of  the  hot  and  cold  heat  exchangers  are  recorded.  There  are 
two  thermocouples  for  each  heat  exchanger,  one  being  at  the 
center  and  one  at  the  wall.  Due  to  an  intermittent  failure 
of  the  thermocouple  at  the  center  of  the  cold  heat 
exchanger  we  have  only  recorded  the  other  three 
temperatures.  Then  we  record  the  drive  frequency,  the  mean 
pressure,  the  amplitude  of  the  dynamic  pressure  signal  and 
the  amplitude  of  the  lock-in  velocity  signal.  Then, 
depending  on  which  method  of  operation  was  used,  ue  record 
the  DC  heater  voltage  and  current,  or  just  the  DC  heater 


96 


current.  For  the  letter  case  we  use  the  known  heater 
resistance  versus  temperature  character ist ic  to  calculate 
the  heater  load  (instead  of  voltage  times  current).  After 
the  data  are  recorded  we  change  either  the  dynamic  pressure 
amplitude  or  the  heat  load,  and  wait  for  the  system  to 
equilibrate.  A  standard  data  run  for  a  given  gas  type  and 
stack  involves  8  to  Id  data  points  and  takes  about  a  ueek, 
from  start  to  finish.  The  difference  in  the  number  of  data 
points  is  based  on  the  different  temperatures  achievable  at 
different  pressure  amplitudes,  while  staying  within  the 
piston’s  displacement  amplitude  limit. 

In  addition  to  the  above  data,  we  also  record  the 
vacuum  in  the  system  and  the  temperature  and  resonant 
frequency  of  a  gas  analysis  tube,  similar  to  that  described 


by 

Polturak,  et. 

al . 

( 1986) ,  and 

Garrett,  et.  al.  (1986), 

that 

we  attached 

to 

the 

system. 

These  are  monitored  just 

for 

completeness 

and 

are 

not  directly  involved  in  any 

performance  evaluation.  With  the  exception  of  one  data  run, 
all  of  our  gases  had  analysis  certifications  for  purity, 
and  these  values  were  used  for  calculations.  In  one  case  ue 
mixed  our  own  helium-argon  gas  and  used  the  gas  analyser, 
built  by  D  .  Hofler,  to  determine  the  mixture 
concentration . 

The  final  bit  of  data  that  is  recorded  involves  the 
rate  at  which  the  engine  warms  up,  after  being  fully  cooled 
down  and  turned  off.  This  data  is  analyzed  and  together 
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with  the  total  heat  capacity  of  tha  cold  and  of  tha  systea, 
is  usad  to  detarama  a  valua  for  tha  haat  laak  out  of  tha 
refrigerator.  Tha  haat  laak  datarainad  by  this  aathod  is 
than  raducad  by  2.4  aw/*c,  which  is  tha  portion  that  has 
baan  datarainad  to  ba  dua  to  tha  tharaal  conductance  of  tha 
stack  and  gas,  which  is  internal  and  intrinsic  to  tha 
stack.  This  haat  laak  is  than  included  in  tha  total  haat 


load. 


v.  pm 

In  this  chapter  we  Mill  discuss  the  Methods  end  equa¬ 
tions  that  Mere  used  m  analyzing  our  data,  this  Mill  be 
followed  by  a  nuaerical  and  graphical  presentation  of 
various  aspects  of  our  experieent.  It  should  be  kept  in 
mind  that  the  purpose  of  these  expennents  was  to  iaprove 
the  coefficient  of  perforeance  (COP)  of  the  prototype 
refrigerator,  built  by  Hofler  (1986). 

A.  DATA  ANALYSIS 

All  of  our  data  analysis  was  done  through  the  use  of  a 
Microsoft  Multiplan  spreadsheet  prograe  and  graphical 
display  of  the  nuabers  that  were  thereby  generated.  This 
Method  was  used  because  mo  had  voluemous  anounts  of  data 
where  Many  of  the  required  calculations  were  repetitive  in 
nature.  Taking  into  account  all  of  the  quantities,  both 
recorded  and  calculated,  we  ended  up  with  a  spreadsheet  of 
25  coluans  by  122  rows.  It  should  be  noted  that  Many  of 
these  quantities  were  calculated  for  inforwat ional  purposes 
and  are  not  directly  related  to  our  final  evaluation  of  the 
refrigerator's  coefficient  of  perforeance. 

Our  analysis  consists  of  coeparing  a  norealized 
(norealized  to  the  Carnot  COP  for  the  sane  teeperature 
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difference)  coefficient  of  performance,  for  the  refrig¬ 
erator  with  modif ication*  (COM) ,  to  the  same  quantity  for 
the  original  Hofler  refrigerator.  Since  we  have  Modified 
both  the  thermoacoustic  stack  and  the  gas  type,  tie  have 
three  basic  areas  for  comparison.  He  Mill  compare  the 
results  with  different  stacks  and  the  same  gas  type,  with 
different  gases  and  the  same  stack,  and  then  tilth  different 
|  stacks  and  gases  combined.  In  order  to  make  these 

comparisons  tie  have  created  plots  of  COPR  versus  total  heat 

i 

load  (Qtotal).  In  addition  to  these  we  have  created  plots 

I 

of  temperature  ratio  (Tc/Th)  versus  Qtotal,  to  validate  our 

I 

j  data.  Based  on  the  experiments  of  Hofler  (1986),  these 

i 

latter  plots  are  expected  to  be  linear. 

In  order  to  create  these  plots  we  need  to  know  the 
values  for  COPR,  Qtotal,  Tc  and  Th.  The  values  for  Tc  and 
Th  are  simply  raw  data  taken  during  the  experiments.  We 
have  used  the  temperature  at  the  center  of  the  hot  heat 
exchanger  as  Th,  and  the  temperature  at  the  wall  of  the 
cold  heat  exchanger  as  Tc  (due  to  the  intermittent  failure 
I  of  the  cold  heat  exchanger's  center  thermocouple). 

I 

S  The  value  for  Qtotal  has  two  components.  The  first  is 

the  column  labelled  QHEAT  in  Appendix  C.  This  is  the  power 
applied  to  the  heater,  located  below  the  cold  heat 

i 

[exchanger,  in  watts.  It  is  calculated  by  either  multiplying 
the  heater  DC  voltage  times  the  heater  DC  current,  or  by 
i  multiplying  the  heater  resistance  by  the  square  of  the 

i 

» 

s 
i 

f 
? 
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heater  DC  current  (P*VI  or  P=* I  * R  for  DC).  The  voltage  and 
current  values  are  raw  data.  The  resistance  values  were 
taken  from  a  plot  of  Tc  vs.  R  created  from  the  two  sets  of 
data  where  we  recorded  both  heater  voltage  and  current 
(R=V/I).  This  plot  is  shown  as  Figure  V-l.  Both  of  these 
methods  were  used  as  discussed  in  Chapter  IV.  The  second 
component  of  Qtotal  is  the  heat  leak,  which  is  labelled 
HEATLK  in  Appendix  C.  As  discussed  in  Chapter  IV,  this 
value  is  determined  by  analysis  of  the  rate  at  which  the 
refrigerator  warms  up  once  the  power  is  turned  off.  Based 
on  this  warmup  data  and  using  the  heat  capacity  of  the  cold 
portion  of  the  resonator  a  value  for  the  total  heat  leak  to 
the  cold  end  of  the  refrigerator  is  determined  (in  mW/#c). 
This  number  is  then  reduced  by  2.4  mW/°c,  which  is  the 
thermal  conduction  of  the  stack  and  the  gas.  The  resulting 
quantity  is  the  conduction  external  to  the  stack  and  is 
accuarate  to  within  8%,  The  value  for  Qtotal  is  then  the 
sum  of  QHEAT  +  HEATLK.  The  error  in  QHEAT  is  negligible, 
but  the  error  in  HEATLK  is  not.  The  error  in  Qtotal  is 
about  8%  when  QHEAT  is  zero  and  is  about  1.5%  to  2%  for 
values  of  Qtotal  near  heat  loads  of  optimum  efficiency. 

The  remaining  quantity  necessary  for  the  calculation  of 
COP  is  the  acoustic  power  (or  work).  We  use  a  relation 
similar  to  that  derived  by  Hofler  (1986).  The  difference 
being  that  we  no  longer  use  the  velocity  signal  from  the 
bandpass  filter;  using  the  velocity  signal  from  the  lock-in 
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instead.  Also,  we  have  added  two  correction  terns  to  the 
equation.  We  are  using  an  equation  of  the  form: 

u 

1  Wfvv 

where  Hp  is  th^  dynamic  pressure  transducer  sensitivity, 
My  is  the  volume  velocity  sensitivity,  Vp  is  the  measured 
pressure  signal,  is  the  measured  velocity  signal,  (  1- 

i.)  is  an  accelerometer  non-linearity  correction  term,  and 
is  a  volume  velocity  calibration  correction  term. 
<A=(7.4xl0  ^  sec2/cm)a,  where  a.  is  the  acceleration  in 

_3 

cm/sec2  ( x  10  ),  as  shown  in  column  13  of  Appendix  C. 

«vv=0. 05( f-550) ,  where  f  is  the  frequency  in  hertz,  as 
shown  in  column  4  of  Appendix  C.  The  ratio  My/Mp  is  the 
result  of  the  small  cavity  volume  velocity  calibration, 
discussed  previously,  and  is  therefore  just  a  number  that 
we  put  into  the  equation  prior  to  calculation.  Since  the 
refrigerator  was  disassembled  on  several  occasions  during 
our  experiments  this  calibration  was  performed  several 
times.  The  calibration  value  changes  over  time  as  the 
driver  is  used  and  therefore  this  number  My/Mp  is  corrected 
in  our  equations  for  all  data  taken  after  each  calibration. 
The  quantity  Mp  is  also  a  calculated  number  that  is  input 
into  the  equations,  since  it  is  just  the  dynamic  pressure 
transducer  sensitivity.  The  actual  equations  that  we  have 
used  are  of  the  following  form: 

W=598 10 (  l-7.4xl0-5a)Vp(0.003VLI )/( 1839. 8-0. 05 ( f-550 ) ) , 
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where  59810  is  (l/Mp)4  including  unit  conversions,  (  1- 
7.4xl0~5a)  is  a  correction  term  for  accelerometer  non- 
liearities,  as  discussed  above,  (0.003)/^)  is  the  measured 
velocity  signal  (with  unit  conversions  in  the  coefficient), 
1839.8  is  (My/Hp)  the  result  of  our  volume  velocity 
calibration,  and  (0. 05(f-550) )  is  a  correction  to  the 
volume  velocity  term  due  to  the  calibration  being  frequency 
dependent,  as  discussed  above.  The  accelerometer  correction 
(l-(ft)  is  about  2%  at  most,  and  the  frequency  correction 
( *Vy)  is  about  IX  at  most.  The  frequency  correction  was 
obtained  for  the  case  of  pure  helium.  However,  given  the 
much  larger  frequency  range  of  the  gas  mixtures,  the  amount 
of  correction  in  some  cases  is  too  small,  giving  a 
calibration  error  of  about  1%.  The  total  worst  case  error 
for  W  is  about  5%  with  3.5%  being  a  more  typical  value. 

Using  the  Qtotal  and  H  values  calculated  by  the  above 

methods  w«  get  COP-Qtotal/W.  To  get  COPR  we  need  to  divide 

this  by  the  Carnot  COP  for  the  same  temperature  difference. 

The  Carnot  COP  is  known  to  be  COP_ _  . =Tc/(Th-Tc ) ,  which 

carnot 

means  COPR  can  be  derived  as  COPR=COP( (Th/Tc)-1) .  This  is 
the  equation  we  have  used  to  determine  COPR.  On  our 
spreadsheet,  in  Appendix  C  this  is  the  C0PR2  column.  The 
worst  case  error  in  the  absolute  values  for  COPR  may  be  as 
high  as  15%  however,  typical  values  are  about  5%.  Also,  the 
relative  error  involved  in  comparing  data  taken  under 
different  conditions  is  probably  considerably  less  than  5%. 
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We  now  have  all  of  the  values  necessary  to  evaluate  the 
effects  of  our  modifications  on  the  refrigerator’s 
coefficient  of  performance. 

B.  RESULTS  OF  EXPERIMENT 

As  stated  earlier,  we  will  break  this  discussion  into 
three  major  sections.  First  we  will  discuss  the  results  of 
our  stack  modifications.  This  will  be  followed  by  a 
discussion  of  performance  with  different  gases  and  the  same 
stack.  We  will  conclude  with  the  overall  improvements  using 
our  best  stack/gas  combination.  Also,  a  brief  discussion 
concerning  the  results  with  different  pressure  amplitudes 
and  different  resonator  lengths  will  be  provided. 

Since  we  have  data  from  two  new  stacks  and  two  gas 
mixtures,  as  well  as  two  different  pressure  amplitudes  and 
two  resonator  lengths,  available  for  comparisons,  we  muj;t 
be  selective  in  our  presentation,  so  as  not  to  have  an 
excessive  number  of  plots  that  show  the  same  result.  To 
fulfill  this  goal  we  will  present  one  or  two  sample  sets  of 
curves,  with  analysis,  for  each  major  category  of  interest, 
and  provide  numerical  results  for  the  other  cases  within 
that  category,  where  apprpriate.  In  general  a  set  of  curves 
includes  a  plot  of  COPR  versus  Qtotal  and  a  plot  of  Tc/Th 
versus  ototal. 

Prior  to  presenting  any  of  our  figures  we  shall  intro¬ 
duce  the  notation  that  is  used  in  the  titles  and  legends. 
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All  of  our  plots  are  designed  to  provide  data  comparisons 
and  as  such  each  contains  nore  than  one  curve.  The  legend 
lists  the  specific  curves  using  the  following  notation: 

GAS  TYPE — STACK — RESONATOR  LENGTH— PRESSURE  AMPLITUDE, 
where  abbreviations  have  been  used  in  all  cases.  For  gas 
types  we  have  helium,  helium-argon,  and  helium-xenon  which 
are  represented  by  He,  HeAr,  and  HeXe,  respectively.  For 
the  stacks  we  have  old,  new,  and  dual  which  are  represented 
by  0,  N,  and  D,  respectively.  The  resonator  length  is 
designated  by  s  or  1  for  short  or  long  and  the  pressure 
amplitude  is  represented  by  the  percent,  relative  to  the 
mean  pressure  (Po/Pm),  by  the  number  2  or  3. 

1 .  Stack  modification  results 

To  illustrate  the  effects  of  our  stack  modifications 
we  have  elected  to  compare  the  data  for  pure  helium,  at  2% 
pressure  amplitude,  with  the  short  resonator  installed. 
Figure  V-2  shows  the  temperature  ratio  (Tc/Th)  versus  total 
heat  load  (Qtotal)  curves  for  all  three  stacks  under  these 
conditions.  It  is  fairly  clear  from  the  curves  that  the 
data  are  essentially  linear  for  all  three  cases,  as 
expected.  The  curves  also  show  that  for  a  given  temperature 
ratio  the  old  stack  can  support  more  heat  load  than  the 
dual  stack  and  the  new  stack  can  support  more  heat  load 
than  either  of  the  others.  This  statement  can  be  reversed, 
such  that  for  a  given  heat  load  the  temperature  spanned  is 
greatest  for  the  new  stack  and  lowest  of  the  dual  stack. 
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Q  total  (watt*) 


Temperature  Ratio  vs.  Qtotal,  For  Helium,  Short  Resonator 


Those  curves  also  show  that  our  data  are  following  the 
expected  trends  for  the  refrigerator,  in  general.  Figure  V- 
3  shows  COPS  versus  Qtotal  for  all  three  stacks.  It  is 
clear  from  this  Figure  that  the  COP  is  significantly 
improved  with  the  new  stack  and  vastly  degraded  with  the 
dual  stack.  Evaluating  the  difference  at  the  peaks  of  the 
curves,  we  see  that  the  dual  stack  has  a  peak  COPS  of  0.034 
while  the  old  stack  has  a  peak  at  0.06?  and  the  new  stack 
peaks  at  0.078.  This  shows  a  16X  improvement  for  the  new 
stack  and  a  51X  reduction  for  the  dual  stack.  Conducting 
this  same  analysis  for  the  remainder  of  our  data  sets 
yields  an  average  improvement  of  14X  with  the  new  stack  and 
an  average  reduction  of  57X  with  the  dual  stack. 

It  is  clear  from  these  results  that  the  dual  stack 
has  not  performed  as  predicted  while  the  new  stack  has 
provided  significant  improvements.  It  is  believed  that  the 
basic  theory  used  to  design  the  dual  stack  is  valid,  in-so- 
far  as  the  concept  is  concerned.  However,  it  is  obvious 
that  we  have  underdesigned  the  stack  by  either  omitting  an 
important  factor  or  over  simplifying.  At  the  time  of  design 
and  construction  of  the  stack,  it  was  our  belief  that  plate 
spacing,  so  long  as  it  was  not  too  small,  would  not  have 
any  detrimental  effects.  Based  on  this  belief,  and  an 
initial  evaluation  of  the  new  stack’s  performance,  we 
decided  to  build  the  dual  stack  with  the  hot  end  having  the 


same  fishing  line  size  as  that  used  in  the  new  stack. 
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Figure  V-3,  COPR  vs.  Qtoial  for  Hollum,  Short  Resonator 


However, 

witn 

the 

dual  layer  at 

the 

hot 

end  and  the  thick 

fishing 

\  me 

at 

the  cold  end, 

we 

ended  up  with  a  20% 

increase 

in 

plate 

spacing  at 

the 

hot 

end  and  a  163% 

increase 

in 

plate 

spacing  at  the 

cold 

end 

(relative  to  the 

spacing  m  the  new  stack).  Based  on  -the  extreaely  poor 
results  with  the  dual  stack  it  appears  that  the  magnitude 
of  the  plate  spacing  is  aore  critical  than  we  had 
originally  thought.  The  theory,  presented  in  Chapter  II, 
shows  that  a  plate  spacing  of  2  to  4  thermal  penetration 
depths,  with  the  entire  resonator  cross-sect  ion  filled  with 
plates  is  desired  for  the  best  heat  puaping 
characteristics.  However,  our  analysis  of  the  theory  did 


not  predict 

this 

massive 

drop 

in  COP 

with  a  larger  plate 

spacing.  To 

some 

extent 

this 

plate 

spacing  effect  is 

i 1  lustrated 

by 

the  new 

stack 

*s  performance,  as  well. 

Although  our  intuitive  reasoning  for  the  design  of  the  new 
stack  predicted  improved  performance,  the  results  are 
better  than  expected.  Combining  this  with  the  extremely 
poor  performance  of  tho  dual  stack,  and  the  fact  that  the 
new  stack  has  the  smallest  plate  spacing  of  the  three 
stacks,  ue  can  see  that  plate  spacing  is  much  more 
important  than  previously  thought.  Rough  calculations  for 
helium  show  that  the  average  plate  spacing,  in  number  of 
thermal  penetration  depths  is  3.3  for  the  new  stack,  3.5 
for  the  old  stack,  4.0  for  the  hot  end  of  the  dual  stack, 
and  8.3  for  the  cold  end  of  the  dual  stack. 
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2.  Results  for  Gas  Mixtures 


For  this  analysis  we  have  elected  to  present  two 
sets  of  data,  both  of  which  »re  for  the  new  ttack,  This  is 
primarily  due  to  the  fact  that  the  new  stack  is  the  only 
one  for  which  we  have  any  helium-argon  data.  Also,  based  on 
the  stack  analysis,  the  new  stack  is  the  best  one 
available.  For  this  analysis  we  will  discuss  the  results 
for  2x  pressure  amplitude  and  both  the  short  and  long 
resonator.  The  long  resonator  consists  of  a  one  inch  spacer 
that  lengthens  the  resonator  and  thereby  changes  the 
position  of  the  stack  realtive  to  the  acoustic  standing 
wave . 

Figures  V-4  and  V-5  show  the  plots  of  Tc/Th  versus 
Qtotal  for  the  short  and  long  resonators,  respectively.  It 
is  clear  from  these  curves  that  all  of  the  data  are  linear, 
as  expected.  It  is  interesting  to  note  that  in  Figure  V-2 
all  of  the  curves  for  pure  helium  were  essentially 
parallel,  even  though  the  data  was  for  different  stacks. 
Now,  in  Figures  V-4  and  V-5  we  see  that  the  curves  are  no 
longer  parallel  for  the  different  gases,  even  though  the 
data  are  all  for  the  same  stack.  This  is  believed  to  be 
related  to  the  fact  that  the  heat  pumping  capacity  is 
related  to  the  PU  product  (pressure  times  volume  velocity). 
For  helium  gas  at  constant  pressure  amplitude  the  PU 
product  will  vary  in  a  consistent  manner,  with  respect  to 
temperature,  regardless  of  stack  design.  For  different 


Q  total  (watts) 

Figure  V-4.  Tampareture  Ratio  vs.  Qtotai,  Now  Stack,  Short  Resonator,  All  Qasaa 

\\z 
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Figure  V-5.  Tamprtature  Ratio  vs.  Qtotal,  New  Stack,  Long  Resonator,  All  Gases 
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gases  however,  the  veloc^ies  Mill  vary  because  of 
different  sound  speeds,  end  therefore  even  et  constant 
pressure  aaplitude  the  PU  product  Mill  vary  differently, 
producing  curves  Mith  different  slopes. 

The  plots  of  COPR  versus  Qtotal  for  the  short  and 
long  resonator  are  shoun  in  Figures  V-6  and  V-7, 
respectively.  It  is  clear  froa  all  of  thes?  cur\es  that 
hel iua-argon  produces  a  certain  degree  of  iaproveaent  in 
COP,  over  pure  heliua,  and  that  hel iua-xenon  produces  the 
highest  COPS.  For  the  short  resonator  curves  of  Figure  V-6 
we  see  that  the  peak  COFR  values  are  approximately  $.078 
for  heliua,  $.102  for  hel iua-argon,  and  0.129  for  heliua- 
xenon.  Froa  these  values  Me  see  that  hel iua-argon  gives  a 
3 lx  increase  over  heliua,  Mhile  hel iua-xenon  provides  a  65X 
iaproveaent  over  heliua  and  a  27X  iaproveaent  over  hel iua- 
argon.  It  should  be  noted  that  the  peak  COPR  values  are 
shifted  to  louer  heat  load  values  for  the  gas  mixtures. 
This  is  believed  to  be  due  to  a  reduction  in  this  PU 
product  that  uas  discussed  above.  Similar  calculations  for 
other  stacks  and  other  pressure  amplitudes  show  an  average 
improvement  in  the  peak  COPR  values  of  61X,  when  comparing 
helium-xenon  to  helium.  There  is  no  other  short  resonator 
data  available  for  helium-argon.  The  long  resonator  curves 
of  Figure  V-7  have  been  included  primarily  because  they 
illustrate  the  peak  COPR  value  that  has  been  measured  to 
date  of  0.197  for  helium-xenon.  Although  the  helium-argon 
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Figure  V-6.  COPR  vs.  Qtotal,Naw  Stack,  Short  Resonator,  All  Gases 
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Figure  V-7.  COPR  vs.  Qtotal,  New  Stack,  Long  Resonator,  All  Gases 


curve  appears  somewhat  incomplete  we  feel  safe  in  assuming, 
based  on  all  of  our  other  curves  and  data,  that  the  curve 
is  within  a  few  percentage  points  of  its  peak  value  at  the 
last  plotted  data  point.  We  will  therefore  use  that  datum 
for  our  peak  value.  These  curves  show  the  same  general 
trends  as  those  of  Figure  V-6,  however  the  magnitudes  of 
the  improvements  are  somewhat  reduced.  For  peak  values  of 
0.132  for  helium,  0.1595  (from  data  of  Appendix  C)  for 
helium-argon  and  0.197  for  helium-xenon,  we  see  that  helium- 
argon  only  produces  a  21%  improvement  over  helium  and 
helium-xenon  produces  a  49%  improvement  over  helium  with  a 
24%  improvement  over  helium-argon,  tt  is  believed  that  the 
improvement  percentages  with  the  long  rresonator  are 
reduced  because  the  long  resonator  itself  produces  about  a 
62%  improvement  over  the  short  resonator,  thereby  placing 
the  refrigerator  closer  to  a  maximum  possible  or  limiting 
COPR  value.  Thus  the  absolute  improvements  we  achieved  have 
a  smaller  effect,  percentage  wise. 

3.  Overall  Improvement 

In  this  section  we  present  the  overall  result  of  the 
thesis.  That  is,  the  overall  improvement  in  the  COP  of  the 
prototype  refrigerator,  taking  into  account  both  the  stack 
modifications  and  the  binary  gas  mixtures.  To  illustrate 
this  case  we  have  elected  to  compare  the  results  with  the 
old  stack  and  pure  helium  at  2%  pressure  amplitude  to  those 
achieved  with  the  new  stack  and  helium-xenon  at  2%  pressure 
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amplitude.  These  sets  of  data  were  selected  because  the 
original  refrigerator  used  only  the  old  stack  and  pure 
helium,  while  our  best  results  were  achieved  with  the  new 
stack  and  helium-xenon.  We  have  not  used  the  long  resonator 
data,  even  though  that  gives  our  highest  COPR,  because  we 
did  not  take  any  long  resonator  data  with  the  old  stack. 
However,  for  a  number  for  the  old  stack,  long  resonator, 
pure  helium  setup  we  can  use  the  0. 126  peak  COPR  given  by 
Hof ler( 1986)  in  his  doctoral  dissertation.  Comparing  this 
to  our  best  case  peak  COPR  of  0. 197  we  get  a  56% 
improvement. 

For  the  comparison  of  our  experimental  data  we  will 
present  the  COPR  versus  Qtotal  plot  only,  since  the  Tc/Th 
versus  Qtotal  plots  have  already  been  discussed  for  this 
data.  Figure  V-8  shows  our  plot  for  overall  improvement. 
With  peak  COPR  values  of  0.129  for  the  helium-xenon  and 
0.067  for  the  helium  we  see  that  the  overall  improvement  in 
the  refrigerator’s  coefficient  of  performance,  due  to  our 
modifications,  is  93%. 

4.  Other  Effects  of  Interest 

There  are  two  other  areas  of  interest  that  have  been 
mentioned  briefly  during  our  discussions  of  results.  These 
are  the  effects  produced  by  changing  resonator  length  and 
those  produced  by  changing  the  pressure  amplitude. 


In  Figure  V-9  we  present  the  data  for  all  three 
gases  in  the  new  stack  with  both  short  and  long  resonators. 


COPR 


Q  total  (watts) 

Figure  V-9.  COPR  vs.  Qtotal,  Showing  Resonator  Effects 
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From  this  we  can  see  that  just  the  shift  from  short  to  long 
resonator  produces  increases  in  COPS  as  follows:  70X  for 
helium,  59%  for  helium-argon  and  52%  for  helium-xenon.  This 
rather  large  increase  in  COP  is  due  primarily  to  the  shift 
in  the  position  of  the  stack  relative  to  the  standing  wave. 
With  the  long  resonator  installed,  the  hot  end  of  the  stack 
is  moved  closer  to  the  velocity  antinode  and  the  heat 
transfer  power  is  increased  significantly,  thereby 
producing  the  large  increase  in  the  COP  and  the  observed 
shift  of  the  COPS  peak  towards  higher  heat  loads. 

In  Figures  V-10  and  V— 1 1  we  present  a  comparison 
between  both  gas  types  and  pressure  amplitudes.  In  Figure  V- 
10  we  have  helium-xenon  with  the  short  resonator,  both 
pressure  amplitudes  and  all  three  stacks.  In  Figure  V-ll  we 
have  all  of  the  same  curves  for  pure  helium.  In  all  of 
these  curves  it  is  obvious  that  at  higher  pressure 
amplitudes  the  peak  COPR  value  has  shifted  toward  higher 
heat  loads.  Based  on  this  and  the  discussion  of  Figures  V-4 
and  V-5  it  is  clearly  shown  that  heat  transfer  power  (or 
heat  pumping  capacity)  increases  with  pressure  amplitude 
(due  to  the  PU  product  increasing).  It  is  interesting  to 
note  that  the  peak  COPR  value  is  in  fact  lower,  at  higher 
pressure  amplitudes,  for  all  of  the  cases  for  which  we  have 


data, 

with  the  exception 

of 

the 

pure 

helium,  short 

resonator,  old 

stack 

case. 
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Figure  V-10.  COPR  vs.  Qtotal  for  HeXe,  Short  Resonator,  All  Stacks  and  Amplitudes 
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Figure  V-11.  COPR  vs.  Qtotal  for  He,  Short  Resonator,  All  Stacks  and  Amplitudes 
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pressure  amplitude  for  all  of  his  data.  In  fact  our  data 
with  his  stack  and  helium  also  follow  this  trend. 
Additionally,  Hofler  (1986)  reported:  "there  is  also  some 
indication  that  the  magnitude  of  these  maxima  approach  some 
constant  value  at  higher  amplitudes."  Our  results  for  gas 
mixtures  and  other  stack  designs  tend  to  contradict  this, 
and  lead  toward  a  rather  alarming  conclusion  that  there  is 
some  relatively  low  pressure  amplitude  at  which  the  COPR 
maxima  peak  and  above  which  the  COPR  maxima  drop  off  again. 
This  seems  to  indicate  that  there  may  be  some  non-* 
linearity,  such  as  turbulence,  that  limits  the  best  OOP’s 
to  relatively  low  amplitudes. 
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VI.  CONCLUSIONS  AMD  RECOMMENDATIONS 


The  primary  purpose  of  this  thesis  is  to  provide 
modifications  for  the  Space  Thermoacoustic  Refrigerator 
(STAR)  that  will  improve  its  thermodynamic  efficiency, 
specifically  its  coefficient  of  performance  (COP).  Through 
the  use  of  a  helium-xenon  binary  gas  mixture  and  a  new 
stack  design,  with  less  string  and  reduced  plate  spacing, 
we  have  produced  a  93X  increase  in  the  peak  value  of  COP. 

Along  with  this  improvement  we  have  raised  several 
questions  that  are  as  yet  unanswered.  Based  on  our  results 
with  a  long  resonator,  relative  to  those  with  a  short 
resonator,  it  would  appear  that  there  is  some  optimum 
resonator  length  and  design  that  must  be  found.  Further 
testing  with  different  resonator  spacer  thicknesses,  must 
be  performed  to  resolve  this  issue.  The  unexplained  drop 
m  peak  COP  relative  to  Carnot  (COPR),  with  increasing 
pressure  amplitude  opens  the  door  for  more  testing  with 
different  setups  and  many  pressure  amplitudes.  It  is  clear 
that  the  binary  gas  mixture  of  helium-xenon  is  superior 
under  all  of  the  conditions  that  we  have  observed.  However, 
the  failure  of  the  dual  stack  and  the  better  than  expected 
results  with  the  new  stack  leave  the  final  stack  design  as 
an  open  item.  It  is  clear  that  another  dual  stack  and  a 
different  version  of  the  new  stack  should  both  be  built. 
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with  saallar  pUtt  spacings,  to  deteraine  the  optima  stack 
design.  It  is  the  answers  to  these  question  that  will  yield 
the  design  capable  of  producing  the  aaxiaua  theraodynaaic 
efficiency  for  the  Space  Thoraoacousttc  Refrigerator 
(STAR). 

Although  this  project  deals  specifically  with  the  space 
application  of  our  refrigerator ,  there  is  no  reason  that  it 
should  be  Halted  to  that  application.  Due  to  the  lack  of 
aoving  parts  that  are  present  m  aost  other  cooling  engines 
(i.e.  sliding  seals)  this  unit  should  be  extreaely 
reliable.  It  also  has  reasonable  efficiency  relative  to 
aost  conaercial  refrigeration  devices.  The  fact  that  it 
does  not  use  chlorof lourocarbons  is  a  great  advantage, 
especially  with  the  close  scrutiny  that  this  ozone 
destroying  substance  is  currently  receiving.  Therefore, 
this  type  of  refrigeration  device  could,  and  should  be 
developed  for  coaaercial  use  and  should  be  considered  as  a 
possible  candidate  for  virtually  any  type  of  cooling  need. 
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APPENDIX  A.  PRANDTL  NUMBER  PROGRAM 

This  appendix  provides  a  program  and  documentation  for 
the  calculation  of  Prandtl  numbers  for  binary  gas  mixtures. 
The  program,  written  by  the  author,  is  designed  to 
calculate  Prandtl  numbers  for  all  concentrations  of  a 
binary  mixture  of  Helium  and  one  other  gas.  The  program 
is  written  in  Fortran  for  use  on  the  Naval  Postgraduate 
{  school  mainframe  computer,  with  a  Hatfor-77  compiler.  Two 

I 

output  files  are  created  when  the  program  is  run.  The 
first,  called  "PRANTL  OUTPUT  Al",  is  a  listing  of  Prandtl 
numbers  versus  helium  fraction  for  further  use  in  a 
plotting  program,  to  create  Prandtl  number  versus  helium 
fraction  curves  for  the  gas  mixture  selected.  The  second, 

i  called  "PRANTL  RESULTS  A l"  is  a  tabular  listing  of  Prandtl 

I 

number,  thermal  conductivity  (in  cal/cm-sec-K) ,  and  dynamic 
viscosity  (in  gm/cm-sec)  versus  helium  fraction  for  the 

» 

selected  gas  mixture.  The  listings  go  from  pure  helium  to 
zero  percent  helium  in  one  percent  intervals  to  provide 

j  easy  access  to  any  specific  mixture.  All  of  the  theory  and 

f 

equations  necessary  to  create  the  program,  other  than 
knowledge  of  a  programming  language,  can  be  found  in 
[  "MOLECULAR  THEORY  of  GASES  and  LIQUIDS"  by  Hirschfelder , 

l  Curtiss,  and  Bird,  John  Wiley  &  Sons,  Inc,,  New  York,  1954. 

! 

\  The  primary  source  was  chapter  eight  of  the  text,  hereafter 
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referred  to  as  Hirschfelder,  although  section  and  aquation 
numbers  will  be  specified  where  necessary  for 
clarification. 

As  written,  this  program  was  designed  for  personal  use, 
exclusively  to  derive  the  data  necessary  to  create  the 
previously  mentioned  curves  of  Prandtl  number  versus  helium 
fraction.  It  has  been  modified  to  include  the  second 
tabular  output  for  completeness  and  further  use  by  others 
in  this  field.  The  program  is  currently  written  to  require 
hard  input  and  therefore  requires  editing  of  the  mam 
program  for  each  change  of  gas  type  or  temperature.  In 
addition  to  satisfying  the  specified  needs,  this  method  was 
used  due  to  there  being  thirteen  input  parameters  that  vary 
with  gas  type  and  temperature,  five  of  which  require  table 
lookups  in  Hirschfelder.  This  program  could  be  modified  for 
incorporation  into  a  larger  program  as  a  subroutine, 
provided  the  specific  gas  types  and  temperatures  of 
interest  are  completely  specified  and  all  of  the  input 
parameters  are  looked  up  and  incorporated  into  a  data  array 
for  reading.  This  has  not  been  done  to  date  as  the  current 
useage  did  not  require  these  features. 

As  previously  stated  the  primary  source  for  this  program 
is  chapter  eight  of  Hirschfelder,  which  provides  the 
equations  necessary  to  calculate  various  transport 
coefficients,  such  as  the  coefficient  of  viscosity  and  the 
coefficient  of  thermal  conductivity.  The  equations 
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presented  in  chapter  eight  rely  heavily  on  the  elementary 
theory  of  transport  phenomena  presented  in  section  1.2  of 
the  text  as  well  as  the  rigorous  kinetic  theory  for 
monatomic  gases  presented  in  chapter  seven.  The  intricate 
details  of  this  theory  will  not  be  presented  here,  however 
a  discuss ior  of  the  methods  used,  including  the  main 
assumptions  and  restrictions,  will  bte  provided  as  well  as  a 
definition  of  each  of  the  terms  or  symbols  in  the  program,  >. 

The  development  of  the  kinetic  theory  for  gases  is  based 
on  knowledge  of  the  distribution  function  that-  represents 
the  number  of  molecules  of  a  specific  species,  in  a  unit 
volume  element  about  a  point  in  space,  with  velocities  in  a 
unit  range  about  a  specified  velocity,  at  a  given  instant 
in  time.  Since  the  usual  definitions  for  transport 
coefficients  only  apply  under  conditions  of  equilibrium  (or 
only  slightly  different  from  equilibrium)  this  restriction 
applies  to  the  results.  Under  this  limit  the  distribution 
function  is  nearly  Haxwellian  and  is  solvable  by  a 
perturbation  method  developed  by  Chapman  and  Enskog.  The 
solutions  are  then  used  to  obtain  expressions  for  transport 
coefficients  in  terms  of  a  set  of  integrals,  Omega(l,s), 
which  involve  explicitly  the  dynamics  of  molecular 
encounters.  The  Chapman-Enskog  theory  relies  on  several 
assumptions  and  therefore  has  limited  applicability.  Each 
assumption  will  be  briefly  discussed  and  related  to  the 
specific  situation  of  our  experiment. 
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Only  binary  collisions  are  considered,  which  limits  us 


to  densities 

low 

enough  that  three 

body 

collisions  are 

unimportant. 

The 

theory 

predicts 

that 

viscosity  is 

independent 

of 

pressure 

at  constant 

temperature. 

Experiments 

on  Nitrogen  have 

shown  that 

only 

a  4X  deviation 

in  viscosity 

occurs  as  pressure  is 

varied  form  1  to  60 

atmospheres . 

Our 

operating 

condition 

of  10 

atmospheres  is 

well  within  this  and  should  therefore  closely  follow  the 
theory  in  this  respect. 


The  theory  further  limits  itself  to  Classical  Mechanics 
and  thereby  precludes  low  temperature  situations  where 
quantum  effects  become  significant.  These  quantum  effects 
are  less  than  lx  for  helium  above  200  Kelvin  (and  even 
smaller  for  heavier  isotopes).  The  coldest  temperature  we 
have  achieved  to  date  has  been  193  Kelvin  under  a  no  load 
condition  so  our  normal  operations  are  all  above,  200  Kelvin 
and  we  are  not  concerned  with  quantum  effects. 

The  Chapman-Enskog  theory  uses  a  series  approximation, 
to  the  above  mentioned  distribution  function,  in  solving 
the  Boltzmann  equation  and  then  uses  only  the  first 
approximation,  which  is  only  valid  when  the  gradients  of 
the  physical  quantities  (i.e.  density,  velocity, 
temperature)  are  small.  The  changes  of  these  properties 
over  one  mean  free  path  must  be  small  with  respect  to 
unity.  Based  on  the  mean  free  paths  for  gases  at  pressures 
above  one  atmosphere  this  condition  is  satisfied  except 
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under  conditions  of  extreme  gradients  such  as  those  present 
in  shock  waves.  Since  we  do  not  have  any  shock  waves  we  are 
not  restricted  by  this. 

The  theory  further  assumes  that  the  dimensions  of  the 
containing  vessel  and  any  obstacles  within  it  are  large 
compared  to  the  mean  free  path,  which  provides  no 
limitations  for  our  gases  and  our  pressures.  Although  the 
Chapman-Enskog  theory  applies  strictly  to  monatomic  gases 
the  results  have  been  shown  to  be  good,  even  for  polyatomic 
gases,  provided  that  the  molecules  are  not  to  non- 
spherical.  Finally,  the  development  of  the  theory  results 


in 

the 

set 

of  integrals  Omega(l,s), 

which  are  present  in 

Ki\) 

all 

of 

the 

equations  for  transport 

coefficients.  These 

integrals,  in  turn  are  dependent  upon  the  potential 
function  that  represents  the  molecular  interactions. 
Throughout  the  derivation  in  Hirschfe lder ,  and  for  all  of 
the  equations  used  here,  the  Lennard-Jones  (6-12)  potential 
is  assumed. 

Since  some  of  the  other  terms  and  symbols  depend  upon 


the  Omega 

terms  in 

the 

program,  we 

Will 

start 

with  a 

discussion 

of  where 

the 

Omega  terms 

come 

f  rom . 

These 

Omega ( 1 , s ) 

terms  represent  a  set  of 

integrals, 

,  linear 

combinations 

of  which 

are 

used  to  evaluate  the 

more 

complex 

bracket  integrals  arrived  at  in  solving  for  the  expansion 
coefficients  of  the  Sonine  polynomials.  The  transport 
coefficients  are  expressed  in  terms  of  the  Sonine 
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polynomial  expansion  coefficients  and  can  then  be  expressed 
in  teras  of  these  Omega(l.s)  integrals.  In  Hirschfelder 
there  is  also  an  Onega-star ( 1 , s )  symbol  which  represents 
the  deviation  of  a  particualr  molecular  model  (in  our  case 
the  Lennard-Jones  (6-12)  potential)  from  the  idealized 
rigid  sphere  model  of  molecular  interactions.  All  of  the 
Omegas  in  the  program  represent  an  Omega-star  value.  The 
Omega  integrals  are  functions  of  the  reduced  temperature 
(T*=kT/()  and  are  tabulated  in  table  I-M  of  Hirschfelder. 

Where  more  than  one  variable  of  a  given  type  exists, 
such  as  sigmal,  sigma2,  and  sigm!2,  the  1  and  2  represent 
the  value  of  the  quantity  for  pure  gas  1  or  pure  gas  2 
(which  is  always  helium  for  our  purposes),  respectively. 
The  12  represents  the  value  for  gases  1  and  2  together. 
This  is  true  for  all  cases  except  ETA  12  and  LAH12.  These 
two  variables  represent  artificial  quantities  for  a 
hypothetical  pure  substance  with  molecules  of  molecular 
weight  2 ( M 1 ) ( H2 ) / ( M 1+M2)  that  interact  according  to  a 
potential  curve  specified  by  interaction  parameters  sigmal2 
and  (12.  These  parameters  are  defined  in  terms  of  the 
Lennard-Jones  (6-12)  potential  in  chapter  one  of 
Hirschfelder,  where  sigma  is  the  distance  of  closest 
approach  (in  angstroms)  of  two  molecules  which  collide  with 
zero  initial  relative  kinetic  energy  and  (  is  the  maximum 
energy  of  attraction  of  the  two  molecules.  Table  I-A  of 
Hirschfelder  provides  values  for  sigma  and  (/k  for  various 
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substances,  where  k  is  the  Boltzmann  constant  and  (/k  has 
units  of  Kelvin.  The  parameters  (12  and  sigml2  are  arrived 
at  through  empirical  combining  laws  which  relate  force 
constants  between  unlike  molecules  to  those  between  like 
molecules.  These  laws  are  sigml2=( 1/2) (sigmal+sigma2)  and 

u 

tl2=((l*(2)  .  The  quantities  A12  and  B12  are  merely 
frequently  encounterd  ratios  of  different  Omegas  that  have 
been  tabulated  separately  as  functions  of  T*  in  table  I-N 
of  Hirschfelder . 

The  remaining  terms  or  symbols  in  the  program  are  fairly 
straightforward,  based  on  the  above  explanations  and  each 
type  will  now  be  defined: 

PR=prandtl  number 
HEshelium  fraction 

LAM=thermal  conductivity  coefficient  (in  cal/cm-sec-K) 
ETA=coeff icient  of  dynamic  viscosity  (in  gm/cm-sec) 
M=molecular  weight 

CP=specific  heat  capacity  at  constant  pressure 
(in  cal/gm-K) 

T=absolute  temperature  (in  Kelvin) 

Xl=fraction  of  non-helium  gas 

X2=helium  fraction 

I=counting  variable  for  the  Do  loop 

The  output  parameters  are  Prandtl  number,  helium 
fraction,  thermal  conductivity,  and  dynamic  viscosity.  They 
are  represented  by  PR ( I ) , HE ( I ) , LAM I X ( I ) ,  and  ETAMIX(I) 
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respectively.  The  remaining  parameters  are  all  merely 
mathematically  defined  intermediate  steps  in  the 
calculations  of  the  output  parameters  and  have  no  physical 
meaning.  The  terms  used  for  viscosity  are  XETA,  YETA,  and 
ZETA,  and  the  terms  used  for  thermal  conductivity  are  Ul, 
U2,  UY,  UZ,  XLAM,  YLAM,  and  ZLAM. 

This  section  will  now  provide  the  most  rececnt  version 
of  the  program.  Since  the  explanation  of  the  input  is 
fairly  lengthy  it  will  be  covered  first,  rather  than 
disrupt  the  flow  of  the  program.  This  will  include  a  list 
of  all  required  input  parameters,  their  sources  where 
appropriate,  and  a  sample  set  of  inputs  that  were  used. 
There  are  thirteen  input  quantities,  A12,  B12,  0MEG12,  HI, 
M2,  CPI,  CP2,  SIGMA1,  SIGHA2,  SIGM12,  0HEGA1,  0MEGA2,  AND 
T.  After  selecting  Ml,  M2,  and  T  you  have  specified  your 
gas  combination  and  temperature  of  interest.  Values  for  CP 
may  be  obtained  from  various  sources,  such  as  the  CSC 
tables,  thermodynamics  textbooks,  or  standard  tables  of 
thermodynamic  properties.  The  values  for  SIGMA  as  well  as 
i/ k  are  obtained  from  table  I-A  of  Hirschf e lder .  Using  t/k 
for  each  gas  and  (12=((1*(2)  you  calculate  Tl*,  T2*,  and 
T12*  (where  T*=kT/0.  Now  using  the  appropriate  T*  you  go 
to  table  I-M  of  Hirschfelder  for  the  OMEGA  values  and  table 
I-N  for  A12  and  B12.  This  provides  a  complete  set  of  input 
quantities.  The  following  table  is  a  sample  set  of  inputs 
for  two  gas  combinations  at  two  different  temperatures: 
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GAS  COMBINATION  AND  TEMPERATURE 


QUANTITY 

He-Xe 

He-Xe 

He-Ar 

He-Ar 

OF  INTEREST 

20® 

300 

200 

300 

Ml 

131.3 

131.3 

39.948 

39.948 

M2 

4.0026 

4.0026 

4.0026 

4.0026 

SIGHA1 

4.055 

4.055 

3.418 

3.418 

SIGHA2 

2.556 

2.556 

2.556 

2.556 

( i/K)  1 

229 

229 

124 

124 

( */K)2 

10.22 

10.22 

10.22 

10.22 

OMEGA  1 

1.700 

1.394 

1.275 

1.  104 

0MEGA2 

0.7467 

0.7485 

0.7467 

0.7485 

CPI 

0.0382 

0.0382 

0. 12477 

0. 12477 

CP2 

1.2512 

1.2512 

1.2512 

1.2512 

( €/K) 12 

48.38 

48.38 

35.60 

35.60 

T 12* 

4.  134 

6.201 

5.618 

8.427 

A 12 

1.098 

1. 1034 

1.  102 

1.  108 

B12 

1.093 

1 . 0904 

1.091 

1 . 0905 

0MEG12 

0.9632 

0.8916 

0.9080 

0.8470 

PROGRAM  PRANTL 
C 
C 


C  First  the  array  variables  are  dimensioned  to  allow 
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C  calculating  and  plotting  the  desired  quantities.  Then 
C  the  real  quantities  are  defined  since  some  begin  with 
C  letters  that  would  make  thee  integers  by  default. 

C 
C 

DIMENSION  PR( 15®) ,HE( 15® ), LAMIX ( 15®, ETAHIX ( 15®) 

REAL  M 1 , M2 , OMEGA  1 , 0MEGA2 , OMEG 12 , LAM  1 , LAM2 , LAM  12 , PR , HE 
REAL  LAMIX, ETAHIX 
C 
C 

C  Now  comes  the  input  section  where  the  counting  variable 
C  is  initialized  to  zero  and  all  required  data  is  input. 

C 

C 

A  12=1.  1®8® 

B12=l . 0905 
OMEG 12=0. 847® 

1=0 

Y  39.948 
M2=4 . 0026 
CP1=0. 12477 
CP2=  1.2512 
SIGMA1=3. 4 18 
STnMA2-'  '56 

SIGM 12=0. 5*<  SIGMA 1+SIGNA2 ) 

1 36 


OMEGA  1=1.  1®4 


0HEGA2-Q. 7485 
T=30® 

C 

C 

C  Now  the  output  files  are  created  and  opened  and  headings 
C  are  placed  where  appropriate. 

C 

C 

OPEN ( 67 , F I LE=  * PRANTL  OUTPUT  AIM 
OPEN ( 68, FILE=* PRANTL  RESULTS  AIM 
WRITE(68, *) ’THIS  OUTPUT  FILE  IS  FOR  A  MIXTURE  OF 
*GASES  WITH  ’ 

WRITE (68, 1)H1,M2 

1  FORMAT ( ’ATOMIC  WEIGHTS  OF  »,F1®.7,*  AND* , FI®. 7, * , AND 

*  A  ’) 

WRITE (68, 2 )T 

2  FORMAT ( ’TEMPERATURE  OF  *,F11.7,’  DEGREES  KELVIN.*) 
WRITE(68,*) 

WRITE (68,*)*  PRANTL  HE  THERMAL 

*  DYNAMIC’ 

WR iTE ( 68, * ) *  NUMBER  FRACTION  CONDUCTIVITY 

*  VISCOSITY* 

WRITE(68, *) 

WRITE ( 68, * ) ’  CAL/CM-SEC-K 

*  GM/CM-SEC ’ 

WRITE(68,*) 
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c 

c 

C  Now  the  Do  loop  is  opened  to  perform  the  required 
C  calculations.  The  equations  are  from  chapter  eight  of 
C  Hirschfelder  and  the  equation  numbers  are  as  fellows: 

C  ETA1  and  ETA2  equation  8.2-18 
C  ETA  12  equation  8.2-21 
C  LAN  1  and  LAM2  equation  8.2-31 
C  LAM  12  equation  8.2-35 

C  U 1 , U2 , UY , UZ , XLAM , YLAM ,  and  ZLAH  equation  8.2-36 
C  XETA , YETA ,  and  ZETA  EQUATION  8,2-22 

C  The  equation  for  Prandtl  number  is  ETAMIXftCPMIX/LAMIX. 

C 

C 

DO  10  X 1=0. 0, 1.005,0,01 
1  =  1  +  1 

ETA  1=2 . 6693E-5*SQRT (M1*T)/(SI GMA 1*S I GHA 1 IOMEGA  1 ) 
ETA2=2 . 6693E-5*SQRT ( M2*T ) / ( S I GH A2*S I GHA2*0MEGA2 ) 

ETA  12=2 . 6693E-5*SQRT ( 2*M 1*M2*T/ ( H 1+M2 >  > / ( S I GH 12* 
*SIGM 12*0HEG 12 ) 

LAM  1= 1 . 989 1E-4*SQRT ( T/M 1 ) / <  S I GMA 1*S I GMA ItOMEGA 1 ) 
LAM2= 1 . 989 1E-4*SQRT ( T/M2 ) / ( SIGMA2*SIGMA2*0MEGA2 ) 

LAM  12= 1 . 989 1E-4*SQRT ( T* ( M 1+M2 ) / ( 2*H 1*M2 ) ) / ( ( SIGM 12 
***2 ) *OMEG 12 ) 

U 1= ( 4*A 12/ 15 ) - ( ( (B12/5)+( 1/ 12  )  )  *  ( M 1/M2 ) )+( ( (M1-M2) 
***2)/(2*Ml*M2) ) 
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U2»(4*Al2/15)-< ( (B12/5)+( 1/12) )*(M2/M1 ) )  +  ( < (M2-M1) 
*#*2)/(2*Ml#M2) ) 

UY»( (4*A12/15)*< ( (H1+M2)**2)/(4*N1*H2) )*(LAM12*LAH 
*  12/ ( LAM 1*LAM2 ) ) )-( (B12/S)+( 1/12) )-( ( ( 12*B12/32*A12 
*)-(25/(32*A12) ) )*< ( (H1-M2)**2)/(H1*H2) ) ) 

UZ=( (4*A12/15)*( ( (H1+H2)**9)/(4*M1#M2) )*( ( LAH12/ 

*LAH 1 )  +  ( LAH 12/LAH2 ) ) - 1 ) ) - ( B 12/5 ) - ( 1/12) 

X2= 1-X1 

XETA* C  X 1*X 1/ETA  1 ) + ( 2*X 1*X2/ETA 12 )  +  (  X2*X2/BTA2 ) 

YETA= ( 3*A 12/5 )* ( ( X 1*X 1*M 1/ ( H2*ETA 1 ) )  + ( X2*X2*M2/ (HI* 
*ETA2) )+(2*Xl*X2*( (M1+H2)**2)*(ETA12*ETA12)/(4*M1*N2* 
*ETA1*ETA2*ETA12) ) ) 

ZETA= ( 3*A 12/5 ) * ( ( X 1*X 1*M 1/H2 ) + ( X2*X2*M2/M 1 ) + ( ( 2*X 1* 
*X2)*( ( ( ( ( H 1+M2 ) **2 ) / ( 4*H 1*H2 ) )*( ( ETA 12/ETA 1 )+( ETA 12 
*/ETA2) ) )-l) ) ) 

XLAM= ( X  1*X 1/LAM  1 )  +  ( 2*X 1*X2/LAH 12 )  +  ( X2*X2/LAM2 ) 

YLAM= ( X 1*X 1*U 1/LAM 1 ) + ( 2*X 1*X2*UY/LAM 12 )  +  ( X2*X2*U2/ 
*LAM2) 

ZLAH=  ( X  1*X lSDl )  +  ( 2*X  1*X2*UZ )  +  ( X2*X2*U2 ) 

PR< I )  =  ( ( 1+ZETA ) * ( XLAM+YLAM ) * ( ( X 1*M 1*CP 1 )  +  ( X2*M2*CP2 ) 
*))/(( XETA+YETA ) * ( 1+ZLAM ) * ( ( X 1*M 1 ) + ( X2*M2 ) ) ) 

HE( I )=X2 

LAM I X ( I ) * ( 1+ZLAM ) / ( XLAM+YLAM ) 

ETAMIX( I  )  =  ( 1+ZETA) /(XETA+YETA) 

C 

C 
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C  Non  wo  mil  onto  to  tho  output  filo*  tho  dot*  that  h*s 
C  boon  colculotod. 

C 

C 

«RITB<e7,5)PRU),HEtl) 

UR I TE ( 68 , 6 ) PR ( I ) , HE ( I ) , LAM I X ( I ) , BTAH I X ( I ) 

5  FORMAT (F 10. 7, FI®. 2) 

6  FORMAT (Fl®. 7, F 12. 2, E 15. 7,115.7) 

1®  CONTINUE 

STOP 

END 
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APPENDIX  B.  PLOTTING  PROGRAM 


This  Appendix  provides  a  program,  written  By  the 
Author,  thAt  will  drew  curves  of  Prandtl  number  versus 
helium  f ruction  from  the  output  file,  "PRANTL  OUTPUT  Al", 
creeted  by  running  the  program  "PRANDTL  FORTRAN  Al"  thAt 
was  presented  in  Appendix  A.  This  program  lust  be  compiled 
in  FCRTVS  And  is  specif icully  designed  to  use  DIS3PLA  on 
the  NavaI  Postgraduate  School  mainframe  computer.  Various 
DISSPLA  subroutines  Are  celled  by  the  progrem  end  eech 
will  be  expleined  with  comment  statements,  DISSPLA 
requires  the  detA  thet  is  to  be  plotted  to  be  in  en  arrey 
format.  Each  of  the  arrays  used  here  are  set  to  a  maximum 
of  150  points.  If  modifications  are  made  to  "PRANDTL 
FORTRAN  Al"  such  that  more  than  150  points  are  calculated 
the  array  sizes  must  be  redefined.  Two  character  strings 
are  used  in  this  program.  They  are  NENIO  and  NEN20,  both 
of  which  are  used  to  rename  files.  The  first  changes  the 
input  filename  into  the  format  required  for  DISSPLA  and 
the  second  changes  the  file  back  to  its  original  name  for 
storage.  The  only  parameters  used  are  PR(I)  and  HE 1 1 ) 
which  represent  Prandtl  number  and  helium  fraction, 
respectively.  Both  of  these  quantities  are  real  numbers 
and  are  defined  as  such  in  the  program.  This  completes  the 
introduction  and  explains  the  first  three  lines  of  the 
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program.  The  remainder  of  this  appendix  Mill  bo  a  listing 
of  tho  program  with  explanatory  consents.  Latest  update  of 
prograa,  March  1988. 

Program  QtPPKN 
ch arac t er*8®  NEW  10 , NEW2* 
dimension  PK( 150) ,HE( 15®) 
real  PS, HE 

c 

c  The  next  three  lines  rename  the  input  file  to  the 
c  format  necessary  for  DISSPLA  and  open  it  for  use. 
c 

NEW  10* ‘rename  V/‘  PSANTL  OUTPUT  Al*//‘  file 
*XYZ1  A 1 ’ 
call  excast NSW  10) 
open <  67 , f i le= ‘ XYZ 1  * ) 
c 

c  In  the  next  seven  lines  a  Do  loop  is  used  to  read 
c  all  of  the  data  from  the  input  file  into  the 
c  specified  arrays.  Then  the  input  file  is  closed  and 
c  renamed  it  to  its  original  name.  The  format  line 
c  ensures  the  data  is  read  from  the  file  in  the  same 
c  format  as  it  was  written  into  the  file, 
c 

1  formattf 10.7, f 10.2) 
do  1®  1=1, 10® 
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read(67, 1 )PR( I ) , HEM ) 

10  continue 
20  close  (67)  •  ■ 

New2$= ’ rename  V/*  file  XYZ1  A IV/’  PRANTL 
♦OUTPUT  A 1  * 
call  excms(NEW2$) 

C  .  ..  ...  •r:,'..  .  .  . 

c  The  remainder  of  ther program  cal  Is  DISSPLA  subroutines 
c  to  create  and  plot  the  desired  curves.  Call  teic618 
c  sets  up  the  required  machine  for  the  plotting  and 
c  l^nks  you  to  the  plotter.  You  must  be  on  a  tek618 
c  terminal  to  use  this  program.  Call  blowup  is  used  to 
c  magnify  the  plot  to  near  full  page  si2e,  for  clarity, 
c  Call  page  sets  up  the  plage  size,  in  inches.  Call  area2d 
c  sets  up  a  two  di^  nsional  plot  area  of  the  specified 
c  size,  in  inches.  Call  xname  and  yname  allow  you  to 
c  label  the  axes.  Call  headin  allows  naming  the  plot, 
c  Call  graf  sets  up  the  physical  origin  and  the  maxima 
c  for  the  axes.  The  ’scale’  portion  of  this  call  allows 
c  DISSPLA  to  select  appropriate  increments  for  the  axes 
c  based  on  the  minima  and  maxima  specified.  Call  thkfrm 
c  and  call  frame  place  frames  around  the  inier  plot  area 
c  and  the  entire  plot,  of  a  specified  thickness.  Call 
c  curve  tells  the  plotter  to  plot  the  specified  arrays  of 
c  data  points.  The  arrays  must  be  called  in  the  proper 
c  order  (y-axis,  x-axis)  to  get  the  correct  curve.  Call 
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c  endpl  ends  the  specific  plot  end  Sets  up  a  new  page  for 
c  another  plot.  Call  donepl  tells  the  devide  you  are  done 

c  plotting.  . . 

c 

call  tek618 
call  blowup(  1.5) 

.•v  /  ■  ■  ...  _  ‘  ( ' 

<•  cal  1  page (1 1,8.5)  . ■  ■ 

°  'call  area2d(8. ,6. )  ..  .  •.• 

call  xname( ’helium  fractions*, 100) 
cal  1  yname(  ’ Prandtl  numbers  * ,  100)  •  /; 

call  headin ( ’Prandtl  number  vs.  helium  fractions*, 

1  *  *  f '  . .  .  •  • 

*100,  1. ,  1) 

call  graf (0.0, ’scale*, 1.0, 0.4, ’scale’, 0.7) 
call  thkfrm( .02) 
call  frame 

call  curve (HE, PR, 100,0) 

call  endpl (0) 

cal  1  donepl 

stop 

end 
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APPENDIX  C.  RAW  DATA 

1 '  This  Appendix  provides  a  listing  of  our  spreadsheet, 
discussed  in  Chapter  V,  including  all  recorded  data  and 
calculated  values  that  we  have  used.  Many  of  the  values  are 

t  * 

not  related  to  our  coefficient  of  performance  evaluation, 
but  Were  calculated  for  informational!  purposes. 

The  Appendix  is  broken  down  into  three  sets  of  three 
pages  each,  in  order  to  fit  the  24  column  by  122  row  spread¬ 
sheet  within  the  desired  margins.  Each  set  of  three  pages 
covers  the  entire  24  columns  for  the  specified  data  set3, 
with  nine  columns  per  page.  We  will  discuss  each  of  the 
columns,  omitting  the  first  which  specifies  the  data  set. 

Columns  1,2,3,4,6,8,9,10  and  11  are  all  raw  data  and  are 
self  explanatory.  Column  5  is  the  measured  mean  pressure  of 
column  4,  after  conversion  from  milli-volts  to  psia  (at  1.5 
psia  per  mv).  Column  7  converts  the  measured  dynamic 
pressure  signal  of  column  6  to  percent  pressure  amplitude 

( P>*/P  =1586xC6/C5) .  Column  12  is  an  acceleration  term  that 
0  m 

is  calculated  for  use  in  the  work  equation  as  a  correction 
for  accelerometer  non-linearities  (acc=0.2848xC3xC8/170.5) . 
Column  13  is  the  absolute  temperature  ratio 
(TC/TH=(C2+273. 15)/(Cl+273. 15) 1 .  Column  14  is  the  applied 
heater  load  in  watts  (QHEAT=C9xC 10x0. 001  or 

C 

QHEAT=C 10xC 10xC 1 1/ 10  ,  as  appropriate).  Column  15  is  the 
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work  in  watts  (H=598l0x( 1-7. 4xl0“5xC12)xC6x0.003xC8/C 1839.8- 
0.05(03-550 ) 3 ,  as  discussed  in  Chapter  V.  Column  16  is  the 
coefficient  of  performance,  neglecting  the  heat  leak  from 
the  cold  end  of  the  refrigerator  (C0P=C 14/C  15) .  Column  17 
is  the  coefficient  of  performance  relative  to  Carnot, 
neglecting  the  heat  leak  (C0PR=C16( ( 1/C 13)-1 ) ) .  Column  18 
is  the  heat  leak  in  watts,  as  discussed  in  Chapters  IV  and 
V  (HEATLK=0.0064(C1-C2) ) .  Column  19  is  the  coefficient  of 
performance,  taking  into  account  the  heat  leak 
(COP2=(C14+C18)/C15) ...  Column  20  is  the  coefficient  of 
performance  relative  to  Carnot,  taking  into  account  the 
heat  leak  (dOPR2=C 19( ( 1/C 13>-1 ) ) .  Column  21  is  the  total 
heat  load  in  watts  ( Qtot=C 14+C 18 ) .  Column  22  is  the  overall 
temperature  difference  (DELTA  T=C1-C2).  Column  23  is  an 
approximate  viscous  penetration  depth  for  the  given 
temperature,  pressure,  frequency  and  gas  type  and  an 
approximate  value  for  dynamic  viscosity  (VISC 
PEN=SQRT( (C 1+273. 15) (0.005014692852) /(C3xC5) ),  where  the 
coefficient  takes  into  account  all  constants  and  the  value 
used  for  dynamic  viscosity).  Column  24  is  an  approximate 

_  U 

thermal  penetration  depth  calculated  from  6  =6  <r  (THERM 
PEN=1, 219820636xC23) .  The  equations  for  viscous  and  thermal 
penetration  depths  use  a  single  value  for  and  v  within 
each  data  set,  respectively,  instead  of  a  value  calculated 
for  each  data  point.  This  is  done  for  ease  of  calculation. 
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The  resulting  values  are  therefore  not  precise  but  merely 
reasonable  approximations. 

The  spreadsheet  of  raw  data  and  calculated  values  is 
provided  on  the  next  nine  pages,  and  does  not  follow  on 
this  page  for  formatting  reasons. 
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TH 

TC 

FREQ 

PflEAN 

PtCAN 

VPRESS 

Po/Pe 

VLQCK1N 

VHEAT 

deg  C 

deg  C 

Hz 

•Vres 

psia 

Vres 

X 

aVres 

V 

25.1 

-58.6 

522.03 

100. 17 

150.26 

0.28391 

2.997 

367.3 

0.000 

22.4 

-57.2 

523.20 

99.98 

149.97 

0.18870 

1.996 

241.6 

0.000 

old 

22.6 

-44.7 

537.10 

100.00 

150.00 

0.18890 

1.997 

262.1 

6.870 

stack 

22.9 

-34.0 

548.70 

99.97 

149.96 

0.18860 

1.995 

277.0 

9.714 

pure  He 

23.2 

-24.7 

558.70 

99.98 

149.97 

0.18940 

2.003 

289.1 

11.897 

short 

20.3 

-13.5 

570.40 

101.06 

151.59 

0.18870 

1.974 

307.0 

13.780 

resonator 

23.1 

-13.4 

570.40 

100.40 

150.60 

0.18900 

1.990 

309.0 

13.780 

24.7 

-43.0 

539.60 

99.93 

149.90 

0.28280 

2.992 

391.5 

13.780 

22.8 

-62.2 

518. 10 

99.96 

149.94 

0.18960 

2.006 

228.2 

0.000 

24.2 

-71.3 

508.40 

100.05 

150,08 

0.28390 

3.000 

326.3 

0.000 

22.6 

-67.5 

507.94 

99.54 

149.31 

0.18902 

2.008 

233.1 

25.9 

-68.6 

506.90 

100.20 

150.30 

0.28332 

2.990 

349.8 

23.2 

-44.7 

533.70 

104.17 

156.26 

0.18873 

1.916 

256.1 

23.5 

-34.0 

544.90 

102.50 

153.75 

0.18903 

1.950 

280.2 

new 

23.6 

-24.9 

554.49 

103.20 

154.80 

0.18947 

1.941 

293.1 

stack 

24.1 

-13.5 

562.52 

102.00 

153.00 

0.18939 

1.963 

319.8 

pure  He 

26.9 

-43.1 

535.55 

100.10 

150.15 

0.28340 

2.993 

416.8 

short 

25.3 

-73.5 

501.57 

100.04 

150.06 

0.28367 

2.998 

337.2 

resonator 

22.6 

-68.6 

507.10 

99.97 

149.96 

0.18900 

1.999 

230.1 

20.7 

-34.0 

544.50 

100.09 

150.14 

0.09465 

1.000 

147.7 

23.7 

-24.5 

554.62 

100.05 

150.08 

0.18886 

1.996 

309.0 

23.7 

-13.4 

565.44 

100.11 

150.17 

0.18922 

1.998 

331.7 

23.0 

-44.6 

533.81 

100.07 

150.11 

0.18899 

1.997 

272.5 

23.7 

-6B.9 

235.50 

98.67 

148.01 

0.28390 

3.042 

169.3 

0.000 

24.1 

-60.5 

239.86 

98.13 

147.20 

0.23310 

3.050 

181.4 

6.871 

old 

24.2 

-51.9 

244,21 

97.84 

146.76 

0.28310 

3.059 

194.1 

9.706 

stack 

24.3 

-42.7 

248.82 

97.54 

146.31 

0.28380 

3.076 

208.6 

11.904 

He-Xe 

22.4 

-38.4 

250.88 

97.94 

146.91 

0.18870 

2.037 

142.2 

8.416 

short 

23.2 

-70.1 

234.76 

98.84 

148.26 

0.28330 

3.031 

167.9 

0.000 

resonator 

21.0 

-62.5 

238.85 

99.14 

148.71 

0.19000 

2.026 

118.6 

0.000 

22.2 

-17.5 

260.29 

98.60 

147.90 

0.18950 

2.032 

159.8 

11.898 

22.9 

-77.6 

229.24 

100.99 

151.49 

0.28380 

2.971 

160.3 

24  1 

-60.5 

070 
iw'Qt  iO 

101.12 

151.68 

0.28340 

2.963 

181.7 

new 

o?  q 

-51.9 

242.48 

101.03 

151.55 

0.28291 

2.961 

193.5 

stack 

24.5 

-42.8 

246.82 

101.73 

152.60 

0.28291 

2.940 

202.8 

He-Xe 

22.1 

-38.3 

248.93 

101.63 

152.45 

0.18910 

1.967 

141.5 

short 

00  7 

-50.0 

243.33 

101.71 

152.57 

0.18870 

1.962 

130.0 

resonator 

21.9 

-67.5 

234.57 

102.03 

153.05 

0.18910 

1.960 

113.2 

23.5 

-77.8 

229.07 

100. 16 

150.24 

0.28340 

2.992 

161.6 

23.3 

-80.3 

227.13 

89.60 

134.40 

0.28421 

3.354 

183.8 
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I  HEAT 

RHEAT 

ACC 

TC/TH 

(MAT 

WORK 

cop..  . 

.core- 

HEATLK 

fl)A 

ohms 

cgsE-3 

watt 

Mtt 

Q*hutar 

cl7/cirn 

watt 

0.00 

320.28 

0.7194 

ft  OftOft 

utUVJCw 

9.9215 

0.0000 

0.0000 

0.5692 

0.00 

211.14 

0.7307 

0,3000 

4.3736 

0.0000 

0.0000 

0,5413 

old 

72.81 

235.15 

0.7725 

0.5002 

4.7429 

0.1055 

0.0311 

0.4576 

stack 

102.82 

253.88 

0.8078 

0.9988 

4.9991 

0.1998 

0.0475 

0.3869 

pure  He 

125.82 

269.80 

0.3384 

1.4969 

5.2348 

0.2860 

0.0551 

0.3257 

short 

145.60 

292.51 

0.8848 

2.0064 

5.5386 

0.3628 

0.0472 

0.2298 

resonator 

145.60 

294.41 

0.9768 

2.0064 

5.5747 

0.3599 

0.0506 

0.2432 

146.00 

352.87 

0.7676 

2.0119 

10.5129 

0.1914 

0.0580 

0.4740 

0.00 

197.49 

0.7128 

0.0000 

4.1544 

0.0000 

0.0000 

0.5780 

0,00 

277.10 

3.6788 

0.0000 

8.8393 

0.0000 

0.0000 

0.6494 

0.00 

94.00 

197.77 

0.6954 

0.0000 

4.2636 

0.0000 

0,0000 

0.5676 

0.00 

94.00 

296. 18 

0.6840 

0.0000 

9.5190 

0.0000 

0.0000 

0.5954 

89.50 

94.40 

228.31 

0.7709 

0.7562 

4.6697 

0.1619 

0.0481 

0.4278 

116.90 

94.48 

255.04 

0.8062 

1.2911 

5.1085 

0.2527 

0.0608 

0.3623 

new 

1 wO • 20 

94.55 

271.47 

0.8366 

1.6775 

5.3509 

0  •  o  1  j5 

0.0612 

0.3056 

stack 

158.30 

94.64 

300.49 

0.8735 

2.3716 

5,8244 

0.4072 

0.0590 

0.2369 

pure  He 

153.00 

94,40 

372.86 

0.7667 

2.3566 

11.2885 

0.2088 

0.0635 

0.4410 

short 

0.00 

94.00 

282.51 

0.6690 

0.0000 

9.1956 

0.0000 

0.0000 

0.6224 

resonator 

0.00 

94.00 

194.91 

0.6916 

0.0000 

4.2091 

0.0000 

0.0000 

0.5746 

0.00 

94.48 

134.34 

0.8139 

0.0000 

1.3606 

B  mm 

II. 

0.0000 

0.3446 

138.00 

94.55 

286.27 

0.8376 

1.8006 

5.6168 

0.3206 

0.0621 

0.3037 

156.80 

94.64 

313.29 

0.8750 

2.3268 

6.0303 

0.3859 

0.0551 

0.2337 

99.50 

94.40 

242.98 

0.7717 

0.9346 

4.9701 

0.1880 

0.0556 

0.4259 

0.00 

66.60 

0.6881 

0.0000 

4.6249 

0.0000 

0.0000 

0.5926 

72.97 

72.68 

0.7154 

0.5014 

4.9399 

0.1015 

0.0404 

0.5414 

old 

102.94 

79.18 

0.7441 

0.9992 

5.2838 

0.1891 

0.0650 

0.4870 

stack 

126.09 

86.70 

0.7748 

1.5010 

5.6900 

0.2638 

0.0767 

0.4288 

He-Xe 

89.11 

59.59 

0.7943 

0.7499 

2.5844 

0.2902 

0.0752 

0.3991 

short 

0.00 

65.84 

0.6852 

0.0000 

4.5772 

0.0000 

0.0000 

0.5971 

resonator 

0.00 

47.32 

0.7161 

0.0000 

2.1716 

0.0000 

0.0000 

0.5344 

125.77 

69.48 

0.8656 

1.4964 

2.9152 

0.5133 

0.0797 

0.2541 

8.00 

94.00 

cl.38 

0.6605 

0. 3800 

4.4137 

3.0000 

0.8000 

0.6131 

95.70 

94.16 

0.7154 

0.8624 

4.9930 

0. 1727 

0.3687 

fl.Slol 

new 

119.33 

94.29 

78.37 

0.7448 

1.3420 

5. j06o 

0.2529 

0.0866 

0.4624 

stack 

141.10 

94.41 

33.61 

0.7739 

1.8796 

5.5598 

0.3381 

0.0983 

0.4105 

He-Xe 

96.80 

94.45 

58.84 

0.7954 

0.8850 

2.5979 

0.3407 

0.0876 

0.3684 

short 

73.50 

94.32 

52.84 

0.7553 

0.5095 

2.3824 

0.2139 

0.0693 

0.4410 

resonator 

0.33 

94.00 

44.35 

0.6970 

0.0000 

2.0797 

0.0000 

0, 0000 

0.5453 

0.00 

94.00 

61.83 

0.6585 

0.0000 

4.4430 

0. 0000 

0.0000 

0.6179 

8.30 

94.00 

69.73 

0.6505 

0.0000 

5. 0646 

0.0000 

0, 3000 

0.6320 
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cop; 

C0PR2 

Qtot 

DELTA  T 

Vise  PEN  THERM  PEN 

M  d9 

c20/carn 

watt 

dtq 

cat 

C* 

0.0574 

0.0224 

0.5692 

33.70 

0.008307 

8.010133 

0.1238 

0.0456 

0.5413 

79.60 

0.008268 

0.010086 

old 

0.2020 

0.0595 

0.9578 

67.30 

0.008162 

0.009957 

stack 

0.2772 

0.0659 

1.3357 

56.90 

0.008081 

0.009857 

pure  He 

0.3482 

0.0671 

1.3226 

47.90 

0.008012 

0.009773 

short 

0.4043 

0.0526 

2.2362 

33.30 

0.007848 

0.009573 

resonator 

0.4044 

0.0568 

2.2546 

36.50 

0.007911 

0.009651 

0.2365 

0.0716 

2.4858 

69.70 

0.008203 

0.010006 

0.1391 

0.0561 

0.5730 

35.00 

0.008315 

0.010143 

0,3735 

0.3348 

0.6494 

95.50 

0.003410 

0.010259 

£• 1331 

0.0585 

0.5676 

90.10 

0.006413 

0.010262 

0.0625 

0.0289 

0.5954 

94.50 

0.008440 

0.010296 

0.2535 

0.0754 

1.1839 

67.90 

0.008031 

0.009796 

0.3237 

0.0778 

1.6534 

5".  50 

0.008017 

0.009779 

new 

0.3706 

0.0724 

1.9831 

48.58 

0,307921 

0.009663 

stack 

0.4478 

8.0649 

2.6085 

37.60 

0.007917 

0.009658 

pure  He 

0.2478 

0.0754 

2.7976 

70.00 

0.008229 

0.010038 

short 

0.0677 

0.0335 

0.6224 

98.80 

0.008483 

0.010348 

resonator 

0. 1365 

0.0609 

0.5746 

91.20 

0.008402 

0.010249 

0.2533 

0.0579 

0.3446 

54.70 

0.038077 

0.009853 

0.3746 

0.0726 

2.1043 

48.20 

3.008045 

0.009814 

0.4246 

0.0606 

2.5606 

37.10 

0.007966 

0.009717 

0.2737 

0.0810 

1.3605 

67.60 

0.008190 

0.009991 

0.1281 

0.0581 

0.5926 

92.60 

0.006535 

0.012596 

0.2111 

0.0843 

1.0428 

84.60 

0.006498 

0.012524 

old 

0.2813 

0.0967 

1.4862 

76.10 

0.006450 

0.012432 

5tSCK 

0.3392 

0.0986 

1.9298 

67.00 

0.006401 

0.012338 

Hg-xe 

0.4407 

0.1141 

1.1391 

60.30 

0.006341 

0.012223 

short 

0, i -0J 

0.0599 

3.5971 

93.30 

0.006554 

0.012595 

resonator 

0.2461 

0.0975 

0.5344 

33.50 

0.006444 

0.''",421 

8.6305 

3.0932 

1.7505 

39.70 

3.306203 

0..  .  ;55 

3. 1389 

8.0-14 

e.cni 

168.50 

0.00o538 

C.612it'2 

j « J  6 1 

Ci.  1093 

1.3784 

34.00 

8. 806422 

■6.012373 

new 

3.3400 

0. 1 165 

1.8044 

/  0 . 80 

8.806367 

0.012272 

stack 

3.4H9 

0. 1203 

2.2902 

Ll  •  00 

0.006295 

0.012134 

ne-xa 

3.4825 

0.1241 

1.2535 

60.40 

0.886246 

3.012040 

short 

Q, 3990 

0.1293 

0.9526 

72.30 

0..  23o317 

0.012177 

senator 

ill  OL'T-! 

U.  4.Q4.L 

0.1140 

0.5453 

3°.  48 

0.006420 

0.312374 

0.1391 

u.3721 

0.617? 

101.33 

0.006575 

0.012672 

8. 1246 

0,0670 

8.6320 

103.00 

0.006978 

0.013451 
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TH 

TC 

deg  C 

deg  C 

23.7 

-72.5 

new 

24.4 

-60.4 

stack 

24.2 

-51.8 

He-Ar 

25.  0 

-42.9 

(22%) 

22.3 

-38.2 

short 

22.5 

-50.1 

resonator 

22.2 

-65.0 

24.  a 

/  £.«*. 

22.9 

-60,3 

23.9 

-44.3 

new 

23.9 

-34.0 

stack 

24.6 

-25.1 

pure  He 

24.9 

-13.7 

long 

22.9 

-61.2 

resonator 

27.5 

-43.1 

28.8 

-32.3 

28.8 

-32.3 

26.1 

-62.0 

24.4 

-64.0 

new 

24.7 

-60.3 

stack 

25.2 

-51.8 

He-ftr 

25.8 

-43.0 

( 19%) 

22.9 

-38.4 

long 

22,9 

-50.1 

resonator 

22.1 

-60.9 

24.4 

-64.9 

21.8 

-63.1 

new 

22.6 

-50.1 

stack 

??  A 

mi  T 

-38.2 

He-Xe 

23.2 

-26.2 

long 

23.0 

-14.0 

resonator 

23.8 

-2.3 

21.7 

-63.0 

23.0 

-41.0 

23.7 

-34.1 

dual 

23.5 

-24.9 

stack 

24.1 

-13.1 

pure  He 

24.2 

-2.7 

long 

27.9 

-13.3 

resonator 

27.2 

-24.1 

27.1 

-33.4 

26.0 

-38.6 

24.8 

-55.5 

FREQ 

PMEAN 

PMEAN 

Hi 

eVres 

psia 

291.33 

180.84 

151.26 

299.42 

100.49 

150.74 

305.07 

100.88 

151.32 

310.85 

101.84 

152.76 

313.68 

100.82 

151.23 

386.06 

100.25 

150.38 

296.28 

100.25 

150.38 

291.54 

100.25 

150.33 

470.43 

109.40 

150.60 

485.56 

99.98 

149.97 

495.98 

99.95 

149.93 

504.18 

99.80 

149.70 

514.63 

100.38 

150.57 

469.90 

100.10 

150.15 

488.19 

99.92 

149.88 

498.54 

100.00 

150.00 

498.52 

99.93 

149.90 

468.88 

99.98 

149.97 

284.28 

99.94 

149.91 

286.60 

99.82 

149.73 

291.71 

99.70 

149.55 

296.91 

99.64 

149.46 

299.55 

100.09 

150.14 

292.58 

99,99 

149.99 

286.15 

100.02 

150.03 

283.80 

99.98 

149,97 

207.92 

99.78 

149.67 

213.39 

99.94 

149.91 

218.27 

99.91 

149.87 

222.90 

99.83 

149.75 

227.38 

99.91 

149.87 

231.63 

100.01 

150.02 

207.98 

99.85 

149.78 

490.81 

100.01 

150.02 

497.57 

99.34 

149.76 

506.39 

99.83 

149.75 

517.38 

100.56 

150.84 

526.79 

100.35 

150.53 

517.66 

100,27 

150.41 

507.60 

100.16 

150.24 

498.35 

100.18 

150.27 

493.36 

100.20 

150.30 

472. 10 

59.68 

89.52 

VPRESS 

Po/Pe 

VL0CK1N 

VTM 

% 

tVrts 

0.28311 

2.968 

190.1 

0.28301 

2.973 

208.1 

0.28321 

2.968 

220.3 

0.28303 

2.939 

229.3 

0.18950 

1.987 

168.6 

0.18878 

1.991 

149.8 

0.18921 

1.996 

135.1 

0.28333 

2.988 

191.7 

0.18916 

1.992 

212.8 

0,18902 

1.999 

249.2 

0.18920 

2.001 

276.7 

0.13902 

2.003 

298.7 

0.18934 

1.994 

326.3 

0.18932 

2.000 

212.2 

0.28343 

2.999 

386.0 

0.28365 

2.999 

425.6 

0.28359 

3.001 

426.2 

0.28335 

2.997 

316.7 

0.28315 

2.996 

183.5 

0.28300 

2.998 

191.7 

0.28310 

3.002 

208.1 

0.28268 

3.000 

226.7 

0.18942 

2.001 

159.6 

0.18923 

2.001 

142.0 

0.18982 

2.007 

126.1 

0.28274 

2.990 

182.1 

0.18940 

2.007 

94.1 

0.18902 

2.000 

108.7 

0.18948 

2.005 

124.2 

0.18916 

2.003 

138.6 

0.18957 

2.006 

154.2 

0. 18923 

2.001 

167.8 

0. 18941 

2.006 

94.3 

0. 18892 

1.997 

198.8 

0. 13833 

1.995 

208.8 

3. 13978 

2.010 

226.1 

0. 18866 

1.984 

241.1 

0.18900 

1.991 

260.2 

0.28320 

2.986 

372.1 

0.28346 

2.992 

345.3 

0.28369 

2.994 

322.8 

0.28325 

2.989 

309.7 

0.17804 

3.154 

353.6 

VHEAT 

V 


15! 


I  HEAT 

RHEAT 

ACC 

TC/TH 

aA 

ohftt 

cgsE-3 

0.90 

94.00 

92.51 

0.6759 

new 

80.19 

94.16 

104.08 

0.7150 

stack 

107.70 

94.29 

112.26 

0.7444 

He-Ar 

128.40 

94.41 

119.32 

0.7723 

(227.) 

91.20 

94.45 

84.13 

0.7952 

short 

67.10 

94.32 

76.58 

0.7544 

resonator 

0.00 

94.00 

66.86 

0.7048 

0.00 

94.00 

93.35 

0.6763 

0.00 

94.16 

167.22 

0.7190 

123.00 

94.40 

202.12 

0.7687 

new 

159.00 

94.48 

229.24 

0.3051 

stack 

183.50 

94.55 

251.56 

0.8331 

pure  He 

212.60 

94.64 

280.50 

0.8705 

long 

0.00 

94.14 

166.56 

0.7159 

resonator 

190.50 

94.40 

314.77 

0.7652 

237.60 

94.49 

354.42 

0.7977 

233.10 

94.49 

354.90 

0.7977 

0.00 

94.12 

248.04 

0.7056 

0.00 

94.08 

87.14 

0.7029 

new 

70.80 

94.16 

91.77 

0.7146 

stack 

130.60 

94.29 

101.40 

0.7419 

He-Ar 

163.60 

94.40 

112.43 

0.7699 

(197.) 

121.50 

94.44 

79,86 

0.7929 

long 

31.70 

94.32 

69.40 

0.7534 

resonator 

0.00 

94.15 

60.27 

0.7189 

0.00 

94.06 

86.33 

0.6999 

0.00 

94.15 

32.68 

0.7122 

new 

83.40 

94.29 

38.75 

0.7542 

stack 

120.20 

94.42 

45.28 

0.7950 

He-ke 

146.40 

94.54 

51.60 

0.3333 

long 

171.70 

94.67 

58.57 

0.8751 

resonator 

190.60 

94.80 

64.92 

0.9121 

0.00 

94.15 

32.76 

0.7127 

0.00 

94.39 

162.98 

0. 7839 

72.73 

94.46 

173.54 

3.8053 

dual 

114.50 

94.56 

191.25 

0.8369 

stack 

143.33 

94.67 

i08. u6 

0.8749 

pure  He 

175.70 

94.80 

228.96 

0.9095 

long 

207.00 

94.67 

321.75 

0.8632 

resonator 

154.80 

94.57 

292.77 

0.3292 

35.70 

94. 4o 

268.71 

0.7985 

0.00 

94.41 

255.22 

0.-841 

0.00 

94.23 

278.84 

0.7305 

QHEAT 

WORK 

COP 

CQPR 

HEATLK 

watt 

watt 

Q«hwt«r 

cl7/carn 

watt 

0.0000 

5.2182 

a  attoa 

0.0000 

0.5588 

0.6041 

5.7066 

0.1059 

0.0422 

0.4918 

1.0937 

6.0*i  77 

8.1810 

0.0621 

0.4408 

1.5565 

6.?9,9 

0.2472 

0,0729 

0.3938 

0.7856 

2.9544 

0.2659 

0.0635 

0.3509 

0.4247 

2.7462 

0.1546 

0.0583 

0.4211 

0.0000 

2.4835 

0.0000 

0.8080 

8.5058 

0.0000 

5.2659 

0.0000 

0.0000 

0.5580 

0.0000 

3.9001 

8.0000 

0.0000 

0.5242 

1.3594 

4.5538 

8.2935 

0.0898 

0.4328 

2.3385 

5.0523 

0.4728 

0.1145 

0.3643 

3.1837 

5.4409 

0.5851 

0.1172 

0.3131 

4.2776 

5.9424 

0.7198 

0.1071 

0.2432 

0.0000 

3.8926 

0.0000 

0.0000 

0.5298 

3.4258 

10.4881 

0.3266 

8.1002 

0.4448 

5.3343 

11.5415 

0,4622 

0.1172 

0.3849 

5.3568 

11.5549 

0.4636 

0.1176 

0.3849 

0.0000 

8.6416 

0.0000 

0.0000 

0.5558 

0.0000 

5.0388 

0.0000 

0.0000 

0.5746 

0,4720 

5.2597 

0.0897 

0.0358 

0.5525 

1.6082 

5.7084 

0.2817 

0.0980 

0.5005 

2.6834 

6.2051 

0.4325 

0.1293 

0.4472 

1.3941 

2.9346 

0.4751 

0. 1241 

0.3985 

0.6296 

2.6099 

0.2412 

0.0789 

0.4745 

0.0000 

2.3261 

0.0000 

0.0000 

0.5395 

0.0000 

4.9933 

0. 0000 

0.0000 

0.5805 

0.0000 

1.7318 

0.0000 

0.0000 

0.5349 

0.6558 

1.9959 

0.3286 

0.1071 

0.4580 

1.3642 

2.2852 

0.5970 

0.1540 

0.3818 

2.0263 

2.5450 

0.7962 

0.1593 

0.3112 

2.7910 

2.8365 

0.9840 

0.1405 

0.2331 

3.4439 

3.0800 

1.1182 

0.1077 

0. 1644 

0.0000 

1.7356 

8.0200 

0.0000 

0.5336 

8.0000 

3. 6356 

0. 3080 

0.0808 

0. o7 i 0 

0.4992 

3.8053 

0.1312 

0.0317 

0.3410 

1.2397 

4. 1467 

3.2990 

0.0583 

0.2S56 

2.0821 

4.3914 

0.4741 

0.0673 

0.21*5 

2.9265 

4.7417 

0.6172 

8.0614 

0. 1587 

4.0565 

10.0871 

0.4021 

0.0633 

0.2431 

2.2662 

9.3872 

3.2414 

0. 0497 

0.3027 

0.6938 

8.7963 

0.0789 

0.0199 

0.3570 

0.0000 

8.4337 

8.0000 

0. 0800 

0.3811 

0.0000 

6.0381 

0.0000 

0.8080 

0.4738 

152 


C0P2 

C0PR2 

Qtot 

DELTA  T 

Vise  PEN  THERM  PEN 

w  Cl9 

c20/carn 

watt 

deg 

Cft 

Cft 

3.1869 

0.0513 

0.5588 

96.20 

0.007043 

0.010759 

new 

0. 1921 

0.8765 

1.0960 

84.80 

0,00696? 

0.010644 

stack 

0.2539 

0.0872 

1.5345 

76.00 

0.006887 

0.010521 

He-Ar 

8.3097 

3.8913 

1.9503 

67.90 

0.006799 

0.010337 

Q21) 

3.3847 

0.8990 

1.1365 

60.50 

0.006773 

0,010347 

short 

8.3K30 

3.1082 

0.3457 

72.60 

0.006877 

0.010507 

resonator 

0.2036 

0.3853 

0.5058 

87.20 

8.006936 

0.010673 

0. 1360 

3. 8507 

0.5530 

96.20 

0.007364 

0.010792 

v* 

0.0525 

0.5242 

33.20 

0.008709 

0.010623 

0.3936 

3.1134 

1.7922 

a3 .  /  8 

0.088604 

8.013496 

new 

0.5450 

0.1319 

2. 7533 

57.90 

0.008515 

0.010387 

stack 

fl.6427 

8.1238 

3.4968' 

49.70 

0.008462 

3.010322 

pure  He 

0.760S 

0.1132 

4.5208 

38.60 

8.008355 

0.010192 

long 

8.1361 

8.0540 

8.5298 

84.10 

0.00872? 

0.010645 

resonator 

0.3693 

8.1133 

3. 8736 

70.60 

0. 003636 

0.010534 

3. 4955 

3.1257 

5.7192 

61.10 

0.008561 

0.013442 

0,49a9 

3.12&1 

5.7417 

61.10 

0.008564 

0. 010446 

0.8642 

3.3263 

0.5550 

83.10 

0.008789 

6.010720 

8.1140 

8.8482 

0.5746 

88.40 

0.007475 

0.011314 

new 

0. 1943 

0.0773 

1.0245 

85.00 

0.007453 

0.011281 

stac/ 

0.3694 

3.1295 

2.1087 

77.00 

0.307398 

8.011198 

He-Ar 

0.5045 

8.1503 

o. 1 J06 

68.80 

0.007343 

8.011114 

;  i  av . 

\  4  11%} 

0.6109 

0.1595 

1.7926 

61.30 

0.007259 

0.81G926 

long 

0.4230 

Q. 1384 

1.1041 

73.00 

0.007348 

0.811122 

resonator 

0.2319 

0.0907 

0.5395 

33.00 

0.007419 

0.011229 

8.1162 

0.0498 

8.5805 

89.30 

0.007480 

0.011322 

0.3089 

0.1248 

0.5349 

84.90 

0.006S94 

0.0132B8 

new 

0.5581 

0.1319 

1.1138 

72.70 

3.006809 

0.013124 

stack 

0.7640 

0.1971 

1.7460 

60.60 

0.006731 

0.012974 

.  He-Ae 

0.9135 

0.1337 

2.3375 

.49.40 

0.006673 

0.012861 

long 

1.3661 

0. 1522 

3.0241 

37.00 

0.006602 

0.012724 

resonator 

1.1715 

0.1129 

3.6084 

26.10 

0.006546 

8.012618 

0.3075 

0.1239 

0.5336 

84.70 

0. 806890 

0.013279 

0. 1339 

0.0286 

0.377a 

64.00 

8.808544 

0.013422 

V «  i-*-  *-‘Q 

U.0534 

u.  Sh0.< 

57.30 

0. 006503 

8.0103:: 

Juai 

3. ~a7s 

'0.3/1  ■ 

*  C1C7 

i  ■ 

S3. 40 

0.88842s 

.  %.  i  7  7 

stack 

fci.ji.-tl 

0.0/01) 

2.3015 

37.20 

0.  003-314 

3.010142 

purs  He 

8.  oOfcj  •' 

0.0o47 

3.0852 

26.90 

8.008250 

0.810Q63 

.or.g 

0.4262 

8.8676 

4.2996 

41.20 

0.008377 

0.01821? 

resonator 

0.2737 

0.0564 

2.5689 

51.30 

0.308455 

0.010313 

Q. 1194 

8.0381 

1.0507 

60.50 

8.30S53G 

0.310406 

0.0452 

0.0124 

0,3311 

64.68 

0.808557 

0.010438 

0,0783 

0.8289 

0.4733 

80.30 

0.311312 

0.013798 
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TH 

TC 

FREQ 

PtCAN 

RICAN 

VPRESS 

Po/fta 

VL0CK1N 

deg  C 

deg  C 

Hz 

aVras 

put 

Vres 

X 

*Vre* 

23.7 

-41.8 

217.92 

188.12 

158.18 

0.28365 

2.996 

144.7 

dual 

24.6 

-31.2 

222.87 

188. 73 

151.18 

0.28388 

2.971 

154.9 

(tick 

24.8 

-21.1 

226.16 

108.63 

150.95 

8.28335 

2.977 

168.1 

He-Xe 

25.  6 

-11.4 

238.13 

181.49 

152.24 

8. 28289 

2.947 

177.9 

long 

22.2 

-48.6 

217.91 

97.26 

145.89 

8.18935 

2.858 

99.8 

resonator 

23.8 

-28.0 

223.30 

188.87 

151.31 

8.18898 

1.988 

184.0 

22.5 

-14.8 

228.81 

99.64 

149.46 

8.18928 

2.089 

118.6 

24.4 

-41.0 

217.90 

99.70 

149.55 

8.28332 

3.085 

144.7 

23.3 

-^o,6 

238.56 

99.70 

149.55 

8.28328 

3.004 

170.3 

dual 

24.0 

-27.9 

242.61 

108.30 

158.45 

8.28351 

2.989 

177.1 

stack 

24.2 

-16.9 

247.39 

108.56 

158.84 

8.28333 

2.979 

186.6 

He-Xe 

24.4 

-7.0 

251.73 

181.53 

152.30 

8.28349 

2.952 

193.8 

short 

22.0 

-34.6 

239.37 

99.43 

149.15 

0.18893 

2.809 

113.6 

resonator 

22.3 

-23.8 

244.41 

100.44 

158.66 

0.18864 

1.986 

118.5 

22.1 

-11.9 

249.70 

101.03 

151.55 

0.18987 

1.979 

125.9 

24.1 

-36.3 

238.77 

99,40 

149.10 

0.28244 

3.004 

170.6 

25.9 

-31.0 

545.83 

100.49 

150.74 

0.28347 

2.983 

350.1 

26.6 

-22.9 

554,24 

98.71 

148.07 

0.28264 

3.028 

374.1 

dual 

26.6 

-15.0 

562.48 

99.23 

148.85 

8.28305 

3.016 

386.1 

stack 

27.2 

-7.1 

570.22 

99.49 

149.24 

0.28270 

3.004 

400.1 

pure  He 

23.1 

-34.6 

541.73 

100.89 

150.14 

8.18935 

2.000 

224.7 

short 

23.4 

-22.9 

554.27 

101.09 

151.64 

0.18919 

1.979 

235.3 

resonator 

23.4 

-12.9 

564.27 

99,59 

149.39 

8.18939 

2.011 

254.4 

23.9 

-3.0 

574.07 

100.11 

150.17 

0.18897 

1.996 

264.1 

23.6 

5.0 

581.68 

99.83 

149.75 

0.18951 

2.007 

276.7 

23.4 

-32.2 

543.70 

99.83 

149.75 

U.lUWf 

2.003 

227.8 

154 


1HEAT 

RHEAT 

ACC 

TC/TH 

QHEAT 

WORK 

COP 

CQPR 

HEATLK 

•A 

OhM 

cg*E-3 

Hitt 

Matt 

0*  htat.tr 

cl?/carn 

Matt 

e.N 

94.38 

52.67 

0.7821 

8  nee 

3.9758 

0.9181 

9.0000 

8.3921 

dual 

iae.aa 

94.48 

57.46 

8.8126 

8.9448 

4.2465 

0.2225 

8.0513 

0.3381 

stack 

147.60 

94.40 

63.58 

8.3468 

2.  '8689 

4.6112 

8.4469 

0.0814 

3.273* 

He-X* 

179.58 

94. 7® 

68.39 

8.3762 

3.0513 

4.3788 

9.6264 

8,0335 

8.2242 

long 

8.90 

94.39 

36.84 

0.7874 

0.0000 

1.8181 

ft  ftftftft 
Wi,  OTPS 

ft  ftftftft 
lit  vwft 

0.3306 

rascnator 

77.10 

94.52 

38.79 

8.3278 

0.5619 

1.9052 

8.2949 

8.0613 

9.3091 

121.10 

94.08 

45.33 

8.8765 

1.3882 

2.1763 

0.6379 

8.8898 

9.2212 

8.99 

94.33 

52.6? 

9.7382 

8.0000 

3.9711 

0.0000 

8. 9880 

3.  ^  **i  v 

0.80 

94.43 

67.36 

9. 7?io 

a. sees 

4.6/04 

8.0080 

8.0000 

8. vOU0 

dual 

69.00 

94.52 

71.77 

0.3253 

0.4588 

4.3600 

9.0926 

0.0190 

siack 

106.40 

94.64 

77.11 

0.8618 

1.0714 

5.1161 

0.2094 

0.0336 

0.2459 

He-Xe 

130.30 

94.75 

31.49 

0.3943 

1.6087 

5.3154 

9.3026 

0.0357 

0.1871 

short 

0.08 

94.46 

45.42 

8.8082 

0.0000 

2.0813 

8.0000 

0.3373 

'asonator 

57.90 

94.57 

48.38 

8.3440 

0.3178 

2.1876 

0.1463 

0.0270 

0.2743 

33.70 

94.70 

52.51 

0.8348 

0.6634 

2.3873 

0,2375 

0.0374 

0. 2026 

0.93 

94.43 

63.04 

8.7963 

0.0000 

4.6648 

8.0000 

0.3600 

0.00 

94.49 

319.20 

0.0000 

9.5091 

0.0080 

0.0000 

0.3653 

36.50 

94.5? 

346.34 

9.8349 

0.7076 

10.1127 

Q.07N 

0.0138 

0.3178 

5ual 

119.70 

94. 6o 

362.71 

0.8612 

1.3563 

10.4416 

0.1299 

0.0209 

0.2671 

stack 

143.60 

94.74 

331.09 

0.8858 

1.9536 

10.7941 

0.1810 

0.2202 

pure  He 

8.00 

94.45 

203.33 

8.3052 

8.0000 

4.1120 

0.0000 

0.0000 

3. 3704 

short 

69.40 

94.5? 

217.35 

8.8439 

8.4555 

4.2992 

0.1059 

0.0196 

0.2972 

resonator 

96.68 

94.67 

■  239. 78 

0.3776 

0.3834 

4.6467 

0.1901 

0.0265 

0.2330 

115.00 

94. 3 

253.25 

0.9094 

1.2541 

4.3096 

8.2603 

8.0260 

a. 1727 

129.40 

94.90 

268.35 

0.9373 

1.5890 

5.0486 

0.3147 

8.0218 

0.1194 

0.00 

94.4? 

206.88 

0.8125 

0.0080 

4.1621 

0.0000 

0.0000 

0.3570 

155 


:op2 

C0PR2 

Qtot 

DELTA  T 

Vise  PEN  THERM  PEN 

M  Cl9 

c28/carn 

Matt 

MQ 

CM 

CM 

it.  9986 

8.8275 

8.3921 

64.78 

8.886744 

8.812999 

dual 

8.3821 

0.0697 

1.2829 

55.88 

8.886671 

8.812856 

stack 

8.3873 

6.8924 

2.3391 

45.98 

8.006616 

8.812752 

rt#-x* 

6.6725 

6.8951 

3.2755 

37.88 

8.686539 

8,812684 

Ion? 

6.2893 

8.6565 

6.3886 

62.88 

8.686826 

8.813156 

escnator 

6.4571 

8.8951 

8.8709 

51.68 

6.86m38 

8.812779 

8.7395 

8.1842 

l.o894 

36.58 

8.886584 

8.812691 

6.6993 

6.8281 

6.3963 

65.48 

3.L'367o? 

6.613843 

6.6771 

6.8i9~ 

6.3.068 

50.40 

0.86o4a; 

0.612453 

dual 

6. .562 

6,2331 

8.7593 

51.90 

8.086389 

6.012315 

stack 

8.2573 

8.8413 

1.31a4 

41.13 

8.886321 

8.012134 

He-X* 

6.3379 

8.8399 

1.7958 

31.48 

8.686239 

0.812825 

short 

8.1621 

8.8385 

8.3373 

56.68 

8.086439 

8.812411 

resonator 

0,2738 

0.8585 

8.5913 

46.18 

8.806343 

6.812226 

6.3753 

6.6488 

0.8661 

34.80 

0.006255 

6.612857 

0.8771 

0.8197 

6.3088 

60.48 

0.606471 

3.812472 

8.8334 

0.0898 

8.3653 

56.98 

8.89S122 

8.889987 

6. 1814 

8.8281 

1.8254 

49.58 

8.886142 

8.889932 

dual 

6.1555 

8.8251 

1.6234 

41.&8 

0.083862 

8.889834 

stack 

0.2814 

0.0260 

2. 1 >  a8 

34.38 

8.083804 

0.809763 

pure  He 

6.8961 

3.8213 

Ci 

Kt 

cz* 

4- 

57.70 

8.888131 

0.389918 

snort 

a. 1751 

8.0324 

0.7527 

46.38 

8.088882 

0.089761 

resonator 

0.2403 

8.8335 

1.1165 

36.30 

8.887991 

0.889747 

0,2967 

0.0295 

1.4268 

26.98 

8.887988 

0.889647 

3.3384 

8.8226 

1.7835 

13.68 

3.087363 

8.089592 

0.6358 

8.0198 

8.3570 

55.68 

6.088131 

3.089918 
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